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FLOOD CONTROL ASPECTS OF CAUCA VALLEY DEVELOPMENT 


Phillip Z. Kirpich,! F. ASCE and 
Carlos S. Ospina, 2 F. ASCE 


SYNOPSIS 


The Cauca Valley project applies two techniques relatively new in Latin 
America: integrated water-resource development and administration by a 
: decentralized, TVA-type agency. The flood-control aspects are greatly af- 
: fected by the other aspects, mainly power, irrigation and drainage. 
: Parallel matters discussed herein are the hydrology of the region, which 
is in the equatorial “doldrum belt,” and the rate of development as affected 
: by economic and political considerations. 


INTRODUCTION 


Flanked by the two westernmost ranges of the Colombian Andes, the Cauca 

Valley is a flat, verdant strip of land 200 km long, an average of 20 km wide 

| and almost 400,000 ha in extent; see Figs. 1 and 2. ( Abbreviations and equi- 

valents are listed at the end of this paper.) 

As the elevation is 950 m above sea level, it enjoys a pleasant, balmy 

Climate, despite the equatorial latitude. Long-term averages are: tempera- 

ture 24°C, humidity 63%, and annual rainfall 1,200 mm. 

Although well settled for over four centuries, only a third of the valley is 

| cultivated at present. This low degree of development results from a com- 
plex of sociologic and technologic reasons, some of which are discussed in 

this paper. 

| _ The Corporacidn Autédnoma Regional del Cauca (Autonomous Corporation 
of the Cauca Region) was established in 1954 by the government of Colombia. 


Jote: Discussion open until February 1, 1960. To extend the closing date one month, a 


written request must be filed with the Executive Secretary, ASCE. Paper 2147 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 9, September, 1959. 
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1. Associate, Tippetts-Abbett-McCarthy-Stratton and Chf. Engr., Corporacion 
Auténoma Regional del Cauca, Cali, Colombia. " 
2. Partner, Olarte, Ospina, Arias y Payan, Ltda., Bogota, Colombia. 
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David E. Lilienthal, former chairman of the TVA, assisted the government ir 
preparing the basic legislation; for this reason, the Cauca-Valley Develop- 
ment is often referred to locally as the “Plan Lilienthal.” The area of inter- 
est of the Corporation was defined as comprising the upper Cauca River 
drainage basin (above Cartago) and the adjacent watershed of the Pacific. 
Water-resource developments are planned within three departamentos ; 
(provinces): Cauca, Valle and Caldas: Hence, the abbreviation “CVC” in ~ 
common local usage. The CVC is charged with the integrated development 
of water resources including flood control, electric power, drainage, irriga- 
tion, pollution abatement and related uses. The CVC in 1955 engaged the 
consulting engineering firms of Olarte, Ospina, Arias & Payan, Ltda. of 
Bogota, Colombia, and Gibbs & Hill, Inc., and Tippetts-Abbett-McCarthy, bot! 
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of New York, to carry out the planning and engineering required for the Cauc: 
Valley development. This paper is based on the studies made by these three 
firms, together with the technical staff of the CVC. 

Within CVC’s development area reside 2,600,000 persons. The center of 
the valley is Cali, a booming city of 400,000, with an annual population in- 
crease of 8%. A hundred kilometers to the west is the Pacific deepwater port 
of Buenaventura. 

Almost a third of the valley is in the flood zone. The river overflows al- 
most every year, in Some years as many as three or four times. A major 
flood occurs once in 7 to 10 years when about 85,000 ha are inundated by the 
river (see Figs. 3 and 4), plus about 30,000 ha by tributary streams. 

Because the banks of the Cauca River are higher than the adjoining land on 
each side, damages from Cauca River floods extend well beyond the 85,000-h: 
area subject to direct flooding in two ways: (1) Floodwater trapped behind the 
river banks causes high groundwater levels over extensive adjacent areas (a 
condition that would be materially relieved if there were adequate drainage 
outlets, even in the absence of general flood control). (2) Owing to unfavorabl 
outlet conditions, flooding extends over a wide area at the mouths of many of 
the tributaries. 

Despite the large area flooded, actual damage is small because the fre- 
quent flooding has discouraged significant developments within the extensive 
low-lying areas along the river. The early Spanish settlers placed their 
towns on high ground for sanitary reasons. The present use of land in the 
flood zone and adjacent areas is limited to cattle grazing; when floods occur 
the cattle are simply moved to higher ground and brought back after the land 
dries out. Thus most of the benefits from flood protection are those that will 
result in the future from intensification in use of the land (for agriculture, in- 
dustrial sites and for urbanization) rather than from avoidance of damage to 
existing crops and improvements. 

Natural valley storage is high because of the large area flooded. During 
the February 1950 flood, 31,000 ha upstream of Cali were submerged beneath 
740 million m3 of water, the average depth being 2.4 m. The resulting 
natural regulation caused the peak discharge at Cali to be considerably less 
than at the head of the valley (La Balsa), as follows: : 


Drainage area Discharge 
La Balsa 5,480 km2 1,225 m3/s 
Cali 9,060 km2 1,000 m?/s 
It is evident that flood control in the upper portion (say above Cali), if 


limited to levees or to channel improvements, would cause substantially 
greater flooding downstream. 


Flood Characteristics 


Regional Meteorology 


Table 1 gives average monthly rainfalls at three representative stations. — 
Popayan and Piendamé are in the “upper basin” (that portion of the drainage 
area upstream of the Timba damsite at La Balsa), while La Manuelita is 
representative of the valley. The valley has two distinct dry seasons; the up: 
per basin has a single dry season extending from June to September. | 
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TABLE 1 


MEAN MONTHLY PRECIPITATION 


(millimeters) 
Station Popaydn Piendamé La Manuelita 
University 
Period of 1929-35 
record 1944-56 1932-55 1900-54 
Month Total % Total 4% Total 4% 
January 179 9.7 216 10.2 72 6.4 
February 172 9.3 221 10.4 79 7.1 
March 146 7.9 178 8,4 111 9.9 
April 145 7.8 194 9.1 145 12.9 
May 132 Tek 207 9.8 131 11.7 
July 49 ae f 39 1.8 37 Bed 
August 34 1.8 34 1.6 43 3.9 
September 98 5e3 103 4.9 70 6.3 
October 262 14.2 284 13.4 142 12.7 
November 314 16.9 313 14.8 128 11.4 
- December 236 12.7 230 10.8 91 8.1 
Annual 1,853 100.0 2,121 100.0 1,120 100.0 


- Storm rainfall in the upper Cauca basin is spotty and erratic. There are 
harp variations with respect to time and place even during short periods of 
ss than a day. The precipitation records show that heavy rainfalls at the 
arious stations seldom coincide. 
_ Another characteristic of storm rainfall in the upper Cauca basin is its 
»w depth as compared with the annual rainfall. Thus, for the 18 years of 
aily records at Popayan, the maximum 24-hr rainfall has never exceeded 

of the average annual rainfall. This low percentage can be contrasted with 
.e much higher percentages observed in other parts of the world, outside the 
ititudes of the equatorial trough (or “doldrums”)in which the upper Cauca 

sin is located. In these latitudes, winds are mild and tropical hurricanes 
re unknown. (Tropical hurricanes universally originate at about 10° to 15° 
outh latitude and move respectively northward and southward.) Moreover, 
ld fronts, such as those that cause major storms in the Canal Zone, do not 
ach eee latitudes as they never penetrate further south than about 8°N 
titude. : 
The available records for the upper basin show that depths of storm rain- 
l are surprisingly low and do not increase with frequency as much as might 
‘expected; see Table 2. The 10-year 24-hr rainfall at Popayan of 100 mm 
about the same as at Washington, D. C. while the 100-year 24-hr rainfall at 


ae 
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TABLE 2 


FREQUENCY OF STORM RAINFALL 


(millimeters) 


l-year 10-year 100-year 


24-hr_ point rainfalls 


At Popaydn 64 100 136 
At Piendamé 56 88 120 


Basin rainfall (5,480 
km< above La Balsa 


24-hr 28 40 56 
5-day 80 114 150 


Popayan (estimated from an 18-year record) is about 10% less than at 
Washington; the corresponding basin rainfalls are only about 40% of the poin 
rainfall at Popayan. Admittedly, the Popayan record is not long enough for < 
firm conclusion; on the other hand, that the infrequent, say 100-year, storm 
are not much worse than the 10-year storms is confirmed from the records 
of floods, described below, which indicate that since 1916 all major floods 
have been of about the same order of magnitude; nor is there evidence of an 
earlier and greater floods from the time of original settlement of the Cauca 
Valley in the sixteenth century. 

A long-term record (1900 to date) for La Manuelita indicates a rough cor 
relation between annual rainfall in the valley as against the upper basin; thi: 
record also indicates that there is a 6 to 7-year cycle of unusual wet and un. 
usually dry years. 

The apparent explanation for this cycle is a periodic shifting of the equa- 
torial trough of low-pressure (inter-tropical convergence or “ITC”) and of 
the accompanying pattern of atmospheric circulation off the Pacific coast of 
the northern part of South America. There is a parallel shifting of currents 
in the Pacific ocean off the coast of Colombia and Ecuador; this is the 
general area where the cold northern-flowing Humboldt current meets a 
warm current coming from Panama; both currents then move westward to 
the Galdpagos Islands. See Fig. 1. 


Flood Frequency 


On the basis of gage records (available only since September 1945), the — 
frequency of the maximum recorded flood, that of 11 February 1950, is abou 
once in 10 years. Preceding the bewinnine of gage records, there were four 
“historic floods”: in 1916, 1932, 1934 and 1938. The available evidence ind 
cates that at least three of these were about equal (in peak discharge) to the 
flood of February 1950, confirming that a flood of the magnitude of the 
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February 1950 has a frequency of about 10 years, at least with respect to 
peak discharge. 

The worst flood of record at La Balsa had a discharge of 1,225 m/s from 
a drainage area of 5,480 km“. The corresponding Myers rating C, in the 
Jarvis-Myers formula Q=CVA , with Q and A in English units (cubic feet 
er second and square miles), is 1,050. The available records indicate that 
this value applies to the tributary streams as well. 

Compared with values observed in the United States, C is low in the Cauca 
Valley region; the apparent reason is the absence of both tropical hurricanes 
and storms involving cold fronts. Moreover, as precipitation in the form of 
snow is practically non-existent in the Cauca drainage basin, floods produced 
9y a combination of rain and snow-melt are unknown. 


Degree of Protection 


When the protected area consists of cities or towns, a high degree of pro- 
ection is justified to avoid catastrophic losses that may result from great 
loods occurring once in say 100 to 200 years. In the case of the Cauca Val- 
ey, however, and with the exception of the Aguablanca project described be- 
Ow, only agricultural land is involved. Following protection, the land in the 
lood zone will be devoted to crops (such as rice, corn, beans and cotton) and 
o intensive pasture, for which the average growing season (including harvest- 
ng) is only 5 to 6 months. 

Sugarcane, for which the growing season is 15 to 18 months, is widely 
rown in the Cauca Valley; harvesting takes place the year round—an impor- 
ant advantage in that the refineries can operate almost continuously—and its 
ultivation will surely be extended. However, it is estimated that not more 
nan 10 to 20 percent of the flood zone will be devoted to cane in the future: 
.rst, because the heavier soils located therein are better suited to other 
rops, mainly rice for which the high irrigation water-requirement will be 
elatively easy to supply in the flood zone; second, ample area, mostly free 
f flood risk, is available for cane cultivation outside the flood zone. 


_ Thus, for the expected cropping pattern, 6 months may be assumed to be 
ne average growing season. If protection is provided against a flood with 
verage frequency of once in 10 years (once in 20 growing seasons), 19 crops 
ut of 20 will be protected or 95% of the crops; if against a 20-year flood, 39 
ut of 40 or 97.5% of the crops. Thus, for a flood of half the frequency (with 
u€ Corresponding expense of protection on the order of 50% greater), the 
snefit gain is only 2.5%. 

_ The foregoing assumes that the area and consequently the crops damaged 
7 the 20-year flood, would be about the same as by the 10-year flood. This 
generally true for the largest part of the flooded area, located upstream 
 Buga. 

| Following the reasoning presented above, protection against a 5-year flood 
uld constitute 90% protection so that the benefit gain of 10-year over 5-year 
otection is theoretically only 5%; however, it is necessary to consider what 
ect the degree of risk will have on the attitude of farmers. Conversations 
ith them have indicated that 10-year protection corresponding to 5% risk 

le in twenty) is sufficiently great so that they will not be discouraged from 
inting crops, but that 5-year protection providing 10% risk (one in ten) 

ald have a discouraging effect. It has been concluded, considering the 


10 September, 1959 Hy § 


various factors, that the desirable degree of protection of agricultural lands 
is that obtained by providing control of the 10-year flood. 

There ‘are no important towns or villages within the Cauca River flood zon 
Some towns do suffer floods from tributary streams but these involve local 
projects generally requiring only minor works to achieve protection. 

The city of Cali, although not suffering serious flood damage at present, 
has expanded into the flood plain to the extent that further expansion has been 
limited to narrow flood-free zones to the north and south; see Fig. 9. The 
Aguablanca project, described below, is being constructed to protect 5,000 ha 
of land adjacent to the city against a flood having an estimated frequency of 
100 years, under present conditions. Following completion of the Timba 
reservoir and the Cauca River channel improvement, the degree of protectior 
will be against a flood of at least 200-year frequency. 


Methods of Obtaining Flood Control 


Upstream of the upper end of the valley at La Balsa, there are two reser- 
voir sites; Timba and Salvajina. The drainage areas controlled are, respec. 
tively, 5,480 km? and 3,960 km2, Other significant drainage areas are: 


At Juanchito bridge near Cali 9,060 km? 
At Buga 12,930 °° 
At Cartago 17,820 ™ 


The Timba reservoir is more valuable as a flood-control structure be- 
cause it controls 38% more drainage area than Salvajina. Moreover, various 
additional factors, which cannot be discussed here because of space limitatio 
have required that, of the two reservoirs, Timba be given a higher priority. 

The control of floods on the tributaries must be considered jointly with th 
drainage problem, as discussed in a later section. 


Optimum Combination 


The desired degree of flood protection is obtainable by storage in the Tim 
reservoir and improvement of the channel of the Cauca River, the optimum 
combination being that giving the lowest overall cost, considering the value o 
salable electric energy lost because of storage reserved for flood control. 

The channel capacity of the river is about 550 m3/s at La Balsa and about 
650 m3/ s at Cali. The amount of reservoir storage necessary to reduce pas 
floods of record down to these limits was determined, and this study resultec 
in the development of a “rule curve,” similar to that shown in Fig. 5, giving 
the reservoir level which may not be exceeded on a given calendar date. 
Other rule curves were developed for lesser degrees of flood reduction. 

At the time that Timba is constructed, there will already be in existence 
two major hydroelectric stations (Anchicaya and Calima—See Fig. 2) and one 
coal-burning thermal station with an installed capacity of 53,000 kw located. 
at Yumbo, a suburb of Cali. From a study of system power -output, using as 
basic data streamflow records from 1945 to date, it was determined that the . 
firm power obtainable from Timba, if the rule curve in Fig. 5 is followed, 
will be 60,000 kw. The average annual generation will be 346,000,000 kwh. 

Besides flood control and power, Timba will provide the means of divert- 
ing water to two principal irrigation canals, one on each side of the valley. 


HY 9 CAUCA VALLEY 11 


to] 

1} 

te 
” 
3 
oa 
=] 

= ¢ 
| 

T 


{ 


Pool elevation-m 


+ 


1 


November, 15 


Months 


OPERATING RULES 
Reservoir level Action required 


1. Higher than Rule Curve Moke the total outflow equal to 
but lower than El. 1025.0 the lowest of the following two 
quantities: 


(a) 790 m/s 
(b) 870 m¥g minus the 
tributary flow between 
: Lo Balsa and Cali 
2. Higher than El. 1025.0 Open all outlets to prevent a 
| further rise in reservoir level 


FIGURE 5. RULE CURVE, TIMBA RESERVOIR 


he total net irrigable area in the valley amounts to 280,000 ha, of which 
3,000 will be irrigated by water diverted at the Timba damsite. The un- 
gulated flow is sufficient at all times to meet irrigation water requirements. 
‘Taking the power output corresponding to the rule-curve operation as a 
rm, the optimum combination of Timba reservoir storage and the channel 
provement will be that giving the lowest overall cost of obtaining the fol- 
wing objectives: 


_ 10-year flood protection 

_ Irrigation diversion 

_ Firm capacity of 60,000 kw 

Average annual generation of 346,000, 000 kwh 


; 
L 
4 
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Table 3 presents cost estimates for the following four cases: 


1. No reservoir storage. For this case a low diversion dam for irriga 
tion would be constructed at an estimated cost of 22,000,000 pesos. Be- 
cause of the low head available and poor streamflow regulation, it would 
not be economical to construct a hydro station. Instead, a thermal station 
of 60,000 kw capacity would have to be constructed with costs as indicated 
in Table 3. 


2. Timba reservoir operated for maximum power output with flood con 
trol incidental. The firm capacity being 67,000 kw or 7,000 kw higher tha 
the norm, a credit of Ps 11 million is given in line 11 representing the co 
of 7,000 kw of thermal capacity; in this case, the Ps 11 million actually 
represents cost savings at the already-existing Yumbo thermal station. 


3. Timba reservoir operated according to the rule curve in Fig. 5. Th 
case is for the “norm” described above. 


4. Timba reservoir operated for maximum flood-control effect with 
firm power incidental. With the firm hydro capacity at 43,000 kw, 17,000 
kw of thermal capacity is required to meet the norm. From line 11, 
Table 3, this would cost Ps 26 million. As in Case 2, this cost is taken 
as the incremental cost for 17,000 kw at the pre-existing Yumbo thermal 
station. 


Peak discharges at La Balsa and Juanchito for these four cases would be 
as given in Table 4. Figure 7 presents the results of detailed computations « 
the effect of the Timba reservoir on the 1949-50 flood (Fig. 6) if operated ac 
cording to the rule curve in Figure 5. 

The capital cost of the channel improvements required for each of the fou 
cases, with design capacities at La Balsa and at Cali as given, respectively, 
in the second and fourth columns of Table 4, are given in line 8 of Table 3. 

The total capital cost (lines 12 or 23, Table 3) is increased in each case 
by the capitalized value of the annual cost of operation and maintenance, 
depreciation and fuel consumption. The overall cost thus obtained provides 
the basis for comparison. It is the lowest in Case 3, which therefore repre- 
sents the optimum solution. 


Flood Control and Drainage 


Groundwater levels are high during the wet season (and for a long period 
thereafter) within the 85,000-ha area subject to direct flooding from the 
Cauca River and in an additional area of about equal size immediately ad- 
jacent to the flood zone on each side; flood control of the Cauca River will _ 
thus benefit directly about 170,000 ha or over 40% of the valley. Drainage 
works, principally canals and pumping stations, will of course be necessary 
to supplement the main flood-control works (Timba reservoir and Cauca 
River improvement). | 

Flood control and supplementary drainage works within this 170,000-ha — 
low zone will also benefit the remaining 230,000 ha in the valley but in an in 
direct way. Irrigation of these 230,000 ha is proceeding on a gradual basis. 
(The rate will be accelerated once the Timba dam and the main irrigation — 
canals leading from it are completed.) A deterrent to more rapid develop- 
ment of irrigation has been lack of drainage outlets and high groundwater 


1. Flood control storage 
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TABLE 3 


COST OF FOUR ALTERNATIVE SOLUTIONS 


Case l. Case 2. Case 3. 
Storage Storage Storage 
negli- used for used 
gible maximum according 
power to rule 
output curve 


(million m3) 0 73(a) — 378(b) 
POWER OUTPUTS 
Firm Capacity (kw) 
2. Timba hydro 0 67,000 60,000 
3. Thermal 60,000 (-) 7,000 = 
4. Total 60,000 60,000 60,000 
Firm Generation (million kwh/yr) 
5. Timba hydro 0 381 346 
6. Thermal 346 «=(-) 35 -- 
7. Total 346 346 346 
CAPITAL COSTS (million pesos) (d) 
8. Channel improvements 283 x7 70 
9. Diversion dam for irrigation 22 _ - 
10. Timba multipurpose project - 203 198 
1l. Thermal station 92 (-) 11 - 
| 2. Total 397 363 268 
ANNUAL COSTS (thousand pesos) 
| Operation & Maintenance 
13. Channel improvements 6,000 3,500 2,000 
14. Diversion dam 90 -- ae 
15. Timba multipurpose project — 590 560 
16. Thermal station 1,300 (-) 150 = 
17. Total 7,390 3,940 2,560 
Depreciation 
18. Diversion dam 145 = am 
19. Timba multipurpose project = 1,340 1,310 
20. Thermal station 2,210 (-) 260 -- 
" Total 2,355 1,080 I, 310 
22. Fuel (f) 5,200 (-) 520 — 
OTAL COSTS (million pesos) 
3. Capital costs (line 12) 397 363 268 
4. 0 &M capitalized (g) 92 48 32 
5. Depreciation capitalized (g) 37 17 21 
6. Fuel capitalized (g) 65 (-) 7 - 
591 421 321 


Te Total 


t 
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maximum 
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effect 
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TABLE 3 (Continued) 


(a) Surcharge storage. 
(b) Includes 73 million m> surcharge storage. 
(c) Reservoir drawn down to El 1000. 


(d) Include interest during construction. Foreign-currency cost 
converted to pesos at rate of 1 US dollar = 7 pesos. 


(e) Depreciation assumed negligible for the channel improvements 
Depreciation computed on sinking-fund basis with interest 
at 4%. Resulting annual rates are? 


Life in Annual 
years rate 
Diversion dam and Timba 
multipurpose project 50 0.66% 
Thermal station 25 2.40% 


(f) Coal consumption 0.5 kg per kwh; coal price 30 pesos per 
metric ton; results in cost of 1.5 centavos per kwh. 


(g) The CVC can borrow local funds at an interest rate of 8% and 
foreign funds at 53%. Operation and maintenance, and fuel 
costs are in local currency and are therefore capitalized 
at 8% As 2/3 of the cost of the diversion dam and of the 
Timba multipurpose project are in foreign currency the 
weighted average interest rate is 6.33% For thermal 
stations 4/5 of the cost is in foreign currency, giving a 
weighted average interest rate of 6.20%. As these rates 


are close to each other, for simplicity a uniform rate of 
6.3% was used. 
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TABLE 4 


FLOOD REDUCTIONS BY TIMBA RESERVOIR 


(Cubic meters per second) 


La Balsa Cali 


No channel With channel 
improvement improvement 
1949-50 Flood (Desi Flood 


Case 1 1,300 1,000 1,500 
— 3x 950 1,200 1,210 
ase 790 870 870 
Case 4 620 700 700 


Average Flood 


Case 1 860 870 1,050 
Case 2 460. 670 670 
Case 3 660 870 870 
Case 4 300 510 510 


(x) Reservoir operated according to rule curve in Fig.5. For reservoir 
levels at or above the rule curve, release from reservoir is such 
that discharge at Cali equals 870 m?/s (see Fig.7). 


(xx) Controlled discharges for the average flood are higher for Case 3 
than for Case 2 because of rule-curve operation whereby reservoir 
is evacuated in anticipation of possible occurrence of the design 
flood (that of i eaeoeaye but note that Case 3 provides greater 
reduction for the design flood. 


levels in adjacent downslope areas. In many cases, a rather anarchic situa- 
tion has developed in which upslope landowners have thrown their irrigation 
wastes on their downslope neighbors; or, drainage-ways have been con- 
structed in the higher areas outletting to inadequate natural channels or per- 
manent swamps in the lower areas. 

It was stated above that drainage works alone would materially relieve 
damage caused by Cauca River floods. The Cauca River rarely remains 
above flood stage more than two weeks (many crops, if sufficiently mature, 
can withstand flooding for this long) but the low areas behind the riverbanks 
are now flooded for considerably longer periods, owing to the absence of 
drainage outlets. However, even if drainage outlets are constructed as a first 
stage (providing sizable benefits), general flood control will nevertheless 
soon become necessary to protect the high-value crops that will be grown, as 
well as the drainage facilities themselves. Moreover, most farmers will not 
risk planting of valuable crops in the absence of a sufficient degree of control 
of main-river floods. 

With respect to the tributaries the situation is different because the floods 
are of short duration with little volume in the peak. Studies of several typi- 
cal areas have shown that the most economical design is to provide a rela- 
tively low degree of protection against the tributary floods with occasional 


_ overflows (once in 4 or 5 years) removed by a system of drainage canals. A 


typical case (the Roldanillo-La Union-Toro project) is described below. 


tT 


18 September, 1959 HY 9 
The Timba Project 


The Timba dam will feature an earth embankment 46 m high and 1,200 m 
long; three spillway gates, each 13 m wide by 11.5 m high; a 60,000-kw 
powerhouse; and miscellaneous facilities including a flood-control outlet and 
irrigation intakes. The discussion herein will be limited to the spillway 
capacity and the method used for determining the economic justification of the 
project. 


Spillway Design 


A difficult problem has been establishment of a design flood for the spill- 
way. On the one hand, it is known that the upper Cauca basin is in a region of 
mild climate, completely free of tropical hurricanes and storms involving 
cold polar air masses. On the other hand, the period of actual streamflow 
records is short even though the qualitative historic record since the six- 
teenth century indicates the absence of floods of the order of magnitude re- 
corded in the United States. The hydrologic problem is further complicated 
in that it is practically impossible, owing to the paucity of data and the spotty 
and erratic character of storm rainfall, to apply methods now having standard 
acceptance involving “storm transposition” and the “unit hydrograph.” Other 
methods were therefore applied described briefly as follows: 


1. Correlation between basin rainfall and average daily discharge were 
attempted. The basin rainfall was determined by the Thiessen method; a 
more accurate method involving the plotting of isohyetal maps was not 
possible owing to lack of sufficient data. A correlation between the aver- 
age daily flood flow at La Balsa and the daily basin rainfall was found to 
be unsatisfactory; however, a better correlation was found between the 
total basin rainfall for the preceding five consecutive days and the aver- 
age daily flood flow. Considering all recorded daily flows higher than 300 
m3/s, the correlation plotted as shown in Fig. 8. It is thought that the en- 
velope curve represents the maximum daily discharge that can be reason- 
ably expected as a function of the 5-day rainfall. 


2. The 1000-year 5-day precipitation was obtained by extrapolating the 
18-year rainfall record up to a frequency of 1,000 years giving a 5-day 
basin precipitation of 186 mm (7.3 inches), or 62% more than the 10-year 
storm (see Table 2). It is believed that the upper Cauca basin, because of 
its location in the equatorial latitudes, can receive a probable maximum 
depth of rainfall of no greater order of magnitude than the 1000-year rain- 
fall. Such a premise is considered sound based on the following reasoning: 
great floods in the United States have been caused by rare and fortuitous 
“meteorological accidents” in which, however, one or more of the follow- 
ing conditions are always involved: a tropical hurricane, a meeting be- 
tween a cold polar air mass and a warm tropical air mass, Or snow-melt. 
None of these conditions are possible in the Cauca Valley region. 


3. The maximum average daily discharge for the spillway design flood ~ 
was determined from the envelope curve in Fig. 8. 


4. From study of the shapes of hydrographs for past isolated floods, 
which indicate the expected rate of rise and fall of floods, a hydrograph 


was sketched based on the maximum daily average discharge determined 
in Step 3. 
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FIGURE 8. DAILY DISCHARGE VS. 5-DAY PRECIPITATION 
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5. The resulting hydrograph had the following characteristics: 


Peak 6,000 m3/s 

Width of base 3 days 

Base flow 400 m3/s to 480 m3/s 
Volume (in 3 days) 685,000,000 m3 or 125 mm 


Comparison of the spillway design flood with past floods of record indi- 
cates that with respect to peak, the maximum flood of record has a peak dis- 
charge of 1,225 m3/ s; thus the spillway design flood is almost five times as 
great. With respect to volume, the following comparison is made: 


Maximum flood 


of record Spillway design flood 
Maximum 1-day volume 19 mm 65 mm 
vy 2-day * 36 mm 109 mm 
Wy) 3-day "' 53 mm 125 mm 


The peak discharge of 6,000 m3/s corresponds to a Myers rating C of 
4,550 in the formula Q=CVA_ where Q is the discharge in cubic feet per 
second, and A is the drainage area in square miles. 

With the inflow flood established and with the top of dam set at El 1029.0 
(based on topographic considerations), the capacity and dimensions of the 
gated spillway were determined based on the following criteria: 


Condition 1 Condition 2 
Reservoir elevation at start 
of flood 1025.0 1025.0 
Turbines functioning No No 
Outlet functioning Fully One of the two 
intake gates 
closed 
Freeboard 2.5 m Den 


Study of winds (mild in the upper Cauca basin) and reservoir fetch indi- 
cated that a freeboard of 1.5 m is ample. However, owing to the fact that 
Timba is an earth dam and that its safety depends on proper operation of the 
spillway and outlet gates, it was decided to allow an extra meter of free- 
board in Condition 1. In Condition 2, it was assumed that one of the two 
operating gates controlling the outlet conduit would be stuck in a closed posi- 
tion; however, as this is only a remote possibility (the gates will be operated 
frequently) a reduced freeboard of 1.5 m was accepted for this case. 

Conservatism is also involved in the adoption of an initial reservoir eleva- 
tion of 1025.0, which is the level of the top of the spillway crest gates. The 
record of greatest floods, Table 5, indicates that floods of appreciable volume 
(more than 10% of the flood storage of 305 million m3) have never succeeded _ 
each other within a short time; nor have such floods occurred from early 
March to mid-November. Thus, provided that the rule-curve operation is 
adhered to, the spillway design flood is not likely to occur on top of a full 
reservoir. On the other hand, it was felt that pressure from farmers along 
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TABLE 5 


FLOOD VOLUMES ABOVE RESERVOIR RELEASE 
Cauca River at La Balsa 
(Drainage Area 5,480 ba) 
Streamflow records, October 1945 to April 1958 


Volume Above 


Peak Release of 

Discharge Duration 650 m3/s 
Date m3 /s days Mm3 
May 3, 1949 1,225 3/4 20.1 
Nov. 22, 1949 1,075 6 1/2 49.0 
Dec. 16, 1949 1,000 3 45.1 
Feb. 10, 1950 1,225 12 243.7 
Nov. 27, 1953 885 4 32.0 
Dec. 1, 1954 931 5 1/2 67.6 
Jan. 5, 1955 950 23/2 36.0 
Jan. 1, 1956 1,225 9 1/2 132.5 


Notes: l. Only floods exceeding 15 Mm? above base of 
650 m3/s are listed. 


2, Reservoir release of 650 a°/s would be about 
the average maintained while storing the 
February 1950 flood. See Fig.7. 


the river could, in the future, upset the operating rules; they might insist 

that reservoir releases be kept lower than the amounts planned, even though 
the resulting damage would be minor. 

The adopted criteria resulted in a spillway consisting of three radial gates 

-13-m wide by 11.5-m high. 


| Cauca River Improvement 


i The two principal elements of a balanced flood-control plan are the Timba 
reservoir combined with the Cauca River improvement (Case 3 described 
above). Planning and design are well advanced for the Timba reservoir but 
not for the Cauca River improvement. Preliminary design of the latter has, 
however, been completed as summarized in Table 6; in general, design 
capacities are 20% to 25% higher than existing channel capacities. Admit- 
tedly approximate, the preliminary design is still considered good enough for 
the making of economic comparisons, which are based on estimated annual 
costs of alternative long-range plans. ; 

Final design of the Cauca River improvement (in common with all major 
river improvements) will be a highly complicated matter involving the follow- 
ing factors (among others): 


} § 
X 
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TABLE 6 


PROPOSED CAUCA RIVER CHANNEL IMPROVEMENTS 


(kilometers) 


La Balsa Buga to 


Total 
to Buga Cartago 

Length of Cutoffs 15 23 38 

Present River Length 179 194 373 

Future River Length 119 139 258 
River Shortening Due 

to Cutoffs 60 55 115 
Length of Right—bank 

Levees 26 37 63 
Length of Left—bank 

Levees 18 17 35 


1. Within the 400-km reach from La Balsa to Cartago, the Cauca is an 
alluvial river (as defined by Matthes. (3a)), The profile is parabolic; 
slopes are not great, varying from 40 to 15 cm per kilometer. The 
meanders are most prevalent in the upper 100 km where the slope is 
greatest; this is apparently explained by a higher proportion of bedload 
and higher erodibility of the riverbanks in this portion. The riverbanks, 
in general, consist of fine-grained silty sand having relatively low re- 
sistance to erosion; the river bed consists of sand and gravel. 


2. Since records began at Cali (Juanchito Bridge) in 1945, the river has 
been degrading at an average rate of 5.8 cm per year. (This degradation 
may be due, in part, to dredging of the river for sand and gravel to supply 
the local construction industry which has been increasingly active in re- 
cent years. It is estimated that 350,000 m3 are removed annually by hand, 
and about 100,000 m3 by mechanical means, within a reach of river 
several hundred meters long just upstream of the Juanchito bridge near 
Cali.) Such degradation has not, however, been noticeable at points along 
the river other than Cali. 


3. The Timba reservoir will have the effect of removing a major part 
of the sediment (suspended as well as bedload) now transported by the 
river, thereby increasing its erosive power; on the other hand, the re- 
duction of flood peaks caused by Timba will have an opposite effect. . 


4. Sediments will continue to be carried to the main river by the tribu- 
taries. However, some of the largest tributaries, such as the Palo, Fraile, 
Bolo and Amaime, traverse the valley floor for a distance of about 25 km; 
typically, these tributaries overflow and deposit the bulk of their sediment 
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in alluvial fans on the valley floor before reaching the river. 


5. The tributaries vary in their sediment-carrying characteristics ow- 
ing largely to differing soil characteristics and vegetal cover. For ex- 
ample, the Rio Bugala-grande is a notorious silt carrier; a soil- 
conservation program within its watershed would produce benefits to local 
agriculture and, incidentally, would greatly reduce the quantity of sediment 
now reaching the Cauca River from this source. 


6. At present, the river overflows frequently and deposits appreciable 
quantities of sediment on the flood plain. In the future, after construction 
of Timba and the channel improvement, the net sediment-load remaining 
will be confined to the channel instead of being deposited on the flood plain. 


7. Below Cartago, the river changes from an alluvial to an erodible-bed 
stream, flowing for the greatest part of its length through relatively nar- 
row canyons to its mouth at the Magdalena River. There is practically no 
settlement along the river below Cartago; therefore, any treatment of the 
river upstream can have no deleterious effects. 


8. Excessive degrading of the river in the reach from Cali to Yumbo 
would create serious dislocations because of the existence of several im- 
portant water intakes of which the most important are those for the munici- 
pal water supply of Cali, and for several important industries near Yumbo 
including a large steam-electric station. 


9. A vast hydroelectric project (capacity about 1,500,000 kw), for which 
a reconnaissance study has been completed, involves diverting about 130 
m3/s continuously from the Cauca River to the Pacific slope by means of 
a tunnel through the Cordillera Occidental. This diversion would require 
a barrier dam near Vijes. Such a barrier could be designed to prevent 
the degradation mentioned above. 


10. Some channel degradation in other reaches of the river would be 
beneficial as there would result savings in the cost of levee construction 
and improved drainage conditions for land both within and adjacent to the 
flood zone. 


11. For excavation of the river cutoffs the “pilot cut” method is planned 
(used most notably on the lower Mississippi where the original river length 
of 680 miles from Memphis, Tennessee to Baton Rouge, Louisiana, was 
shortened 25% through use of cutoffs). Techniques to be followed in apply- 
ing this method have been described by Matthes. (2a, 3b) 


12. The estimated average height of the levees (Table 6) is 2 m from 
La Balsa to Buga and 1.5 m from Buga to Cartago. Levees will generally 
be set back from the river-bank a distance of 25 to 50 m in order to pro- 
vide a margin of safety against bank erosion and to carry out erosion- 
preventive measures within the 25 to 50-m belt, to consist mainly of vari- 
ous types of local vegetation with deep root systems. 


Economic Justification 


As described in the introduction, the primary flood control benefit will be 
resulting from more intensive use of the land rather than from prevention 
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of damage to existing crops and improvements; this benefit should be meas- 
ured as the difference in its productive capacity under conditions expected 
with and without flood control.(4) To estimate the difference in productive 
capacity it would be necessary to determine the net income expected from the 
land and from associated resources with and without the project. However, 
because of time limitations and lack of specialized personnel, an alternative 
method was used based on estimates of the increased market value of the af- 
fected land areas. Checks on these market values were obtained by analyzing 
the productivity of the land for typical crops with and without flood control. 

As a first step, net farm income per hectare for various land classes was 
estimated in Table 7 on the basis of rental value plus added value produced 
by management. Rental values were taken at 10% of the selling prices for 
land in the Cauca Valley, for which conservative values, somewhat lower 
than currently prevailing prices, were used. (The rate of 10% is the current 
interest rate in Colombia for relatively large private enterprises.) 

Added value produced by management was taken at 30% of rental value. 
Entrepreneurs customarily pay landowners a rental equal to a third of the 
gross profit. The remaining two-thirds is for entrepreneur’s gross profit 
including bank-interest charges, management services, income taxes and net 
profit. The latter two items are conservatively estimated at 10% and 5%, 
respectively, of the entrepreneur’s gross profit; this is equivalent to 20% 
and 10% respectively of the rental value; see note 2, Table 7, Income taxes, 
while not a part of net farm income from the entrepreneur’s point of view, 
are wealth produced by the project and should therefore be included in 
evaluating benefits. 

A check on the assumed selling price for Class G land is provided by con- 
sidering current values and production costs for rice and soybeams, typical 
crops, as follows: 


1. The present average price of unmilled rice per sack of 150 pounds 
(avoirdupois) is Ps 55; for soybeans, the price per sack of 125 pounds is 
Ps 46.25. 


2. In a 3-year cycle, at least four harvests of rice and two of beans cai 
be obtained. The yields per plaza (0.64 ha) are 47.5 sacks per harvest of 
rice and 16 sacks per harvest of beans. The resulting gross annual crop 
value is Ps 4,000 per plaza. 


3. Production costs per harvest per plaza are Ps 620 for rice and Ps 
250 for beans, which gives an average production cost of Ps 1,000 per 
plaza annually. All production costs are included except entrepreneur’s 
bank-interest charges, the value of his management services and income 
taxes; these are considered part of entrepreneur’s gross profit (see item 
5). Gross profit is therefore Ps 3,000 per plaza. 


4. Allowing 10% for area occupied by roads and irrigation canals, the 
foregoing figure reduces to Ps 2,700 per plaza. 7 


5. The annual rent (see above) is a third of the gross profit or Ps 900. 
(The remaining Ps 1,800 is entrepreneur’s gross profit.) . 


6. As use of money in Colombia now costs about 10% per annum, Clas: 
G land is worth Ps 9,000 per plaza. To be conservative, a value of Ps 
7,500 per plaza was assumed. 
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TABLE 7 


NET ANNUAL FARM INCOME PER HECTARE 


(pesos) 

Added Net Net Farm 
Sale Rental Value Farm Income "B" 

Price Value Per Inceme 

Land C 

lass Per Per Plaza man Per Per 
Plaza Plaza (see Per Plaza Ha 

b note 2) Plaza 
A. Permanently flooded 500 0 0 0 0 0 

B,. Flooded 3 to 6 months 

during most years 1200 120 - 120 120 190 


C. Flooded once in 2 te 5 

years; without drainage 

and irrigation works 2400 240 70 310 310 480 
D. Flooded once in 7 to 10 

years; without drainage 

and irrigation works 3500 350 100 450 450 700 
E. Flood-free; without 


drainage and irriga- 
tion works 4000 400 120 520 520 810 


F. Flooded once in 7 to 10 
years; with drainage 


and irrigation works 6500 650 200 850 700 1100 
990(a) 
Se Flood-free; with drainege 
and irrigation works 7500 750 220 970 820 1280 
1150(a) 
Notes: 


l. Rental value taken at 10% of sale price. 


2. Added Value Produced by Management taken at 30% of rental value 
: (see text). 


}. Net Farm Income "B" equals Net Farm Income "A" less annual cost of irri- 
| gation and drainage works, taken at Ps 150 per plaza (applies only te 
| Classes F and G). The latter cost is (i) 10% of the estimated con- 
struction cost of Ps 850 per plaza of these works limited, however, to 
lateral canals and drainage ditches (main canals and pumping stations 
and the Timba project allocation to irrigation are not included) plus 
| (ii) operation and maintenance charge of Ps 65 per plaza. 

ja) Reduced 10% to allow for land occupied by irrigation and drainage 
canals. 


») One plaza equals 0.64 ha. 
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There are 85,000 ha in the flood zone of which 80,000 ha will be devoted to 
agriculture. The latter area is subdivided into land classes, with correspond- 
ing benefits per hectare, in Table 8. Also indicated in this table are the bene- 
fits in the areas adjacent to the flood zone due to lowered groundwater condi- 
tions. The ultimate direct flood-control benefit, in Part I of Table 8, is Ps 26 
million from La Balsa to Buga and Ps 12 million from Buga to Cartago or a 
total of about Ps 38 million. Following construction of Timba and the river 
improvement, the completion of the bulk of which will take 4 years, full culti- 
vation of the benefited lands will require additional time: this has been as- 
sumed to be 5 years. The ultimate direct flood-control benefit has therefore 
been discounted from Ps 38 million to Ps 32 million as shown in Part II of 
Table 8. 

Additional indirect flood-control benefits will include: value added through 
processing of increased agricultural production; reduction in loss of human 
life; improved sanitary conditions; reduction of loss of livestock by drowning 
savings in cost of highway and airport construction and maintenance resulting 
from improved flood and drainage conditions; and elimination of flood damage 
caused by interruptions to commerce and to industrial production. Although 
exact data are not available, it is believed that these benefits are worth at 
least half of the direct benefit. The total flood-control benefit then is 1.5 
times 32 or Ps 48 million. 

The annual cost of flood control was determined by adding the capital cost 
of the river improvement to that part of the Timba project capital cost al- 
located to flood control, as estimated on the basis of the “justifiable- 
alternative-expenditure method.” It was assumed that the foreign-currency 
part of the resulting total capital cost will be amortized in 20 years and that 
interest on the local-currency part will be paid indefinitely. Allowances were 
then made for depreciation and for operation and maintenance. The resulting 
annual costs were: 


Timba project Ps 4,250,000 


Cauca River improvement 6,890,000 


Total Ps 11,140,000 


As the annual benefit of Ps 48 million is about 4.3 times the annual cost of 
Ps 11.1 million, the two projects are clearly justified. 


Pilot Projects 


For planning purposes it has been assumed that complete development 7 
(including irrigation and drainage, as well as flood control) of the entire val- 
ley floor of 395,000 ha will be accomplished during a 30-year period. The 
key structure in the program is the Timba dam. Owing to its cost and multi- 
purpose nature, several years are required to arrange for its financing. 
However, in the meantime, it is possible to proceed with the reclamation of 
certain areas that are not greatly affected by floods from the Cauca River nor 
dependent on the Timba reservoir for irrigation water. Moreover, these : 
initial projects may be regarded as “pilot projects” as they will point the way 
toward efficient development of the bulk of the area in later years. 

Three such pilot projects are Roldanillo-La Unién-Toro (11,200 ha), 
Aguablanca (5,000 ha) and Bugalagrande-Zarzal (16,000 ha). Their selection : 
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TABLE 8 


DIRECT FLOOD-CONTROL BENEFIT 


(thousands of pesos annually) 


I - Ultimate Direct Benefit 


Benefit La Balsa to Buga Buga to Cartago 


Present Land 


aie per Farming F Farming 
Hectare Area Benefit Area Benefit 
ee efayic (ha) pee. = thal). ves ee 
1. Flood Zone 
A 0.81 8,000 6,500 4,000 3,200 
B 0.62 8,000 5,000 3,000 1,900 
C 0.33 23,000 7,600 8,000 2,600 
D 0.11 6,900 700 20,000 2,200 
Subtotals 45,000 19,800 35,000 9,900 


2. Areas Adjacent 
te Flocd Zone 


C 0.22 (b) 28,000 6, 200 9,000 _ 2,000 
Totals 73,000 26,000 44,000 11,900 


Total Ultimate Direct Benefit Ps 37,900,000 


(a) Difference in net farm income between Class E and present class. 
See Table 7. 


(bd) Estimated benefit due to lowered groundwater in area adjacent to 
flood zone. 


M1 - Discounted Direct Benefit 


: =e. 
Assume 5 years will be required to achieve the ultimate direct benefit. 
Use interest rate of 8%. Determine present value of the benefit for the 50 
years following project completion. 


Year % developed Present Value Product 
a 20 2926 0185 
| 2 40 2856 2342 
| 3 60 2794 476 
| 4 80 2736 - 589 
5 100 2681 2681 

| 6 to 50 100 8.247 (a) 28 247s 
10.520 


Present value of an Annuity Certain for 50 years is 12.23 
iscount is therefore: aH = 86% 
irect Flood Control Benefit is therefore 0.86 x Ps 37,900,000 or Ps32,000,000 


a) Present value of an Annuity Certain for 45 years is 12.11; 
12.11 x .681 = 8.247 


1 
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was based mainly on: extent of local interest, technical feasibility and econo- 
mic return. The first two projects are under construction while the third is 
being designed. 


Aguablanca Project 


As the land reclaimed by this project will ultimately be devoted to urban 
use (see Fig. 9), the degree of protection will be much higher than that adopt 
elsewhere in the valley. The river levee on the east is being built to a profil 
2.5 m higher than the highest recorded—during the flood of February 1950. 
The frequency of a flood 1.5 m higher (allowing 1.0 m freeboard) is estimate¢ 
to be 100 years initially; and 200 years ultimately, after completion of the 
Timba reservoir. 

Elimination of flooding within the Aguablanca project area will eliminate § 
million m3 of natural valley storage or 7.5% of the 740 million m3 between 
La Balsa and Cali. The peak discharge of the 1950 flood was 1,000 m3/s. TI 
eliminated valley storage would cause an increase to 1,032 m / s, for which 
the river level would be only 15 cm higher than under present conditions. 

Although, as described above, diking of the entire flood zone in the valley 
is not feasible, it was decided nevertheless to make an exception in the case 
of Aguablanca for the following reasons: (a) The valley storage eliminated b 
this one project has a negligible effect on river levels. (b) Benefits will be 
much higher in the case of the Aguablanca project, per hectare of reclaimed 
land, than elsewhere in the valley. (c) Further urban encroachments on the 
flood plain are not visualized, at least for the next generation or two. (The 
remaining towns are all located well outside the flood zone.) 

The interior drainage system includes (1) a pumping station to remove in 
terior drainage water during flood stages in the river and to lower the pres- 
ently excessively high groundwater level; and (2) a dredged storage pond to 
reduce pumping requirements and, at the same time, provide a permanent 
pond as a nucleus for a park development. 


Pumping station 
Design data with respect to the pumping station are as follows: 


Interior drainage area 56 km2 
Portion of drainage area urbanized, 
initially 20% 
Design Storm 
Frequency 3 years 
Sump level, minimum El 949.5 
Sump level, maximum El pie 
Peak inflow 18.0 m?/s 
Pondage 770, 000 m3 
Peak outflow mere 3/s 
Cauca River Levels 
Minimum level requiring pumping El 950.3 
Average 950.6 
3-year flood (Design ate 953.8 


100-year flood 956.0 
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Pumps 
Number of units, initially 3 
Number of units, ultimatelly 5 
Rated capacity (each unit) 27,000 gpm 
Rated head (static) 3.0 m 


The design conditions assume coincidence of a 3-year storm on the in- 
terior area with a 3-year flood in the Cauca River. The available data indi- 
cate that these two events, if not entirely independent, are almost so and that 
the probability of their coincidence is about once in 7 years. 


Variable-blade pumps 

For design conditions with river level at 953.8, static heads will be 4.3 m 
maximum and 1.9 m minimum. However, pumping will be required also for 
river levels equal to the sump level or as high as 956.0. The minimum stati 
head is thus zero and the maximum is 6.5 m. 

Because of the largely independent nature of the two hydrologic factors 
(river level and rainstorm on the interior area) the normal occurrence will 
be high rainstorm with low flood-stage and vice versa. It is therefore de- 
sirable, as a second requirement, to have a pump with a characteristic curv 
of head vs. discharge in approximately inverse relation. 

To best meet these two requirements, vertical propellor pumps with vari 
able blades are being specified. This type has the further advantage of low 
starting current, which is important at the Roldanillo-La Unién-Toro projec 
where pumps identical to that at Aguablanca will be used and where electric 
power is limited. 

Although unavailable until recently in the United States except in large 
sizes and on special order, variable-blade pumps in relatively small sizes 
(the 27,000-gpm pumps will be about 42 inches in diameter) are now being 
manufactured in the United States and Europe by several well-known firms. 
Bids recently received have confirmed that the variable-blade type provides 
the most economical solution. 


Artificial storage pond 

The pondage of 770,000 m® will be composed of 600,000 m® ina storage 
pond and 170,000 m3 in the drainage canals. If the drainage canals are late 
replaced by conduits or lined channels, thus increasing velocities and roa 
ing storage-in-transit, the storage pond can be enlarged. 

To obtain the 600,000 m3 of usable storage in the storage pond, 840,000 r 
of soil will be removed by means of a hydraulic pipeline dredge at an esti- 
mated cost of Ps 1,430,000. This expenditure is justified for the following 
reasons: 


1. Without the storage pond, the capacity of the pumping station would 
have to be doubled at a cost approximately equal to that of excavating the 
pond. If the capitalized cost of depreciation and electric energy consumr 
tion, as applied to the required additional pumping equipment, is ge . 
the economic comparison is clearly in favor of the storage pond. 


2. As the interior drainage area becomes urbanized over the course ¢ 
the years with storm sewers and streets constructed, runoff characteris 
tics will become progressively worse in the sense that the volume and es 
pecially the peak will progressively increase. The peak inflow could eas 
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increase four times, that is, up to 36 to 54 m3/s, as compared with 

18 m3/s under initial conditions. With the storage pond, little if any in- 
crease in the pumping capacity will be required to meet the needs of in- 
creased urbanization. However, to be safe, space is being provided at the 
pumping station for two additional pumping units, when and if needed, which 
would raise the rated capacity from 81,000 to 135,000 gpm. 


3. The value of the general area in which the pond will be located will 
be greatly enhanced. This area is now a swamp, 62 ha in extent. The 
completed pond will occupy only 46 ha of which the water area will be 31 
ha and 15 ha will be in islands. Moreover, the area surrounding the stor- 
age pond is being reserved for parks and housing developments. The ex- 
cavation is planned so that the shores of the pond will be well defined, 
without flat shallow places. A minimum depth of 80 cm is being provided 
to permit pleasure boating, and to prevent mosquito breeding and unde- 
sirable aquatic growths. 


Project repayment 

A law adopted by the National Government in April 1958 provides that the 
andowners benefited by the Aguablanca project shall repay the project cost to 
SVC in the form of an assessment which can be paid within optional periods 
if either 5 or 10 years. For the shorter term, interest of 8% will be charged 
ind, for the latter, 10%. (The usual interest rate in Colombia for commercial 
-redits is 10%.) In addition, landowners must pay CVC 10% of the enhanced 
ralue of the land benefited by the project, as determined by a committee of 
hree assessors chosen jointly by CVC and the landowners. 

Although many landowners are opposed to this formula, they should have no 
lifficulty in paying the assessment it calls for. Unimproved flood-free land 
.ear the project area has a market value of 15,000 to as high as Ps 50,000 
yer ha with an average of about 20,000. In comparison, the average value of 
he land in the flood zone is about Ps 7,500 per ha. 

To assist those landowners claiming hardship, the law of April 1958 con- 
ains a provision whereby landowners may elect to pay a part of their as- 
‘essment in the form of land. 

The project benefits, as measured by the enhancement value averages 

°s 12,500 per ha or Ps 72,350,000 for the project area of 5,800 ha (5,000 ha 
vithin the project proper plus 800 ha outside that will benefit because of im- 
yroved drainage). The estimated project cost including rights-of-way is 

2s 15,060,000 giving a benefit-cost ratio of 4.8 to 1. 


toldanillo-La Union-Toro Project 

An interesting aspect of this project, Fig. 10, is the dual-purpose nature 
if two of its features: the La Union Pumping Station and the intercepting 
anal 


e La Union Pumping Station will contain three 27,000-gpm units to pump 

ainage water to the Cauca River during flood stages; also, to pump water 

om the principal drain to irrigation canals located on the high bank of the 

ver. As static lifts are similar to that at the Aguablanca project, it has 

n possible to adopt identical requirements for the pumping units at both 
ons. | 

The intercepting canal will be of relatively large capacity in order to 

in dle flood flows from streams draining the east slope of the Cordillera 
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yecidental; this canal will also serve as the main irrigation-supply canal for 
he western half of the area. To accomplish this second purpose, the canal 
7ill be provided with check structures each including a tainter gate operated 
utomatically to maintain water levels sufficient for irrigation; during floods, 
he gates will open completely so as to provide unobstructed passages. 

The natural levee along the left bank of the Cauca River provides protec- 
ion against floods up to an estimated frequency of once in 5 years. When the 
rrigation canals on the Cauca River bank are constructed, the excavated 
aaterial will be placed to form a levee 1 to 2 m high; the protection then 
vailable will be against about a 7-year flood. This degree of protection is 
onsidered adequate in view of the low degree of present land use and intense 
esires of the small landowners in the area to cultivate their land. The future 
‘auca River improvement will increase the degree of protection to about 10 
ears. 
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Abbreviations and Equivalents 


centimeters = 0.394 inches 
degrees Centigrade 
elevation in meters 

U.S. gallons per minute 
hectare = 2.47 acres 
kilometers = 0.621 miles 
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km2 square kilometers = 0.386 square miles 
kw kilowatts 
m meters = 3.28 feet 
m> cubic meters ; 
mm millimeters (100 mm = 3.94 inches) 
Mm? million cubic meters = 811 acre-feet 
m3/s cubic meters per second = 35.3 cubic feet per second 
plaza 0.64 ha 
Ps Colombian pesos 
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FRICTION FACTORS IN CORRUGATED METAL PIPE 


Marvin J. Webster,! F. ASCE, and 
Laurence R. Metcalf,2 M. ASCE 


SYNOPSIS 


A study made to determine friction coefficients for flow in large corrugated 
tal pipe by measurement of head losses in pipes 3, 5, and 7 ft in diameter 

d to determine velocity distribution for several discharges in each test pipe 
that the results could be generalized to show the relationship between fric- 
m coefficients and absolute roughness is discussed in this paper. Data taken 
ring the study are presented together with recommended values of friction 

d roughness coefficients. 


INTRODUCTION 


‘Numerous laboratory and field tests have been made to determine friction 
stors for various types and sizes of pipes and conduits. However, compara- 
‘ely few tests had been made prior to 1950 on corrugated metal pipe of the 
rger sizes commonly used in engineering practice. For this reason, the 
fice, Chief of Engineers, U. S. Army, authorized tests on 3-, 5-, and 7-ft- 
um corrugated metal pipe to be conducted at the Bonneville Hydraulic 
boratory. 
In 1926, Yarnell, Nagler, and Woodward published results of hydraulic 
sts(1) at the State University of Iowa on 12-, 18-, 24-, and 30-in.-diam cor- 
ted metal pipe that varied in length from 24 to 36 ft. The results of these 
are applicable to similar installations of short pipes and culverts, such 
those under highway and railroad fills. However, they are not directly ap- 
cable to long pipelines because entrance, transition, and exit losses cannot 


é: Discussion open until February 1, 1960, To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 2148 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 9, September, 1959. 

Chf. Hydr. Design and Lab. Section, U. S. Army Engr. Dist., Portland, 
Corps of Engrs., Portland, Ore. 

Hydr. Engr., U. S. Army Engr. Dist., Portland, Corps of Engrs., Portland, 
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be isolated from the pipe friction losses when the test pipes are of such sho! 
length in relation to their diameter. Keulegan(2) used velocity distribution 
data from the foregoing tests to compute friction factors for long pipelines 
and for pipes of larger diameter. 

In 1950, Straub and Morris published results on studies(3) conducted at th 
St. Anthony Falls Hydraulic Laboratory, University of Minnesota, of flow 
through 18-, 24-, and 36-in.-diam circular and arched corrugated metal pips 
193 ft in length and laid on a slope of 0.002. These studies showed a definite 
systematic increase in the values of the friction and roughness coefficients 
“f” and “n” with increased Reynolds numbers (Re) over the entire range of 
discharges tested. These findings were contrary to expectations because 
most “f” vs Re curves for both smooth and rough commercial pipe have 
shown constant or slightly decreasing characteristics throughout the experi- 
mental range. Values of “n” from 0.0216 to 0.0252 were obtained for cor- 
rugated metal pipe flowing full and from 0.0228 to 0.0258 for open-channel 
flow. Values recommended for use in design were 0.025 and 0.0235 for full- 
pipe and open-channel flow, respectively. 

The purpose of the tests at the Bonneville Hydraulic Laboratory was to 
determine friction factors for 3-, 5-, and 7-ft-diam corrugated metal pipe a 
indicated by head-loss measurements for a range of velocities up to approxi 
mately 10 fps for 5- and 7-ft-diam pipe and 16 fps for 3-ft-diam pipe. Test: 
were made with both full-pipe and open-channel flow in the 3- and 5-ft-diam 
pipe and with full-pipe flow in the 7-ft-diam pipe. The initial test schedule 
was expanded to include studies on the 5- and 7-ft-diam pipe with 25 per cer 
of the periphery (at the invert) smooth paved and on the 5-ft-diam pipe with 
50 per cent of the periphery smooth paved. 


Description of Project 


The tests were conducted at a site located on the Washington shore of the 
Columbia River at the north end of Bonneville Dam. This site was selected 
because sufficient space was available for the test and metering sections an 
because sufficient head and water supply were available to produce the dis- 
charges and velocities required for the tests. The layout for the 7-ft-diam 
pipe, typical for the three sizes of pipe tested, is shown in Fig. 1. The 3-ft- 
diam pipe in the test section is shown in Fig. 2. Water for the study was 
taken from one of the fingerling by-pass channels through the main dam. 
From the fingerling by-pass, the water flowed through a 6-ft-diam, smooth- 
wall pipe to the test section. From the test section, the flow passed through 
a 7-ft-diam, smooth-wall pipe and then dropped 11 ft to the metering sectiol 
where it was measured by a system of flat-plate orifices and then returned. 
to the river through gate valves and short outlets. The 11-ft difference in 
elevation between the test and metering sections was sufficient to provide — 
adequate head on the meters for the tests with open-channel flow. 

The corrugated metal test pipes were of standard commercial manufac- 
ture and arrived at the test site in 20-ft-long pieces composed of a series 0 
2-ft-long overlapping sections. The sections were riveted together and eac! 
contained lapped and riveted longitudinal joints that were staggered so that | 
two were in line in adjacent sections. Before installation, the riveted joints 
were sealed against leakage by welding on the outside of the pipe. The 20-f 
sections were then butt-welded together to form a pipeline that was set to — 
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grade and alignment along a continuous footing between concrete saddles. Ok 
servations were made when the pipes were flowing full but no sagging or de- 
formation was visible. The average diameters (top of corrugation to top of 
corrugation) of the 3-, 5-, and 7-ft pipe were 2.996 ft, 4.95 ft and 7.05 ft, 
respectively. The depth and longitudinal spacing of the corrugations on all 
pipes were 1/2 in. and 2-2/3 in., respectively. 

The invert paving of the 5- and 7-ft-diam test pipe was accomplished by 
welding 18-gage, 10-ft-long sheets of metal over the corrugations on the bot- 
tom of the pipe (see Fig. 3). The lower sector of the pipe was then coated wi 
bituminous paint to complete the simulation of the conventional asphaltic pav- 
ing in use today. A compromise in the reproduction of the paving was neces- 
sary along the side joints where the sheet metal met the corrugations; these 
joints were left open across the valleys of the corrugations. 


Instrumentation 


The rate of flow was measured in a metering section that was composed o 
a parallel system of flat-plate orifice meters designed to measure a continu- 
ous range of discharges from 8 to 385 cfs and to provide an overlap between 
meters operated with a differential head of from 1 to 9 ft. An aperture-plate 
diameter ratio of 0.7 to 1.0 was selected for all orifices except the one in the 
5-ft section; this orifice had a 0.66-to-1.0 ratio to secure slightly different 
discharge characteristics. The orifices were machined from steel plates th 
varied in thickness from 1/2 in. for orifices up to 3 ft in diameter to 1 in. fo: 
those larger than 5 ft in diameter. Piezometer taps were located one pipe 
diameter upstream from the orifice plate and one half the orifice diameter 
downstream from the plane of the orifice. The meters were assembled prio1 
to rating and were moved intact to the metering section. Periodic inspectior 
were made of the orifice plates and piezometer taps to insure that they re- 
mained clean and free from rust. In the metering section, the orifice meter 
in the 6-ft-diam pipe (Fig. 1) was preceded by 19 and followed by 6 diameter: 
of straight, constant-diameter, smooth-wall pipe; in all of the other meter 
installations, the approach and exit pipes exceeded the length standards pre- 
scribed by the Special Research Committee on Fluid Meters of the American 
Society of Mechanical Engineers.(4) First, two 24.0- by 16.2-in. orifice 
meters were calibrated separately, using a standard calibrated venturi mete 
and then installed in the metering section. Successively larger orifice mete) 
were calibrated as follows: the entire length of constant-diameter, smooth- 
wall pipe to be used upstream from each orifice and the 20-ft-long section 
that contained the orifice were placed in the test section, the orifice piezo- 
meter taps were attached to a standard differential water manometer, and _ 
varied rates of flow were introduced into the system. The quantity of flow, 
as indicated by previously calibrated orifices in the metering section, to- 
gether with the observed head differential on the test orifice, provided the 
basic data from which experimental coefficient curves were plotted. After © 
all data were checked for accuracy and consistency, each orifice meter and 
adjacent pipe were moved to the metering section and installed for further b 
use. “ 

The test sections in the 3-, 5-, and 7-ft-diam pipe were 180, 251, and 350 
ft long, respectively. Pressure observation points for the 7-ft-diam pipe ar 
indicated in Fig. 1. There were 12 stations in the 3- ft pipe at approximatel: 


FRICTION FACTORS 


3 Five-Ft-Diam Pipe With 50% Paving 


39 


40 September, 1959 HY 9 


15-ft spacings and 10 stations in the 5-ft pipe at approximately 20-ft spacings. 
Each pressure observation point consisted of two 1/8-in. -diam piezometers 
that were drilled into the valleys of the corrugations 45 degrees above the pips 
invert on each side. The two piezometers were interconnected by a common 
lead to a single manometer in the bank shown in Fig. 4. Individual pressure 
fluctuations at each piezometer varied from 0.005 ft at low discharges to 
0.015 or 0.02 ft at high discharges, however, it was found that consistent re- 
sults could be obtained from a series of 5, 10, or 20 consecutive visual ob- 
servations. 

It was determined in a preliminary study on 3-ft pipe (Fig. 5) that these 
piezometers were influenced by the partial recovery of velocity head in the 
corrugations and therefore indicated pressures higher than the true static 
pressure as taken at a point 6 in. within the pipe. However, they were con- 
sidered sufficiently accurate to determine the differential pressure between 
successive stations along the pipe axis. The data in Fig. 5 are qualitative in 
nature and are included to show that the uniformity of pressures measured 
along the corrugations depends on exact placement of piezometers. It is be- 
lieved that the nonuniform points shown on some of the grade lines presented 
in Fig. 7 were caused either by abnormal velocity distribution or by imper- 
fections within the pipe that affected localized velocities and pressure adja- 
cent to the piezometers. 

Static pressures during the velocity distribution studies and free water- 
surface levels for open-channel flow were measured by means of tubes that 
were the equivalent of the static leg of a conventional pitot tube (Fig. 6). The 
static tubes were checked to determine the effect of local flow acceleration. 
At first these tubes contained four piezometers that were connected to indi- 
vidual water manometers, and the instrument was aligned with the flow by 
rotating the tube until the pressure, as reflected by the manometers, was 
equal at all piezometers. Prior to the measurement of hydraulic grade lines 
in the 5-ft-diam paved-invert pipe and the 7-ft-diam pipe, the static tubes 
were revised so that the four piezometers in each tube were connected to a 
common manometer that reacted much faster to pressure changes than did 
the instrument of original design. 

Two types of impact tubes were used to measure local velocity heads from 
which velocity distribution curves were determined. The type used first re- 
sembled a pitot tube but lacked the static openings on the sides of the barrel. 
The tube was mounted in a bracket that slid along a strut placed across the 
pipe diameter. The impact tube was satisfactory during tests on the 3-ft pipe 
although at high test velocities it vibrated and was difficult to keep pointed 
into the flow. It was found that the tube would require additional support 
brackets when velocity heads were measured in the lower half (extended posi- 
tions) of the 5-ft pipe. Because it would be necessary to unwater the pipeline 
during each run in order to install the additional brackets, the impact cylinde1 
(shown in Fig. 6) was substituted for the impact tube. This cylinder was sup- 
ported by brackets mounted 12 in. apart on the support strut. It was found to 
be of excellent accuracy as a measuring device and did not turn from the flow. 


Procedure 


Friction losses for pipes flowing full were determined in the following 
manner; after the flow was adjusted to a desired rate through one or more of 
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NOTES 


- READINGS TAKEN AT STA 1 +50. 
2. VELOCITY IN THE PIPE 10 FPS. 
+ ZERO EQUALS STATIC PRESSURE 


(6 IN. INSIDE PIPE). PRESSURE DISTRIBUTION 
4. PIEZOMETERS ~ IN. DIAMETER, ALONG CORRUGATIONS 
PLACED 45DEG ABOVE PIPE INVERT. 
5. DEPTH OF CORRUGATION SIN; CYCLE 3-FT-DIAM PIPE 


LENGTH 2£IN. 


FIG. 5 
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Fig. 6 Velocity and Pressure Instruments 
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ihe orifice meters and respective gate valves, and the average differential 
1eads on the meters were derived from a series of 10 readings, the two ob- 
servers proceeded to the test section and measured the hydraulic grade line 
within the test pipe by means of 10 consecutive readings of the water-surface 
elevations within the manometer tubes. Then the observers returned to the 
metering section and checked the meter readings. In all tests, if the average 
lifferential head on the meters at the beginning and end of a test run agreed 
within 1 per cent, the average of the two heads was used to compute the rate 
of flow; otherwise the entire run was discarded. 

Velocity distribution in the 3-ft-diam pipe was measured by means of a 
-ombination of one of the previously described static tubes and two impact 
ubes at each observation station. Leads from the three tubes were connected 
0 a 3-leg differential manometer where velocity heads were recorded in 
sroups of 10 simultaneous readings. The static tube was mounted 6 in. within 
he pipe and 45 degrees above the invert, whereas the two impact tubes were 
raversed through the horizontal and vertical diameters on separate struts. 

n order that the proximity of the horizontal strut might not affect the ac- 
‘uracy of observations, the horizontal strut was removed while the vertical 
raverses were being made. The velocity distribution in the standard 5-ft- 
iam pipe was measured by means of three impact tubes and one static tube 
hat were used simultaneously at each station. One impact tube was used on 
ach half of the horizontal traverse and the other on the vertical traverse. 
nstrumental data were observed 20 times in succession to determine the 
verage velocity for each position of the tubes. Two impact cylinders and a 
tatic tube were used to measure velocities in the 5-ft-diam pipe with a paved 
avert and in the 7-ft-diam pipe. 

_ For tests of open-channel flow in the 3- and 5-ft-diam pipes, the approxi- 
nate rate of flow was regulated at the intake gate; and the depth in the test 
ipe was balanced by adjusting the valves in the metering section until the 
ydraulic grade line, as recorded by the same piezometers or static tubes 
Sed for the studies of full-pipe flow, was nearly parallel to the pipe invert 
lope. Then the depth of flow was measured by means of a static tube located 
ear the longitudinal center of the test section and the rate of flow was deter- 
rined in the metering section. It was difficult to maintain a constant slope of 
ue water surface in the test pipes, and continuous fluctuations of the forebay 
‘vel made it impossible to conduct tests of open-channel flow during certain 
eriods when Bonneville Dam was used to absorb peaks in the regional power 
emand. Tranquil flow was required to avoid the occurrence of standing waves 
. the test section; therefore, it was necessary to decrease the slope of the 
-ft-diam, paved-invert pipe from 0.005 to 0.002. The same change and pos- 
bly revisions of the intake structure would have been required to permit 
ists of open-channel flow in the 7-ft-diam pipe. 


Conputations 


ll-Pipe Flow 


The average velocity of flow for both standard and paved pipes flowing full 
based on the measured discharge and on the area corresponding to pipe 
eters of 3.00, 4.95, and 7.05 ft. In the case of the paved pipe, the area 
ductions effected by the paving are reflected in the values of “f” and “n.” 

e values of “S” were then determined from straight-line portions of the 
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hydraulic grade line, typical examples of which are shown in Figs. 7 to 9.~ it 
will be noted that the grade lines for all tests are fairly straight with the ex- 
ception of the 3-ft pipe tests shown in Fig. 7. In this case, only the straight 
portions of the grade lines between measuring points 3 and 10 were used to 
compute the slope. The other points, near the ends of the test section, were 
apparently affected by entrance and exit conditions. Values of “S” obtained 
from these grade lines were plotted against measured discharges to form the 
rating curves shown in Fig. 10 and were used to derive numerical equivalents 
for “f,” “n,” and Re that are shown in Figs. 11 and 12. The friction factor 
“f” was computed from the Darcy formula, in which 


the roughness coefficient “n” for the Manning formula, 
V = 1.486 R2/3 31/2 
n 


was computed from the friction factor “f” by the relationship 


10.8 


and Re was determined from the formula 


R, = VD 
Vv 
Typical curves of velocity distribution within the test pipes are shown in 
Figs. 13 and 14. Local velocities were computed from the formula V, = v2gh 
where “h” was the mean difference between readings obtained from the pre- 
viously described impact and static tubes. In the vertical traverse in the 3- 


ft-diam pipe, where the impact tube was displaced downstream from the stati 
tube, a correction of 


2 
He =f x s (where L equals the horizontal 
8 displacement) 


was applied to the value of “h” before the local velocity was computed. Also 
shown in Figs. 13 and 14 are curves of the relationship V1/Vx; as computed 


from the formula: V, =VRSg = VV£/8. This relationship was one of the 
means used to eaalvae the adequacy and consistency of the velocity data. 


Open-Channel Flow 


The first step toward determination of “n” values for open-channel flow in 
the 3- and 5-ft-diam pipes was the computation of cross-sectional areas and 
hydraulic radii for observed depths and rates of flow in the test pipes. Then 
slopes were determined from straight-line portions of plotted hydraulic ered 
lines, and values of “n” were computed from the Manning formula. The 
curves presented in Fig. 15 show the variation of the roughness coefficient — 
with the depth of flow. In order to plot rating curves for open-channel flow, 
the observed measurements were reduced to equivalent values for a slope of 
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0.005. The discharge quantities, as derived from the Manning formula, are 

proportional to the square root of the slope; therefore, measured discharges 
0.005 1/2 eek 

measured slope ) to obtain discharge 

quantities for the rating curves. Average flow depths, when plotted against 

these discharges, provided the curves shown in Fig. 16. 


were adjusted by the factor ( 


Test Results 


The tests showed that the value of Darcy’s “f” for standard corrugated 
metal pipe varied inversely with the pipe size (from 0.077 for 3-ft pipe to 
0.054 for 7-ft pipe; Fig. 11) and that the value of Manning’s “n” was nearly 
constant (0.024; Fig. 12). The values of both coefficients for each pipe size 
varied with Reynolds number and reached a maximum within the test range. 
The invert paving reduced the resistance to flow and caused lower values for 
both coefficients (“f” values of 0.048 and 0.042 for the 5- and 7-ft pipes 25 
per cent paved and 0.035 for the 5-ft pipe 50 per cent paved; “n” values of 
0.021 for the 25 per cent paved 5- and 7-ft pipes and 0.018 for the 5-ft pipe 

50 per cent paved). Values of the plotted points shown on the foregoing 
figures are presented in Tables Ato G. Typical curves of velocity distribu- 
tion within the respective test pipes are shown in Figs. 13 and 14. 

The value of Manning’s “n” for open-channel flow in standard pipe was 
found to be almost independent of the depth of water over the range tested 
and was approximately equal to the value (0.024) obtained for full-pipe flow 

(Fig. 15). The values of “n” increased as the depth of flow increased in the 
/paved-invert pipes. The curves shown in Fig. 15 were extrapolated down to 
an assumed value of *n” = 0.011 to compute the lower parts of the curves in 
Fig. 16 when the flow was contained within the paving. This value is com- 
monly used when computing the flow of water through welded steel or cast 
‘iron pipe in good condition. 


Discussion 


The over-all results of the study confirmed the findings of the tests made 
at St. Anthony Falls Hydraulic Laboratory(3) that friction factors for a given 
size of corrugated metal pipe are a function of Reynolds number (Figs. 11. 
‘and 12) and were in general agreement on values of “f” and “n” within the 
/scope of that study. However, the tests at St. Anthony Falls were concluded 
‘within the range where “f” and “n” increased with Reynolds number and im- 
plied that the coefficients would continue to increase indefinitely, whereas 
this study determined maximum values of the coefficients for each pipe size. 
From a practical standpoint, and except in the most critical situations, these 
maximum values are adequate for design purposes. In a theoretical applica- 
tion, however, such as extrapolating these values to-other sizes of pipe, it 
should be remembered that flow in corrugated metal pipe does not neces- 
Sarily follow the laws established for flow in other types of pipe, Morris(5) 
has established that “wake-interference” flow is obtained on corrugated sur- 
faces of commercial conduits, whereas “isolated-roughness” flow occurs over 
the surfaces of concrete, riveted steel, wood-stave, cast-iron, and brick con- 
duits. 
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TABLE A 


SUMMARY OF TEST RESULTS 
FULL-PIPE FLOW 


3-FI-DIAM CORRUGATED PIPE 
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69 


. e 
Fee eh ORO 20 
NO3II3 


REEL 
pat pee 


(CES) (2) Temp no73 f 
(Deg. Fah) 

(2) (3) ( (5) (6) 
6.72 9.0331 66.1 258 0.0706 
6.81 20351 66.1 261 20730 
6.89 .0360 66.1 26) 20732 
6.95 -0360 70.3 281 -0720 
6.95 -0360 66.1 267 -0720 
7.05 -0370 68.2 270 -0723 
33 -9370 60.1 257 .0663 
7.75 -O419 60.4 273 .0671 
8.20 ~087 59.5 285 0698 
8.54 -0516 60.4 301 .0682 
8.75 -0560 59.8 305 .0706 
8.88 2958) 66.44 3h0 -0715 
9.26 -0672 55.0 301 -0758 
9.63 9701 60.) 339 0729 
9.75 -0711 59.8 31 -0720 

10.45 -0818 58.3 357 0723 
10.45 -0799 66.1 01 20705 
10.45 -08)7 66.1 01 -07k7 
10.45 -0808 66.1 ol -0713 
10.45 0828 66.5 399 20731 
10.45 0808 66.7 o2 20713 
11.96 21091 59.0 413 .0735 
12.26 .1208 55.0 399 -077h 
12.20 1130 66.7 469 0732 
12.35 .1169 66.2 75 -0739 
12.35 -1159 66.4 475 -0733 
12.51 -1169 66.4 479 -0720 
13.14 1334 59.1 yoy 2076 
13.78 1461 67.0 532 -07h2 
13.93 -1568 68.5 Sug -0730 
13.9) .1568 66.1 535 .0778 
13.9h 1519 66.1 533 -0754 
13.96 «1509 70.3 56h, .0748 
15.52 .18h0 55.0 505 -0736 
16.02 .2016 60.1 562 .0758 
w7eee 237 66.3 662 -0763 
L7h2 +2435 66.2 666 077 
17 lh - 2hbh 66.7 671 -0775 
174k 205 66.5 670 -0762 
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27.88 
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TABLE B 


SUMMARY OF TEST RESULTS 
FULL=-PIPE FLOW 


S-FT-DIAM CORRUGATED PIPE 


vV So Water Re 3 
(FPS) (CFS) (2) Temp Dake) 
(Deg. Fah) 
3.07 59.1 0.18 hho 978 
4.06 78.1 “By hh.o 1291 
56 87.8 2 Te) 1452 
5.08 97.8 52 hh.o 1618 
5.16 105.1 é Who 1753 
5.61 108.0 -63 4h.o 1786 
6.12 117.9 75 hh.0 19h9 
6,81 131.0 093 Mh.5 2181 
7.39 12.2 1.08 Uh.s 2368 
75k LiSe2 1.15 hh.5 216 
7.40 1h2.h 1.11 hh.o 2356 
8.28 159.) 37 Ae) 2636 
8.76 168.6 1.53 LO 2788 
9.2h 7 7eo 1.70 TRRe) 292 
: 189.8 1.93 hh.o 3139 
10.28 19Te7 2.10 45.0 332 
10.28 197.7 2.07* 45.0 332 
11.01 211.9 2.37 45.0 3562 
11.01 211.9 2.35% 45.0 3562 
3.02 58.0 o17* 43.5 953 
(3.96 76.2 31 43.5 1252 
Deane) 8.7 038% 43.5 1390 
5.52 106.2 61% 36 173 
291 11368 - 70% [335 1868 
9.30 178.9 1.70# 43.5 2936 
> 9.53 183.3 1.80% 43.5 3009 
10.39 199.9 2.13% 43.5 3282 
12.15 21.5 2 ox 43.5 3521 
*% Head loss measured with static tubes. 
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TABLE C 


SUMMARY OF TEST RESULTS 
FULL-PIPE FLOW 


5-FT-DLAM CORRUGATED PIPE 
WITH 25 PER CENT OF PERIPHERY (INVERT) PAVED 


V NS) Water Re_, 

(FPS) (CFS) (4) Temp x10 f n 
(Deg. Fah) 

.08 3 0.143 71.0 146) 0.0482 0.0211 
3°52 ee 166 70.5 1570 -0479 .0210 
3.4h7 66.8 180 71.0 1651 -0478 .0210 
3.59 69.1 19h 711.0 1709 -04,80 .0210 
3.73 71.8 209 71.0 1770 +0479 .0210 
3.89 7.9 232 711.0 182 -0489 0212 
3.91 75.2 231 71.0 183 .0477 .0210 
4.00 77.0 222 711.0 190) -O481 .0211 
Tea Ke) 8.7 6292 70.5 207k -0480 .0210 
h.71 90.6 33h 70.5 2200 20479 .0210 
5.08 97.8 387 71.0 2h13 0473 .0209 
5.h7 105.3 hh6 70.5 2579 -O474 .0209 
5.80 111.6 498 70.0 2726 -O471 -0208 
6.24 120.1 572 71.0 2967 -0468 .0208 
6.38 122.8 -603 70.0 298), -0472 .0208 
6.91 133.0 -698 70.0 3233 066 .0207 
7.00 13h.7 727 68.0 3193 -04,72 20209 
7.06 135.9 0137 72.0 3393 0471 0208 
7 ll 12.6 80) 72.0 3561 - Oh .0207 
7.48 143.9 2825 68.0 3h12 20469 .0208 
7.98 153.6 90 68.0 360 -0L70 .0208 
8.16 162.9 1.048 68.0 3859 066 .0207 
8.97 172.6 1.168 68.0 090 -063 .0207 


A summary of recommended values of friction and roughness coefficients 
for use with corrugated metal pipes flowing full is presented in Fig. 17, in 
which the “n” curve for unpaved pipe was extrapolated to include pipes of 
from 1 to 10 ft in diameter, and values for the “f” curve were computed from 
values of “n.” Points from the tests at St. Anthony Falls (1.5-, 2-, and 3-ft 
pipes) are included to show the validity of the curves for pipes of smaller _ 
diameter than those tested during the present study. The values for “n” shov 
a slight decrease for pipes of large diameter — the over-all change between 
1- and 10-ft diameters is from 0.025 to 0.023, whereas “f” values decline 
sharply from 0.117 for 1-ft-diam pipe to 0.046 for 10-ft-diam pipe. Parallel 
curves for the 25 per cent paved pipe are shown, in which “n” equals approxi 
mately 0.021 and “f” varies from 0.053 to 0.038 for pipes 4 and 8 ft in diame 
ter, respectively. No curves are shown for 50 per cent paved pipe because — 
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TABLE D 


SUMMARY OF TEST RESULTS 
FULL-PIPE FLOW 


5-FT-DIAM CORRUGATED PIPE 
WITH 50 PER CENT OF PERIPHERY (INVERT) PAVED 


v Q s Water Re 5 
(FPS ) (CFS ) (%) Temp xo f n 
(Deg. Fah) 

3.18 61.2 0.112 56.0 122) 0.0354 0.0181 
3.2 65.9 129 56.0 1318 -0351 -0180 
3.51 67.6 137 55.0 1332 0356 0181 
3.6 70.0 lh) 56.0 102 -0346 -0179 
3.7h 72.0 153 55.0 1418 -0347 -0179 
3.88 7.6 165 56.0 1495 -0349 20179 
4.06 78.0 -180 55.0 1540 -0348 -0179 
lo 86.) .222 55.0 1700 20351 -0180 
5.00 96.1 2275 55.0 1900 -0350 -0180 
5.55 106.8 30 55.0 2110 0352 -0180 
5.56 107.1 3hh 55.0 2110 -0355 -0181 
5.98 115.1 2393 53.5 2220 -0350 -0180 
6.02 115.8 395 55.0 2280 -0348 .0179 
6.49 124.9 453 53.5 2h10 032 -0178 
6.51 125.3 -457 55.0 2470 -03L) -0178 
6.99 134.5 -520 53.5 2590 0338 -0177 
7.52 144.7 602 5365 2790 0339 0177 
7.98 153.6 2675 53.5 2960 -0338 20177 
8.00 154.0 ‘ 5h.0 2960 -030 .0177 
8.48 163.1 «765 53.5 31h0 -0339 20177 
9.07 17.6 869 53.5 3360 -0336 0176 
9.58 184.3 2965 53.5 3550 0335 -0176 
10.50 202.1 1.12 53.5 3880 -0330 20175 


only one size was tested; for 5-ft pipe the values of “n” = 0.018 and 
*f” = 0.035 are recommended. 
It was recognized that several factors might have affected the accuracy of 
the test results. Perhaps the most important of these factors involved the 
practice of calibrating successively larger orifice meters against several 
smaller ones. However, every known precaution was taken to ensure that all 
orifices were rated and used correctly and that they were kept in perfect 
condition throughout the course of the study. In addition, the random nature 
of the test pipes prevented an exact study of flow characteristics. For ex- q 
ample, a 6-ft-long straight-edge laid along the inside of the respective pipes 
often revealed that one or two adjacent corrugations extended as much as 3/: 
in. farther into the pipes than did the rest. Such inequalities within the test 
pipes probably had no significant effect on friction or roughness coefficients, 
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TABLE E 


SUMMARY OF TEST RESULTS 


FULL-PIPE FLOW 


7-FI-DIAM CORRUGATED PIPE 
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TABLE F 


SUMMARY OF TEST RESULTS 
FULL-PIPE FLOW 


7-FT-DIAM CORRUGATED PIPE 
WITH 25 PER CENT OF PERIPHERY (INVERT) PAVED 


V Q 5 Water eh 
(FPS ) (CFS) (4) Temp x10 f n 
(Deg. Fah) 

3.06 120 0.085 57.5 1727 0.0410 0.0206 
3.35 131 102 57.5 1885 20411 .0206 
3.68 1h) 126 ores 2072 O21 .0209 
3.94 154 11 57.5 2216 -012 -0207 
rea? 163 158 58,5 237k 0413 .0207 
hebh 173 178 58.5 2519 010 .0206 
be 175 181 BY AS: 2518 0,09 .0206 
), .87 190 .222 58. 2767 2Oh2h .0209 
5.19 202 2h7 58.5 292 -O16 -0208 
5.60 219 291 58.5 3189 oL20 .0208 
5.79 226 318 58.5 3291 0432 -0211 
6.69 261 21. 58.5 3801 -0L26 .0210 
feats) 277 67 58.5 03h 04.20 .0208 
ene 293 S21 58.5 267 018 -0208 
7.15 303 50 58.5 413 04,08 0205 
7.86 307 Ee) 57.5 4h17 -0396 .0202 
8.11 B17 -596 58.5 4,616 -OW11 -0206 
8.26 323 602 57.5 L648 0400 -0203 
8.58 335 64 bye 4820 -0397 -0202 
8.93 39 -693 57.5 5016 039h -0201 
9.30 363 2750 57-5 5223 0393 0201 
9.62 376 802 57.5 5410 0393 .0201 


but their presence would explain the nonuniform distribution of flow that 
sometimes was evident in the velocity traverses. 


CONCLUSIONS 


The results of the study indicate that a closed correlation existed between — 
Reynolds number and both Darcy’s friction coefficient “f” and Manning’s ; 
roughness coefficient “n” throughout the range of experimental discharges ¥ 
reproduced for study. The following values of Manning’s coefficient “n” were — 
obtained from tests on standard and paved-invert corrugated metal pipes ; 
flowing full: 
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TABLE G 


SUMMARY OF TEST RESULTS 
OPEN-CHANNEL FLOW 


Q D R A V iS) Water 
(CFS) (FT) (FT) -(SQ FT) (FPS) (%) Temp 
(Deg. Fah) 
3-FT=DIAM CORRUGATED PIPE 
19.50 1.90 0.852 71 1h 0.526 63.7 
23,01 © 2.25 0.903 5.67 4.06 0.510 61.2 
25.93 2.59 0.907 6.49 4.00 0.491 66.3 
25.00% 3.00% 0.750% 7.07%  3.5)% 0.500% 62.0% 
S-FT-DIAM CORRUGATED PIPE 
32.7 1.98 1.06 ae yeh 0.513 8.0 
39.9 2.2) 16 8.5 4.69 0.500 48.0 
41.4 2.25 Thea ll 8.5 ), 87 0.519 48.0 
52.5 2.61 1.28 10.3 5.10 0.500 48.0 
60.2 2.85 1.34 11.4 5.28 0.500 49.0 
2 3.18 1.42 SL Peanin 0.50 50.0 
83.2 3.56 1.48 14.8 5.62 0.500 50.0 
91.0 3,82 1.50 16.0 5.69 0.513 50.0 
90.9 3.8h 1.50 16.0 5.68 0.500 53.0 
100.) 42 1.48 18.1 5.55 0.500 53.0 
96.2% h.95t 1.2% 19.2% 5.Ol# 0.500% hh .O% 


S-FT=-DLAM CORRUGATED PIPE 


WITH 25 PER CENT OF PERIPHERY (INVERT) PAVED 


30.6) 2.00 1.07 7.29 4.20 0.200 63.0 
MhO.S. 2.40 1.21 9.25 43h 0.200 63.0 
\eeng.ch «862.70 1.31 10.73. &.h0 0.200 63.0 
| 53.80 3.00 abyss: 12.20 41 0.200 64.0 
i 58.46 3.31 1.45 13.67 4.48 0.200 64.0 
B2ehy <63.650 1.647 -1.53 4.33 0.200 64.5 
66535 3.12 1.50 15.53 4.27 0.200 64.5 
70.65 h.06 1.51 16.88 4.19 0.200 65.0 
Mm 72.61 4.06 1.51 16.88 4.30 0.200 65.5 
We73euO ~~ .22 1.50 17.49 4.20 0.200 64.0 
| 70.2% .95* 1.2h% 19.2h% 3.65% 0.200 71.0% 
| S-FT-DIAM CORRUGATED PIPE 
| WITH 50 PER CENT OF PERIPHERY (INVERT) PAVED 
76.8 3.01 1.38 12.25 6.27 0.192 48.0 
Meera.) «3017. «123.00. G07 0.180 49.0 
eee 3,09) 851 n0. os 12.69°" 6.19 0.190 4.8.0 
ery 3571-1507 65.50". 5.62 0.185 49.0 
87.9 387i or 50ge 15.50 5.67 0.200 49.0 
88.1 SBlyues 42-65], wth. 00x 145.50 0.190 49.0 
Bie e237 ls 17550> “5.1 0.194 50.0 
81.8% «= .95*# 1 619.2% = 0.200% 55.0% 


# Estimated from full-pipe data. 
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Nominal Pipe Diameter Coefficient 
Type of Pipe in Feet hi 
Standard Sp get 0.024 
25 Per Cent Paved 5, 7 0.021 
50 Per Cent Paved 5 0.018 
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List of Symbols 


A cross-section area of flow, sq ft 

D pipe diameter or equivalent diameter, ft 

f Darcy friction factor 

g acceleration of gravity 32.2 ft/sec/sec 

h mean difference between readings from impact and static tubes 
H, head loss due to pipe friction, ft 

L Length of conduit reach, ft 

n Manning roughness coefficient 

Q rate of flow, cfs 


R hydraulic radius, ft 7 
R, Reynolds number 
Ss Slope of hydraulic grade line 
Vv average velocity of flow, fps 
Vy local velocity of flow, fps 
Vx shear velocity, fps 
v kinematic viscosity, sq ft/sec 
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REVISED COMPUTATION OF A 
VELOCITY HEAD WEIGHTED VALUE 


Joe M. Lara,! M. ASCE and Kenneth B. Schroeder,2 M. ASCE 


ABSTRACT 


revised computational procedure is proposed on the determination of a 
hted value of the velocity head. An example of the procedure is illustrated 
|hydraulic computation of the discharge of a stream using the slope-area 
iod. 


Procedure 


his paper presents a suggested change in the computational procedure 
used by hydrographers to determine the discharge of a stream by the 
2-area method. A simpler way of computing the weighted value of the 
city head in channels conveying overbank or subsectional flow is shown 
xample. 
he usual method for determining the weighted value of velocity head is by 
of a velocity head correction factor, a, computed by the following labori- 
formula: 

3D aS te 3 
: Ki/Ay + Ke/ap * °°" Bo/an 
A= 
: BKM /s (ay? 


fe Ky, Ky, and K, are the conveyances for each of the component sections 
reas Aj, Ao, and An, respectively. 


“Discussion open until February 1, 1960. To extend the closing date one month, a 
itten request must be filed with the Executive Secretary, ASCE, Paper 2149 is 

of the copyrighted Journal of the Hydraulics Division, Proceedings of the Amer- 
an Society of Civil Engineers, Vol. 85, No. HY 9, September, 1959. 


dr. Engr., Commissioner’s Office, Bureau of Reclamation, Denver, 


i oO. 
ir. Engr., Commissioner’s Office, Bureau of Reclamation, Denver, 
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m(K) = Ky, + Ko + +++ Ky 
Z(A) =A, + AQ + +++ Ap 


Conveyance or K is defined in terms of Manning’s formula as: 


x = 2486 ARB 
an 


where R is the hydraulic radius, 
A is the cross-sectional area, 
and on is the roughness coefficient. 


The above method assumes that the velocity within each subsection 1, 2, 0: 
n is constant and that the discharge of each subsection is proportional to its 
relative conveyance. 

The equation for a@ reveals that astronomical values can result in the 
computation of K. As an example, for K = 230,000, K3 = 12,200,000,000,000,¢ 
or 12.2 x 1015, 

The following procedure is proposed for computing the weighted value of 
the velocity head by a more convenient method. The proposed procedure is 
similar in principle to that used by the Corps of Engineers’ in backwater 
computations.(1) This procedure can be best illustrated by an example using 
the slope-area data collected at Sections 2 and 3 of a stream. A table of the 
sample computations and the descriptions of the determination of the value in 
each column are shown.on the following pages: 


Column Number Description 

1. Section designation. Sections 2 and 3 are the 
downstream and upstream cross sections, 
respectively. 

2. Observed water surface elevations at each 
section. 

3. Subdivisions of cross-sectional areas; square. 
feet. 
Distance between cross sections; feet. 

5. Assumed discharge of 10,800 cfs. (trial 
discharge.) 

6. Subdivisional conveyances of each cross 
section. 2 

Hf Subdivisional discharges computed as follows: ; 


Ma a es 


ap = LP) @ 
z Kg 
for example on line 1 


= 175,000 x 10,800 _ 
Qp Se epee 10,520 cfs. 


¥ Ps 
ee ae, ee ee me 
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Column Number 
8. 


10. 


11. 


12. 


13. 


14, 


15. 
16. 


*Refer to Appendix. 
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Subdivisional velocities from the formula: 


Vp = - — Column 7 divided by Column 3 
for line 1 
10,520 
= hat Bande = 
Vp 1,895 5.55 feet per second 


Column 8 squared times the corresponding 
values in Column 7. 


This is the weighted velocity head value com- 
puted by dividing the sum of Column 9 by the 


product of the total discharge and twice the 
gravitational acceleration—illustrated alge- 
braically as follows: 


ea (va) 
Qxeg 
for Section 3 
ie 32', 376 = 0.466 feet. 
x 2 x 32.2 


Difference in velocity heads. The algebraic 
sign should be retained. 


Observed fall of the water surface—the differ- 
ence in the values of Column 2; feet. 


Friction loss computed by the following 
formulas*: 


If Column 11 is negative, as in the example, 
hp = f£ +0. - 
rg a AA et 


of Column 12 plus 0.9 of Column 11. 
If Column 11 is positive, 


or Column 12 plus 0.5 of Column 11. 


Either formula accounts for an eddy loss cor- 
rection. 


Energy gradient computed by getietaa Column 

by Column 4, 

The square root of Column 14. ; 
¥ 

Arithmetic average of the total conveyances of 

Column 6. (Sometimes the geometric mean if 

used. ) 


a aa 
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ylumn Number Description 


iT. Discharge resulting from the product of 
Columns 15 and 16. When this value agrees 
reasonably with the assumed discharge (Column 
5), the computations are completed. If the re- 
sults are too extreme, another discharge is as- 
sumed and a new set of computations is made. 


To compare this suggested procedure with the usual one, the following 
Iculation shows how a would be derived for Section 3 in the example: 


(175,000)3 + (3,670)3 + (935)3 
oy = — (2.895)? (176)2 (82.4) 
(179, 605 )3 
(2,153.4) 


1,490,000,000 + 1,600,000 + 121,000 
1, 250,000,000 


=e .19 


2 : 
This value multiplied by (0,800 7 2,188.5) results in 0.466, the weighted 


locity head for Section 3, the same answer as obtained in the example. In 
dition, another like computation must be made for Section 2. 

The above suggested procedure provides additional information not readily 
ailable from the usual method and that is—the velocities and discharges 
curring in each subdivision are readily computed. True, they may be 
poretical, but they indicate how much of the flow is conveyed by each seg- 
pnt. This information is particularly useful for the computations of back- 
ter profiles. 


| 
| 


APPENDIX 


rnoulli Theorem 
| 


: 
: 
i 


Section Section 
| 2 
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hy, and hy» = velocity heads at Sections 1 and 2, respectively 


Z, and Z) = water surface elevations at Sections 1 and 2, respectively 
he = friction head loss 


hy = eddy loss; considered as a percentage, m, of the change in 
velocity heads = m (hy, - hy, ) 


Writing the energy equation: 
ayy, = Eat By hee 
Substituting and transposing 
Ba = Z, = 2 - -m = 
f 1 P) ¥ By, By, (hy h, a 
Let f = fall in water surface = Z1 - Zo 
ge 2 ti Byye hy, - ™ (hy, - hy) 
For an expanding reach 
> 7% = is positive 
Bee hie lay area 


m = .50 for positive values 
Substituting 


hp = f + ( - - 0.50 3 
one ny) 30 (hy, - by.) 
= £ + 0.50 (hy, - hy) 
For a contracting reach 


Siebert 


m = .10 for negative values 


- by, is negative 


Substituting 
hp = f + (by, - hy) - 0.10 (by, - hy.) 
= f + 0.90 (hy, - ny) 
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EDDY DIFFUSION IN HOMOGENEOUS TURBULENCE 


Gerald T. Orlob,1 A.M. ASCE 


SYNOPSIS 


‘Taylor’s theory of diffusion by continuous movements and the Einstein 
uation of diffusion were applied to eddy diffusion of particles in a two 
mensional field of homogeneous turbulence produced by broad shallow chan- 
. with extreme bottom roughness. A simple method of determining the 
grangian eddy size, eddy diffusivity, turbulence intensity, and mixing length 
Ss devised. Application of this technique provided experimental verification 
the Kolmogoroff similarity principle which relates eddy diffusion, the rate 
energy dissipation, and the Lagrangian eddy scale. 


INTRODUCTION 


An interchange of fluid between two neighboring zones in a turbulent field 
fessarily results in the simultaneous interchange of every characteristic 
the fluid masses involved. For example, if a fluid mass labeled with a 
cer, say fluorescein dye, is transported laterally across an imaginary 
mdary into a zone which is free of contamination, it must, in order to sat- 
7 continuity, be replaced by a fluid mass which contains no fluorescein. 

> immediate result is to increase the average concentration of the dye on 
side of the boundary and reduce it on the other, each region acquiring the 
perties of the other. The process, if allowed to continue indefinitely will, 
rourse, result in producing a uniform concentration of dye throughout both 
; a turbulent field, the masses involved in the interchange are called 

ies” and the process is known as eddy diffusion because of its similarity 
lecular diffusion. The primary difference between the two processes 
in the analogy drawn between discrete particles, or molecules, moving in 


: Discussion open until February 1, 1960. To extend the closing date one month, a 
ritten request must be filed with the Executive Secretary, ASCE. Paper 2150 is 
art of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
nerican Society of Civil Engineers, Vol. 85, No. HY 9, September, 1959. 


st. Prof. of Civ. Eng., Univ. of Calif., Berkeley, Calif. 
E2055 
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free space, on the one hand, and the rather ill-defined lumps of fluid which 
are called eddies, on the other. Nevertheless, the final result, complete mi 
ing of all fluid properties, is the same whether it occurs by molecular or 
eddy diffusion. In a turbulent regime, the effect of molecular diffusion is 
usually so small that it may be neglected. 

The rate of change in concentration during eddy diffusion is determined 
fundamentally by two independent variables: the carrying capacity of the 
transverse flow and the extent of difference between the concentrations on 
either side of the arbitrary boundary. The latter variable can be expressed 
mathematically or evaluated experimentally as the lateral rate of change of 
concentration. The former is defined as the lateral coefficient of eddy dif- 
fusion, Dz, and is a characteristic of the turbulent field in which fluid dy- 
namicists are particularly interested. 

The coefficient of eddy diffusion has been defined in several ways. If 
eddies are treated as discrete entities which behave in a manner analogous 
to particles experiencing Brownian motion, the Einstein(1) equation is ap- 
plicable, 


2 
ewihdiog) 
Ugh Tos sqy ( 


where D, is the lateral coefficient of eddy diffusion, t is time the o, is th 


standard deviation of the lateral displacement from the mean position of a 
statistically significant number of eddies. In a turbulent stream, eddy mo- 
tions may be traced by introducing particles of a foreign substance which ci 
be readily identified and the motions of which are equivalent to those of the 
water mass they displace. | 

Fick’s laws of diffusion, which were derived from analogy with equations 
for heat conduction, have also been commonly applied to eddy diffusion. Fic 
second law, 


saibiae gat id te 
dt z dz? ( 


where c is the concentration of diffusing substance, z is the lateral positic 
and t and Dz are as previously defined, has been especially popular. It — 
should be recognized, however, that when D, is not constant but is a functi 
of t or z, solution of (2) is difficult and in many cases cannot be an 
by formal procedures. 

The mixing length concept of fluid mixing, usually credited to Prandtl, 
recognizes an analogy between fluid turbulence and the “mean free path” ¢ 
; cept of gas dynamics. This results in a definition of eddy diffusion 


LV | 


where £ isthe “mischungsweg” or mixing length V w'2 is the root me 
square velocity fluctuation in the z-direction. This relation for D, posse 


the advantage of incorporation of a property of the flow, the velocity fluctu 
tion. However, the analogy to gaseous diffusion is rather loosely drawn 


ee | ee 
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» result that the mixing length idea has gradually given way to statistical 
striptions of mixing phenomena. 

The modern day statistical theories of turbulence stemmed largely from 
ylor’s theory of diffusion by continuous movements presented in 1921. (2) 
ylor described the diffusion of discrete particles of fluid by means of the 


iation 
t 
pee wel Ree 
SL AAS LS (4) 
ie) 


ere w'@ is the mean square velocity fluctuation and é is the time displace- 
nt between positions of the same fluid particle according to the Lagrangian 
1eme and R¢, the velocity fluctuation correlation coefficient is defined 


i] ( 
eg Wt, “hes 

E ae (5) 

he 
° 

Shortly after Taylor’s theory was presented, Richardson(3) noticed the 
lure of the Fickian laws of diffusion to account for the extremely wide 
ige in values of the diffusion coefficient observed in atmospheric turbulence. 
jintroduced the concept of “neighbor distribution” which for an average 
yaration of particle pairs, L, results in the neighbor diffusion equation 


2: ec bq 
a 5c | Ftu) Se (6) 


sre q is the mean number of neighbors per unit length of L, the neighbor 
aration, and F(L) is a diffusion coefficient which relates to the scale of 
diffusion phenomena. Based on some rather crude observations of atmos- 
iric diffusion and the associated scales, Richardson suggested that 


abebe e pin’ 7) 


hich e was taken as aconstant. This relation, although presented with- 
|basis in theory, has been designated as the “Four-Thirds Law.” It has 
mn applied to oceanographic diffusion by Richardson(4,5) and others(6) with 
ying degrees of success. The exponent four thirds was viewed by some in- 
as a universal number despite the lack of satisfactory experimen- 
ipport and the absence of a theoretical basis for this conclusion. 
he four-thirds relationship remained without prospect of theoretical sup- 
until 1941 when Kolmogoroff(7,8) introduced a somewhat new theory of 
spectrum of turbulence. Kolmogoroff’s development is based on the funda- 
tal concept that eddying motions are characterized by a wide range of 
h scales. Energy of eddy motion is then conceived to be passed down the 
trum through a succession of eddies of smaller size ultimately being dis- 
ted by viscous action as heat. In this process, as the scale is reduced, 
luence of the large eddies is felt less and less so that within the middle 


' 
; 
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range of the spectrum eddies are assumed to possess properties common t 
all types of turbulence. By definition, such eddies are isotropic. 

Kolmogoroff’s analysis is built on two hypotheses, both applicable to mo: 
tions with large Reynolds numbers: 


1) The small-scale components of the motion depend only on the viscosit 
v, and the mean rate of dissipation of energy per unit mass of fluid, E. 


2) There is a sub-range of small eddies in which mean properties depen 
only on the mean rate of dissipation of energy per unit mass of the fluid. 


Batchelor(9,10) in reviewing Kolmogoroff’s theory argued on dimensiona 
grounds that the effect of variation in the parameters E and vp could only b 
to change the effective length and time scales of the motion, thus 


sie) of Es go) ( 


If the customary analogy between molecular and eddy diffusion is used, 1 
can be replaced with D. By dimensional analysis, there results 


ay : 
D = Constant E°L” 


where L may be taken as the mean size of the eddies participating in the d 
fusion process. Turbulence theories proposed independently in 1945 by 
Onsager(11) and in 1948 by von Weizsicker(12) lead to identical relationshi 

It appears that the empirical proposal of Richardson was apparently a 
good one and in agreement with the Kolmogoroff similarity theory; howeve: 
the literature fails to reveal an experimental confirmation of the theory un 
carefully controlled conditions. 


Experimental Techniques 


The essential requirements for experimental investigation of the 
Kolmogoroff principle are: 


1) a test facility which will create a statistically homogeneous turbulen 
field in which the diffusion phenomenon can be studied, and 


2) techniques for evaluation of the scale of turbulence, the rate of energ 
dissipation per unit mass of fluid, and the coefficient of eddy diffusion. 


Statistically homogeneous turbulence in three dimensions is virtually in 
possible to produce in a laboratory. Conditions along the center line of a 
large conduit probably come closest to the ideal requirements, but boundai 
effects severely restrict the scales’of turbulence which may be examined. 
Wind or water tunnel grids may be expected to produce homogeneity of ta 
lence in the field stream but only in planes parallel to the grids. Because 
scale of turbulence diminishes with distance from the grid, this form of tu 
lence generator is not suited to a study which attempts to relate eddy scal 
and eddy diffusion at a constant rate of energy dissipation. 

If the requirement of three dimensional homogeneity is relaxed and a 
dimensional field is substituted, the experimental techniques are consider 
simplified. A broad, shallow, open channel of sufficient length and bottom 
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ighness to minimize end effects meets the requirements of two-dimensional 
mogeneity very satisfactorily provided the diffusion phenomena is con- 
‘ained to a single plane, i.e. the surface of the stream. 

The experimental channel constructed for the studies was 4 feet wide, 8 
shes deep and 38 feet in length (Figure 1) and was supported on leveling 
rews which permitted adjustment to uniform slopes ranging from 0 to 0.01. 
7 One-quarter horsepower propellor pumps were capable of delivering 

ws up to 1000 gpm to the stilling basin at the head of the flume. A by-pass 
ive between the stilling basin and the return flume facilitate discharge 
atrol. 

Bottom roughness for the experimental channel consisted of expanded 

‘tal meshes as shown in Figure 2 which were stapled to the plywood bottom 
the channel. Three mesh sizes designated by the nominal dimensions of 

: minor diagonal of the diamond pattern as 0.5, 1.0 and 1.5 inches, were 
ployed in the investigation. 

‘The requirement that the diffusion phenomenon be observed only in a plane 
ere two-dimensional homogeneity could be achieved dictated the use of 
ating particles, the motions of which were influenced by turbulence near 

: stream surface. Use of discrete particle motions to statistically charac- 
ize the turbulent regime is in accord with the theories of Taylor(2) and 
shardson.(3) Also, a simple direct method for statistical interpretation of 
ticle motion in terms of eddy diffusion is available through Einstein’s 
zation (Equation 1). 

‘The particles employed in the studies were polyethylene discs approxi- 
'tely one-eighth of an inch in diameter and one-sixteenth of an inch thick. 
lyethylene has a specific gravity of approximately 0.975 so that, when fully 
t, discs of this material just barely float. The buoyancy and size of the 

‘es used in these investigations were apparently satisfactory to prevent 

ir being entrained in the flow and to insure that their motion was repre- 
: of the upper few millimeters of the flowing stream. Particles were 
eased singly through a small glass funnel suspended about one-eighth of an 
h over the water surface from a point gage assembly which could be moved 
rails along the channel sides. 

In accordance with diffusion theory, the position of each particle must be 
ed at a certain time following the instant of release. Each point is then 
jected to a single axis. The distribution of particle positions along this 

's can be defined then in terms of the standard deviation from the mean 
sition. The procedure used in the investigations reported herein differed 
htly in that particles were all intercepted from the flow at a fixed dis- 

ice downstream from the point of release. However, since the theoretical 
ern would be symmetrical about this arbitrarily selected axis, the total 
neriences, i.e., the total accumulated path length of all particles, would be 
mtical by either technique. It is, therefore, reasonable to conclude that the 
ard deviation calculated from an array of particles along a fixed line of 
srception a distance, x, from the source is substantially equal to that of a 
ern whose centroid is also x distant from the source. 

Particles were collected from the surface of the experimental channel by 
erting a scoop of ordinary galvanized window screen into the top one- 
irter inch of the flow. In order to simplify the statistical analysis of the 
ticle positions, the particle collector was compartmentalized with thin 

al sheets spaced one centimeter apart across the full width of the channel. 
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By counting the number of particles in each compartment and recording thes 
numbers with respect to the position of the compartment from a fixed refer- 
ence point, the statistical array was defined. 

In order to reduce the time required for determination of the standard 
deviation of particle position, ordinarily a very time consuming operation, a 
simple graphical procedure was adopted. It consisted of accumulating the 
particles captured in the collector from one side of the dispersion pattern, 
expressing them as per cent of total numbers, and plotting on arithmetic 
probability paper at intervals corresponding to the spacing of cells. The re. 
sulting plot was invariably of the typical Gaussian form; that is, it could be 
closely fitted with a straight line as shown in Figure 3. The standard devia- 
tion was then estimated graphically from the slope of the line. A compariso: 
of this technique with the formal computation technique for determination of 
o, indicated the standard deviation in the error of the graphical method was 


less than 1.0 per cent. It was concluded that these errors were of such sma 
order compared to other sources of experimental variation and the time sav 
ing so great that the graphical technique was fully justified. Moreover, the 
method facilitated interpretation of experimental results while the experime: 
was in progress thus avoiding the occasional necessity of discarding incon- 
clusive or incomplete results and thus producing maximum benefit from the 
research effort. 

Prior to the actual investigation, the performance of the experimental 
channel was verified in order to determine whether or not two dimensional 
homogeneity was actually attained. Determinations of Oz at a fixed distance 
from the source over the width and along the length of the channel indicated 
that except for points very close to the channel walls, say less than 2 inches 
from the walls and within 5 or 6 feet of the entrance or exit channel, the tur- 
bulent structure was essentially homogeneous. During all subsequent tests, 
observations were confined to a smaller area about 2 feet wide and 15 feet 
long in the central portion of the channel. 


Results 


Using the methods described above, it was possible to determine the 
standard deviation of the particle dispersion pattern at varying distances 
from the point of release thus defining, experimentally, the relationship be- 
tween o, and x, According to the basic Einstein diffusion equation (Equatior 
1), oz is directly proportional to x1/2 when the coefficient of diffusion and — 
stream velocity are constant. Both the theoretical analyses offered by : 
Taylor(2) and by Richardson(3) agree with this relation for large values of — 
x but suggest that very close to the source, 0, is proportional to x. Thus, ; 


one would expect, if the theories are valid, to find that a dispersion pattern 
determined experimentally would be wedge-shaped near the source and wo 
gradually approach a parabolic pattern as the distance from the source is 
creased (Figure 4), Also, the diffusion Coefficient, which may be considere 
a function of x, would be expected to increase near the source in proportion 
to x, ultimately attaining a constant and maximum value distant from the 
source. 

If the Einstein equation is modified by substitution of dx/a for dt, the — 
diffusion coefficient becomes related to distance by 
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CELL SPACING - CENTIMETERS 


| 
IGURE 3 PARTICLE DISPERSION IN RELATION TO 
| DISTANCE FROM SOURCE 
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pliant Bel oa; d(og) 
Z 2 dx (10) 


where u is the mean stream velocity and the other terms are as previously 


lefined. A plot of o2 versus x thus indicates by its slope the region in 


which Dz may be expected to be constant and, in particular, how oz, varies 
with x (or t) close to the source. Experimental results, such as those shown 
n Figure 4 for a smooth bottom channel, could often be fitted with straight 
ines of the form 


2 
OPEC AITIEX (11). 


n which the exponent n represents the slope of the curve and m is a 
constant. 
Differentiating (11) with respect to x produces the general equation for 
the diffusion coefficient at any distance x, from the source 
- | - n-1 12 
D, x UMnx (12) 
in which the constants m and n characterize the tangent to the o 2x curve 
‘Equation 11). 3 
If o, itself may be regarded as a crude measure of the scale of the dif- 
fusion phenomena, then by replacing x in equation 12 by its equivalent in 
rerms of o, there results 


| n-| 


= 2 Cra p94 
: D umn || (13) 


| 
If n = 1, the special case of constant diffusion coefficient results. If n=2, 
both D, and o, are directly proportional to x, the distance from the source, 
as suggested by Taylor and Richardson for small values of x. And finally, if 
n = 3, there results the special case of the diffusion coefficient proportional 
to the four-thirds power of the diffusion scale (07), i.e. a statement of the 
similarity principle of turbulent diffusion for a constant rate of energy dissi- 
pation per unit mass. 
_ In Figure 4 the value of n is seen to range from 1.0 to 1.59 depending both 
n x and on the Reynolds number, iy/v . There were indications in these 
nd all similar series that at sufficiently large values of x the value of n 
eached 1.0 and that n tended to increase as the source was approached. The 
maximum values of n observed for the polyethylene particle experiments 
seldom exceeded 2.0. The increase in n with Reynolds number suggests that 
e scale and intensity of turbulence also increased so that the region of 
onstant D, was virtually pushed out of the experimental region and could not 
s directly observed. = 
‘It was apparent from preliminary results with floating particles that veri- 
cation of a “four-thirds” law or the similarity principle would not be possi- 
le by direct observation of the dispersion pattern in laboratory scale chan- 
1els. The sub-range of eddy sizes to which the similarity concept would apply 
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must have been much smaller even than the Size of particles employed. To 
further support this conclusion, a limited number of studies with fluorescien 
dye streams were conducted. By injecting the dye at the surface of the 
stream through a fine glass capillary the shape of dispersion pattern much 
closer to the source could be examined. Such dye streams were photographed 
with time exposures on high speed film under ultra violet light. 

Assuming that the dispersion pattern produced on the photographs defined 
limits of + 202 7, it was noted, as shown in Figure 5, that the value of n 
frequently exceeded 2.0 and occasionally approximated the value of 3.0 
theoretically required for conclusive verification of Kolmogoroff’s relation. 
It is clear that direct corroboration in the laboratory of the applicability of 
equation 9 by inspection of individual dispersion patterns will require much 
finer measurements than were possible in these studies. 

However, the relationship between D,, the eddy diffusion coefficient, and 
Lg, the mean eddy size, can be determined through a comparison of the 
properties of several different turbulence fields. In order to accomplish this 
comparison it was necessary to relate characteristics of the dispersion pat- 
tern, viz oz, to the size of eddies which brought about dispersion. By 
equating the basic definitions of eddy diffusion by the Einstein and Taylor 
equations, one obtains : 

2 os 
| d(oz) é 
ary Ww Re d 5 


(14) 


fe) 

In order to integrate the right hand side of (14) it is essential to know the 
form of Ré. Many experimental observations, particularly those of 
Dryden, (13) Laufer,(14) Favre, Gaviglo, and Dumas(15,16) and Ippen and 
. Raichlen(17) show that the R§é -& curve may be closely approximated by an 
equation of the form 


sae (15) 


Substituting into (14), integrating, and replacing t with x/ai, produces the 
general equation 


‘ ig He 
o, = aut | ky 1-8) (16) 


With the aid of equation (16) and an experimental knowledge of o,, x, and i, 
it is possible to determine both k and w'2 as well as those other quantities 
derivable from the basic equation. It will be noted that the shape of the 
OZ" -X curve is determined only by the eddy scale reciprocal, k/u, the velo- 
city fluctuation being only a translational term which positions the curve 
along the vertical axis. Thus, k/i can be determined without knowing w12 

by simply fitting a curve of relative 0,2 versus x to experimental data over — 
a suitable range of x. In fact, one need only know two pairs of values of 0,2 
and x to obtain an estimate of eddy scale. Once k/i has been established, 
w'2 can be calculated from Equation 16. Since the velocity is ordinarily 


known, k can be determined thus permitting calculation of the Lagrangian eddy | 
size 
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t 
-kt T 
= U = UU (17) 
Ce U ee Ot k 
) 
the eddy diffusion coefficient 
ye (18) 
WwW 
(co = — 
D {> k 
and the mixing length 
ae 
f = \w' (19) 
k 


Z 


In order to simplify analysis of some 150 sets of data collected in this in- 
vestigation, a group of normalized 0,2 -x curves were prepared, samples of 


which are shown in Figure 6. The ordinates for each of these curves are the 
values of of relative to the value at x = 1. The abscissas cover the normal 


experimental range from 1 to 20 feet. Each curve is labeled with the para- 
meter k/i. After plotting o 2 -x data for a given run, an overlay carrying 


the normalized curves at the same scale was superimposed and the appropri. 
ate value of k/i' was selected. Figure 7 shows a typical set of experimental 
data fitted with the general curve. As might be expected, the formulated 
curve agrees with the experimental data fairly well within the range of 1 to 
10 feet from the source but departs somewhat in the region close to the 
source. 

While it was not possible to determine the relationship between eddy size 
and eddy diffusion for an individual dispersion pattern it appeared feasible to 
attempt relating the ultimate coefficient D, (0) and the Lagrangian average 
eddy size, Lg obtained from a number of independent experiments. Although 
the average eddy size, characterized by the area under the Ré - § curve, is 
_ not that limiting size of eddy, L,,, which would be expected to control D,(«) 
the two sizes must be closely related so that 


Ets OL, (20 


where a is a constant of proportionality. Thus, Kolmogoroff’s similarity 
principle, stated in terms of Lagrangian average eddy size and the ultimate 
diffusion coefficient, becomes 


4 
Do) = Constant E's ie (2: 


The rate of energy dissipation per unit mass of fluid in a broad open Ghanian 
is defined as 


Ets UgS, (2: 
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where U is the mean stream velocity, g is the acceleration of gravity, and 

Se is the slope of the energy grade line. 

In the experimental channel it was virtually impossible to attain precise 
control over either E or L, since both quantities were affected by changes 
in slope, velocity, and channel roughness. To examine the independent vari- 
ability of E with D,(o) for example, it was necessary to relate Dz (o) and 
Lg empirically for a number of channel slopes and roughnesses, cross-plot 
the resulting curves for constant L,, and compute the associated values of E. 
Figure 8 shows typical examples of Dz (o) - La curves of the type usually 
associated with the Richardson relation (Equation 7) in which E is not held 
constant. It is interesting to note that each set of plotted points can be fairly 
well represented with a line of four-thirds slope, an apparent corroboration 
of the so-called “Four-Thirds Law.” Also, it is of some significance to ob- 
serve that the vertical position of the curves is affected Significantly by the 
slope of the channel bottom but apparently is little influenced by the size of 
bottom roughness. The fact that Dz (o) is dependent on slope, i.e. is propor- 
tional to the rate of energy dissipation, generally agrees with the similarity 
concept. 

The computation of E for constant values of Lg also required empirically 
relating the mean channel velocity, U, to Lg and in turn selecting Sg from 
a channel calibration curve of Se versus channel slope, and channel velocity. 
‘Values of E determined from equation 22 and corresponding to eddy scales 
of 0.2, 0.5, and 1.0 feet were plotted against the appropriate values of Dz («) 
from Figure 8 to yield Figure 9. 

Despite the fact that each plotted point in Figure 9 is supported by 50 to 
100 separate determinations of the standard deviation in the dispersion pat- 
tern there is considerable scatter in the plots. This must be attributed, in 
part, to the accuracy with which the mean velocity and the true slope of the 
energy line could be determined and the inherent variations in the statistical 
methods for evaluation of L,. Nevertheless, the general trend established by 
the plots appears to agree fairly well with that expected from theory. 

Figure 10 shows the coefficient of lateral 7 diffusion as a function of 

[Lagrangian eddy size at constant E (0.05 ft2/sec%). The trend line in the 

figure is seen to have a shape similar to those of Figure 8 approximating a 
ur-thirds slope at large eddy sizes and tending to flatten at small sizes. 
Again, the experimental data appear to support the theory. 

A further test of the validity of the similarity principle is obtained by 
lotting values of E1/ 31,4/ 3 against D, without regard to the slope of the 
aannel bottom or the condition of bottom roughness. Figure 11 summarizes 

the data representing 105 separate experiments in which-U, Se, and bottom 

ughness are all independent variables. It is especially significant to note 
nat the plotted points can be fitted with a straight line on a unit slope. The 
uation of the line which best fits the data of these experiments is 


: 
y ¥, 2 
Djo) = 0.0136 E° ae a 


where Dz, is the coefficient of the lateral eddy diffusion, E is the rate of 
nergy dissipation per unit mass, and L, is the Lagrangian eddy size. The 

tant is taken as representative of all degrees of bottom roughness pro- 
luced by the three roughness meshes employed in the experiments. 
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The fact that all the data from the series used to prepare Figure 11 are 
rather equally distributed about the best fit line indicates that the scatter is 
probably due to experimental errors alone. It may also be concluded that the 
independent variables are probably accounted for in the quantity E1 3144 3 
or that the experimental range of the variable was not sufficient to cause a 
separation in the plotted points. The latter conclusion is probably applicable 
in the case of bottom roughness, the absolute value of which was not greatly 
affected by a change in mesh size. There is some indication from exploratory 
tests with a smooth channel bottom that the vertical position of the line of best 
fit is a function of bottom roughness; however, data are insufficient to define 
this effect. 


Discussion 


An investigation of the diffusion properties of a turbulent stream requires 
the definition of the coefficient of eddy diffusion and the development of a tech- 
nique satisfactory for its determination. There are many ways in which Dz 
can be defined but most of these require the separate evaluation of quantities 
such as the turbulence intensity, the mixing length, or the eddy scale in order 
that actual values of the eddy diffusion coefficient may be calculated. Be- 
cause of difficulties involved in the physical measurement of these quantities 
there has been a tendency for investigators to resort to empirical or 
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semi-empirical *laws” which relate two properties of the flow, usually Dz 
and one other characteristic, through constants of proportionality and appro- 
priate exponents. Often the empirical equation is only a convenient device 
useful for rough estimation but bearing no resemblance in form or content to 
the true theoretical relationship. Occasionally the researcher who has pro- 
ceeded with only a partially developed theory can, with the aid of considerable 
intuition and astute observation, hit upon a relationship which comes close to 
the truth but which is not quite complete. This appears to have been the case 
with the semi-empirical relationship first proposed by Richardson and which 
has been referred to as the “Four-Thirds Law.” 

Richardson observed that the coefficient of eddy diffusion was nearly pro- 
portional to the four-thirds power of some measure of the scale of the phe- 
nomenon causing diffusion. He proposed a theory of “neighbor diffusion” 
which involved the statistical averaging of “neighbor separations” of adjacent 
particles. The mean neighbor separation he took to be a measure of “scale.” 
Unfortunately, Richardson’s neighbor diffusion theory is difficult to apply, and 
except for a few studies on ocean surfaces, it has been little used. 

Pearson(6,18) has summarized the relationship between the coefficient of 
eddy diffusion and the scale of the diffusion phenomenon by reporting the ob- 
servations of many independent investigators, including himself, for a wide 
variety of scales. Figure 12 is a typical plot of the type Pearson prepared, 
including all of the data reported by him together with some observations 
collected during the investigations which are reported herein. It is possible 
to see from the figure why one might be inclined to regard Richardson’s re- 
lation as a universal law. A line of four-thirds slope fairly well establishes 
the over-all trend of data from scales of about 0.1 foot to scales of more than 
1000 feet. The considerable scatter in the separate observations can probably 
be attributed, in part, to the variety of definitions used for “scale.” These 
definitions included mean neighbor separation of particles, one-quarter of the 
width of the dispersion pattern (approximately = o7), the Lagrangian eddy 
size, and some arbitrary dimensions such as width or depth of channel. 

It has been shown theoretically according to the Similarity concept of 
Kolmogoroff that the fundamental relation should include a third variable, 
the rate of energy dissipation per unit mass. It may be argued that much of 
the spread in the observations of Figure 12 could have been eliminated had 
the rate of energy dissipation been considered in each case. The most ap- 
propriate abscissa for a plot of a supposed universal law would be E1/31,,4/3 
with a consistent definition for D, E, and L as in Figure 11. 

It is also apparent from consideration of the shape of the dispersion- 
distance curve that the coefficient of eddy diffusion is not a quantity of un- 
limited magnitude whose value is controlled solely by the scale of the phenom- 
enon, but that it is limited by the maximum size of eddy motions which may 
exist in the turbulent regime. The size of these limiting eddies is not only 
controlled by the boundaries of the field but by the intensity of turbulence, 
the supply of energy to the system, the rate of energy dissipation, and proper- 
ties of the fluid, If the maximum eddy sizes are sufficiently large, then, in 
accordance with the similarity principle, there will be a sub-range of eddy ~ 
sizes to which the general diffusion law will apply. Within this sub-range it 
may be expected that the coefficient of eddy diffusion will be proportional to 
the four-thirds power of the eddy size, but on the upper side of this range 
diffusion will become decreasingly dependent on scale until it ultimately 
reaches its maximum constant value. 
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The significance of the concept of a limiting diffusion coefficient can best 
be seen by considering the problem of waste assimilation by a receiving water 
say a large coastal bay. If the waste is discharged at a point on the surface, 
the pattern of dispersion may at first appear to increase rapidly, as in 
Figure 13, fanning outward as the waste diffuses laterally while being carried 
away by the current. In the region near the source the similarity principle 
may apply with D, proportional to the eddy size. When the dispersion pat- 
tern expands through some intermediate range of eddy sizes the principle 
will no longer be applicable. D, may continue to be influenced by eddy size 
but to a lesser degree, and the dispersion curve will gradually change shape 
to become concave toward the axis of flow. Ultimately, the maximum eddy 
size will have become involved in the dispersion pattern and the eddy diffusion 
coefficient will cease to increase. Beyond this point the dispersion curve will 
be parabolic in shape, indicative of a constant D. 


SUMMARY AND CONCLUSIONS 


The Kolmogoroff similarity concept leads to the general conclusion that 
the coefficient of eddy diffusion in homogeneous turbulent flow is proportional 
to the one-third power of the rate of energy dissipation per unit mass and the 
four-thirds power of the scale of eddies participating in diffusion. Using a 
broad open channel with a variety of bottom roughness to produce a two- 
dimensional homogeneous turbulence, the experimental validity of the concept 
was examined. The following conclusions resulted from the investigation: 


1, In laboratory turbulence the sub-range of eddies to which the concept 
can be directly applied has a very small maximum size. The upper size 
limit in the open channel studies reported was probably less than one centi- 
meter. Examination of particle and dye dispersion patterns very close to the 
source indicated the existance of the sub-range. ; 


2. The shape of the dispersion pattern for scales larger than the upper 
limit of the critical sub-range conform closely to the theories of Taylor and 
Richardson. Close to the source the relationship between the standard devia- 
tion of the dispersion pattern and the distance from the source is almost 
linear. Far from the source the dispersion pattern becomes parabolic in 
shape. In the transition range the pattern can be described statistically by 


assuming an exponential form for Taylor’s Lagrangian velocity correlation 
coefficient. 


3. For each turbulent spectrum there is a maximum value for the eddy 
diffusion coefficient which is associated with the limiting eddy size. The 
limiting eddy size is directly proportional to the statistical mean eddy size. 

4. The Kolmogoroff similarity principle adequately describes the relation- 
ship between the ultimate lateral eddy diffusion coefficient Dz, (), the F 


Lagrangian mean eddy size L,, and the rate of energy dissipation per unit { 


mass E, in statistically homogeneous flow in broad open channels. The y 
equation 


y ¥ 
D:ic0}.=. 0, State 
i O36 E* LO 


characterizes 105 turbulence conditions investigated. 
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SOME EXPERIMENTS WITH EMERGENCY SIPHON SPILLWAYS# 
Closure by Warren B. McBirney 


WARREN B. McBIRNEY,1 M. ASCE.—The writer is indebted to Messrs. 
Fred W. Blaisdell and Harold W. Humphreys for their interest in preparing a 
Jiscussion and in offering for consideration the drop inlet type of structure. 

Since much of the discussion centers around the definition of a siphon, it 
might be well to analyze the two definitions given. The first, to which the 
writer subscribes, defines a siphon “as a spillway which requires that the 
liquid be lifted above the headpool level.” This definition is supported by 
Russell,(1) Addison, (2) and Rouse.(3) The Bureau of Reclamation so classi- 
fies such a structure, though this position admittedly is based on applied con- 
siderations and not theoretical ones. 

The second definition, which the authors suggest as being appropriate, is 
“a closed conduit, a portion of which lies above the hydraulic grade line.” 
Three Bureau structures are named as qualifying as siphons, the Heart Butte 
and Shadehill morning-glory spillways in the Dakotas, and the Hungry Horse 
morning-glory in Montana. The writer believes none of these operate as 
siphons, even by the authors’ definition. All three spillways were studied ex- 
tensively by model experiments. In general it can be said that certain por- 
tions of the profiles showed negative pressures for the unsubmerged flow con- 
Hitions, but those pressures prevailed in areas below the crest sections which 
received radial or near-radial flow, and not between certain elevations for the 
full circumference. Thus, the low pressures were similar to those on an ogee 
shape, and were related to tendencies for flow separation from the boundaries 
because of imperfect profiles. Since the upper nappe surfaces were essen- 

ially at atmospheric pressure no siphon head could exist because the flow 

was of the open channel type. 

The possibilities of extremely low pressures occurring with submerged 

flow conditions were indeed recognized. At Heart Butte a combination of nine 
r inlets and deflectors was installed 11.5 feet below the crest and an aerated 

: ts deflector placed on the crown side of the shaft 31.5 feet below the crest. 

rther, the tunnel diameter diverged from 11 to 14 feet in diameter around 

- bend. Pressure tests at submerged flow showed negative pressures no 
eater than about -5 feet of water in certain areas, with no indication that 

Ow pressures existed around the entire circumference, as would have been 

he case if a closed conduit siphon head had been acting. 

There was less opportunity for low pressures in the Shadehill spillway be- 

ause of the large spillway and tunnel diameters and the short drop of 33 feet 

4 / 

. Proc. Paper 1807, October, 1958, by Warren B. McBirney. 

. Hydr. Research Engr., Div. of Eng. Laboratories, Bureau of Reclamation, 


_ Denver, Colo. 
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to the invert of the bend. Negative pressures not exceeding -3 feet of water 
occurred for weir-flow discharges less than 60 percent of design capacity, 
and all measured pressures were well above atmospheric at capacity flow. 
Even so, a flow deflector was installed on the crown side of the shaft 14 feet 
below the crest to insure open channel flow from the deflector downstream to 
the tunnel outlet. 

The Hungry Horse morning-glory was studied in one of the most thorough 
investigations ever conducted by the Bureau. The negative pressures on the 
profile were deliberately created to increase the discharge coefficient and 
were restricted to areas below portions of the crest receiving radial or near- 
radial flow. Other negative pressures were relieved by an air vent. An air 
intake structure of 36 square feet cross-section was installed on the crown 
side of the shaft 31 feet below the crest. This conduit had sufficient air 
capacity to prevent the shaft and tunnel from flowing full for any flow or crest 
head thought possible at the site. The writer as well as the research engi- 
neers responsible for these tests believe there is no possibility of the Hungry 
Horse spillway and tunnel ever flowing full. 

None of these structures in model or prototype has shown any evidence of 
siphonic action in the sense that a full conduit was acting to influence the rate 
of flow at higher elevations in the shaft or tunnel. No priming or depriming 
aids were built since the structures were not envisioned as siphons and every 
effort to prevent siphon action was made. The air demand of each structure 
was more than satisfied with oversized conduits to prevent unsteady rates of 
intake. This is not the type of partialization desired by the writer. 

With the above description of the operation of these spillways completed, 
the writer is now prepared to admit that a morning-glory or drop inlet spill- 
way flowing submerged with an unaerated shaft of sufficient drop in elevation 
will exhibit flow characteristics akin to those in the downstream leg of a sad- 
dle or volute siphon. It is doubtful, however, whether the engineer who pro- 
duces such a cavitation structure would derive much credit for his design. 

The head-discharge curve in Figure A of the discussion, for the drop inlet 
closed conduit spillway, is very similar to a typical morning-glory head- 
discharge curve. Placing an antivortex flat plate over an inlet 3 inches in 
diameter may appear to offer control of the vortex and a uniform proportion- 
ing of air and water. Trial of a similar scheme in a model of the Trinity Dam 
morning-glory, northern California, model size 24 inches in diameter at the 
crest, brought near-disastrous results to the clear plastic of the shaft. Air 
was sucked under the plate intermittently instead of steadily, pressures in the 
shaft varied rapidly between -10 and -30 feet of water (prototype), the forebay 
supplying the structure was soon set in rhythmic motion, and the model vi- 
brated violently. Such action in the prototype would cause serious damage in ~ 
a short time. | 

The writer suggests that the authors employ a high frequency response 
pressure cell and oscillograph to a piezometer centered in an antivortex 
plate over a drop inlet spillway. Lack of visible fluctuations in the head pool — 
for a 3-inch-diameter inlet may not be. adequate proof of stability in the : 
setup. a 
Near the end of their discussion the authors applied a linear multiplication 
factor to compare operation of the Bureau saddle siphons with the drop inlet — 
spillways. Such a comparison by this means is not believed to be valid for 
the following reasons: } 
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1. Even assuming both types of structures to be siphons, scale relation- 
ships for siphons are still very much in doubt and are probably not linear. 


2. The authors compare by the multiplication factor a head on the drop in- 
let which *caused the Illinois spillway to prime” to the heads for maximum 
capacity flows in the standard and proposed designs. The latter heads were 
reached long after both siphons began priming. Actually, the standard siphon 
started to prime with about 0.13 foot of head and the proposed siphon re- 
quired about 0.015 foot of head. Both heads are less than the 0.17 foot shown 
by the authors for the Illinois spillway. 


Neither the authors in their discussion nor the writer in his paper offered 
any information on the effects of changing the rates of rise in the forebays of 
respective spillways. The Bureau test setup was not flexible enough to vary 
the rate of rise for any discharge, but other experimenters in the world have 
been able to do so. Findings are that the lower the rate of rise for a given 
discharge the lower will be the head required to cause the siphon to prime. 
Thus, it is quite probable that the authors’ and writer’s spillways will prime 
at lower heads for lower rates of rise. 

The authors state that extending skirts from the antivortex plate downward 
below the weir crest modifies the structure to the form of a vortex siphon 
spillway and that the head-discharge curve then obtained is similar to those 
shown for the standard and proposed siphons. This would seem to be excel- 
lent evidence that there are fundamental differences in the two types of struc- 
cures and that only one of them should be classified as a siphon. 

The writer does not doubt that the drop inlet closed conduit spillway has 
many uses, but in its current form it does not appear to fit requirements of 
the Bureau of Reclamation for an emergency spillway. 
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EVAPORATION OF LAKE ONTARIO 
Discussion by M. A. Kohler 


MAX A. KOHLER, ! F. ASCE.—Major Hunt states that there are several 
methods of computing evaporation—mass transfer, water budget and energy 
budget—only the first two of which he has applied to Lake Ontario. In addition 
to the methods cited by the author, there is the evaporation-pan approach as 
well as several others which constitute combinations of the four. It should 
also be emphasized that there exists a multitude of mass-transfer evaporation 
equations which have been derived by numerous investigators. 

Considering, first, the mass-transfer approach, it will be noted that the 
author’s equation (Eq. 7) relates evaporation to the vertical gradient of hori- 
zontal wind speed over the lake and that the empirical constant was derived 
from data collected at Lake Hefner, Okla.(1) Although one might expect to ob- 
tain satisfactory estimates of normal monthly wind speed over the water from 
shore anemometers, corresponding estimates of the vertical gradients are 
likely to be appreciably in error. In the Weather Bureau’s analysis of the 
‘Lake Hefner data,(2) several mass-transfer equations were derived empiri- 
‘cally which require wind movement instead of wind gradient. One of the 
equations is 

E = 0.00304 U4 (eg = e€9) (A) 


where E is the lake evaporation in inches per day, Uy is the 4-meter wind 
‘movement in miles per day, eg is the saturation vapor pressure at the tem- 
perature of the water fe eee of mercury), and eg is the 2-meter 
vapor pressure over the water, also inches of mercury. The constant of this 
equation is based on records for 115 days when the water-budget data were 
judged to be accurate within 5 percent. The correlation coefficient is 0.956— 
‘much higher than the author cites for his equation. 
Equation (A) has been applied to Lake Ontario using the author’s data 
(Table 7) with the results shown in Table A. It will be noted that Eq. (A) 
provides a better fit with the water budget than does the author’s Eq. (7), on 
a monthly basis as well as for the year as a whole. 
‘The Weather Bureau(3) has recently completed a series of United States 
maps showing: 


1. Mean annual free-water evaporation 

2. Mean annual class A pan evaporation 

3. Annual class A pan coefficient (ratio of free-water to pan evaporation). 
_ 4. Percent of annual class A pan evaporation occuring during the summer 


(May -Oct). 


a. Proc. Paper 1934, February, 1959, by Ira A. Hunt, Jr. 
4 Chf. Research Hydrologist, U. S. Weather Bureau, Washington, D. C. 
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Table A. LAKE ONTARIO EVAPORATION 


b 
a Vapor pressure Pris 
Month Wind (U,,) Pepe iar gs Evaporation in f 
re hs (eo 3 2) Water 
inches of Hg Eq.(A) Eq. 

an 02 O01 12 TEs! 19 
Pee 389 -OTT 21 19 -10 
Mar 307 -073 .18 16 .06 
Apr 284 .039 .08 209 Ok 
May 192 2001 00 00 Ok 
Jun 163 2013 02 02 00 
Jul 181 ease 19 18 eel 
Aug 190 An yf 23 22 -28 
Sep 220 191 ABE: 27 32 
Oct 231 .150 SPL -23 38 
Nov 292 100 222 20 545) 
Dec 395 .076 -2k voll 34 
Annual eece 2.00 ese 


a. At Lake Hefner, 4-meter wind averaged 0.90 that at 8-meters. This 
factor was applied to the data in the author's Table 7 and converted 
to miles per day. 


b. The author's vapor-pressure difference values are incorrectly labeled 
as "millibars", Strictly speaking, Eq. (A) requires the 2-meter 
vapor pressure over the water rather than the average for several 
stations around the lake, 


5. Standard deviation of annual class A pan evaporation. 


The computations(4) were based on records for the period 1946-55. In the 
course of deriving the free-water map, evaporation was computed from solar 
radiation (or percent sunshine), dewpoint, air temperature and wind for six 
meteorological stations adjacent to Lake Ontario. The computed values 
(ft/yr) are as follows: 


Toronto 2.29 Oswego 2.36 
Buffalo 2.37 Syracuse 2.50 
Rochester 2.56 Canton 2.19 


It would appear from these results that the 10-year mean evaporation from ~ 
the lake is near two and one-quarter feet—a value in close agreement with — 
our mass-transfer and the author’s water-budget results. 

Having essentially substantiated the author’s water-budget value (2.21 
ft/yr) by two independent approaches, it is indeed unfortunate that a question 
must be raised concerning the value used for lake precipitation. Early resul 
from the island network on Lake Michigan have led the author to conclude tha 
annual precipitation over Lake Ontario is only 2.30 ft. as compared to 2.92 ft 
on the tributary drainage. The observational program is still in progress on 
Lake Michigan and final conclusions must await further analysis. Neverthe- 
less, present indications are not entirely in accord with the author’s interpre 
tation. : 

Table B shows record-mean summer, winter and annual precipitation for 
the island stations and comparative data for a group of lake-shore stations. — 
There can be little doubt that the islands (and the lake) receive less summer 
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precipitation that would be indicated by an average of the shore stations. But 
is this deficiency compensated by an excess during the winter months? The 

records from South Manitou (and possibly Beaver Light Station) indicate this 
is the case. 

The author concludes that at least part of the differences among the island 
data result from the fact that the larger islands produce an increase of pre- 
cipitation. It would seem that any heating effect should be greatest at the 
downwind side of a small island (and likely even beyond) whereas St. James 
and the Light Station are on the northern and southern portions of Beaver 
Island, respectively. Some of those concerned with the project have felt for 
some time that virtually all the differences could be explained by relative 
windiness at the gage sites. 

In July 1958, anemometers were installed at four of the island stations 
(gage-orifice height) which were read in the latter part of the following 
October. Mean wind speeds (mph) for the period were as follows: 


St. James, Beaver Island 4.04 
Beaver Island Light Station 2.05 
Tle Aux Galets 7.02 
Shoe Island 7.61 


Without going into details in this discussion, examination of the wind data and 
photographs of the gage exposures show conclusively that the relative catch 
of the island gages is highly correlated with windiness. 


Table B. LAKE MICHIGAN PRECIPITATION, 
(SIX-YEAR MEANS, 1952-1958), IN INCHES 


oe ee ee eee 
STATION Summer Winter Annual 
pen iee Qe ee 


St. James, Beaver Island 13.19 12.87 26.06 
Beaver Island Light Station 13.34 14.91 28.25 
South Manitou Island 13.68 15.04 28.72 
North Manitou Island 12.46 14.56 27.02 
: South Fox Island 13.39 12.21 25.60 
|) Guia stand 13.02 13.29 26.31 
‘Ile Aux Galets 10.79 9.04 19.83 
| ery 11.80 8.14 19.94 
shore stations (weighted)” 14.87 13.86 28.73 


ener en EEE EEE 


: Two years (1956-58) for which the shore-station precipitation averaged 
14.96, 13.22 and 28.18 for the summer, winter and annual periods, 


respectively. 


mate gages in the group are shielded. Unshielded gages at the same 
sites receive a few percent less precipitation, on the average. 


112 September, 1959 HY 


In the light of Wilson’s(5) studies of wind vs. gage catch, it would now appear 
that annual precipitation over Lake Michigan may not be significantly less 
than that observed at a network of shore stations. 

Assuming that the average annual precipitation observed at a network of 
eight stations along the shores of Lake Ontario is equivalent to that over the 
lake, the value in Column 6 of the author’s Table 1 becomes 2.76 ft. The 
water-budget evaporation then becomes 2.67 ft. instead of 2.21 ft. Since this 
revised value is appreciably higher than that derived by any of the other 
methods, the writer has attempted to check the results for a different period 
of record (1946-55) as follows: 


Precipitation (average of 8 shore stations) 2.84 feet 
Runoff (104% of 1934-53 as indicated by average 

of 4 tributary basins) 5.04 
Inflow (Niagara River at Buffalo, N. Y., plus 

Welland Canal) 32.69 
Outflow (St. Lawrence River at Ogdensburg, N. Y.) 38.63 
Change in stage (dropped 1.05 ft. in 10 years) 0.11 
Computed evaporation (1946-55) 2.05 


It is entirely possible that the evaporation was less in 1946-55 than for th 
longer period 1934-53, but much of the difference (2.05 and 2.67 ft.) is beliey 
to be erroneous. Without belaboring possible sources of errors in water- 
budget computations, it should be pointed out that bank-storage effects would 
be of proper sign with respect to the two periods—the lake level dropped 0.1: 
ft. per year during the period 1946-55 and rose 0.11 ft. per year during 1934 
53. Even so, any realistic assessment of the reliability of the water -budget 
for this particular lake will show that an error of one-half foot in the annual 
evaporation is not unlikely. A bias of only one percent in the inflow and out- 
flow items amounts to about 0.7 ft. over the lake. 
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PROBLEMS CONCERNING USE OF LOW HEAD RADIAL GATES ® 
Discussion by J. R. Bowman 


dic. Ft BOWMAN, ! M. ASCE.—Mr. Rhone’s paper constitutes a significant 
contribution toward a better understanding of the hydraulic characteristics 
of tainter (radial) gates. The author and the USBR should be congratulated 
for sharing with the profession the knowledge gained through research sup- 
ported by public funds; in many instances, the results of basic research car- 
ried out by public agencies has been imprudently labeled “Restricted” or “Not 
for General Distribution.” 


Discharge Characteristics 


Quite often it is desirable to derive accurate rating curves for discharge 
under a tainter gate, without recourse to a model test program specifically 
organized for this purpose. This is often the case for smaller regulating 
works whose total costs do not warrant expenditures for model tests. Ina 
few instances, it has been possible to check or modify analytical predictions 
by means of field measurements; more often than not, however, comprehen- 
sive field tests will not be feasible, with the result that the responsibility for 
accurate discharge predictions will rest solely with the designer. 

The problem of how best to determine part gate discharge probably has 
been encountered as long as tainter gates have been in use. Where the gate 
sill is located at the weir crest, the rating curves for the gate will be simi- 
lar to those of a vertical-lift gate, provided that the centerline of trunnions 
is not appreciably higher than the centerline of the gate. Where the trunnions 
ie located at the elevation of headwater, or higher, the approach angle 0 
exerts considerable influence on the value of the discharge coefficient, as 
shown in Figure 3. In cases in which the gate sill is located downstream 
i the weir crest, a shifting control is encountered; The control point 
moves gradually from the gate sill, at small openings, to the weir crest as 
full gate opening is approached. Cases in which the gate sill is located up- 
ream from the crest are not often encountered, inasmuch as such a gate 
i hen usually results in increased costs of associated structures which may 
ore than offset any saving in gate costs. 
In the development of rating curves for a tainter gate whose sill is located 
the weir crest, the writer, as well as many of his predecessors, has 
tilized a method whose validity is largely confirmed by the data presented 
n Figure 4. When the approach angle @ lies generally between 75 and 105 
sgrees, the free discharge coefficient can be applied to all gate openings, 


“Proc. Paper 1935, February, 1959, by Thomas J. Rhone. 
. Civ. Engr., Erik Floor and Assocs., Inc., Chicago, Il. 


114 September, 1959 HY 


with an average error of less than 5 percent. The writer is acquainted with | 
few such cases in which field measurements and predicted values differed by 
not over 5 percent. 

Although no simple and accurate method has yet been devised for comput- 
ing the discharge under a tainter gate whose sill is located downstream from 
the weir crest, Mr. Rhone has presented information which may well serve 
as the basis for the effective utilization of a combination of methods. 

On the basis of analyses of a number of rating curves based on model 
studies, the writer can make the following generalizations: 


1) For gate openings ranging from near zero (gate closed) up to 15 or 20 
percent of full gate, with discharge coefficients related to the head on the 
weir crest, the application of the modified weir formula (first equation) re- 
sults in discharge coefficients that vary from unrealistically high values at 
near zero opening to values nearly equal to the free discharge coefficient at 
somewhat larger openings. (“Full gate” is defined herein as the gate opening 
at which the upper nappe breaks clear of the bottom edge of the gate.) 


2) For the same range in gate openings, but with discharge related to the 
head on the gate sill, application of the orifice formula results ina nearly 
constant range in discharge coefficients. This observation is substantiated 
by the curves of Figure 6. 


3) For larger gate openings, the discharge coefficient, expressed in term 
of the head on the weir crest, is nearly constant and can be closely related tc 
the free discharge coefficient. 


In cases in which extreme accuracy has not been required, the writer has 
made use of an approximation which combines both analytical and graphical 
methods, and which expresses discharge as a function of gate opening; the 
procedure, for a given headwater condition, is as follows: 


1) Determine the depth of water on the gate sill at free discharge, con- 
sidering this depth as “full gate.” For the majority of downstream settings 
on ogee crests, it can be assumed that critical depth will occur in the im- 
mediate vicinity of the gate sill. Any corrections in depth or discharge, due 
to end contraction, velocity of approach, etc., should be made at this time. 


2) Compute and plot discharge versus gate opening from zero opening to 
approximately 25 percent of “full gate” by means of the orifice equation and 
the discharge coefficients of Figure 6. 


3) Draw a straight line, passing through the full gate point of step (1) and 
tangent to the curve of step (2). This completes the rating curve for a given 
headwater condition. a 


In conjunction with step (2) above, it should be stated that, although the 
definition of gate opening shown in Figure 6 represents good academic prac- 
tice, a good deal of manipulation is required for the precise determination of 
any gate opening. It is suggested that computations could be simplified by _ 
expressing gate opening as a vertical distance above the gate sill, and by ad- 
justing discharge coefficients accordingly; any differences in gate settings — 
and/or crest shape then would be accounted for by the approach angle 6, 
which could be determined graphically with sufficient accuracy for any 
purpose. 
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In reference to the author’s apparent bewilderment concerning the results 
own in Figure 7, it should be apparent that, for gate sills located down- 
tream from the crest, the approach angle and, consequently, the coefficient 
f discharge should be less than for a weir crest setting. This trend is amply 
llustrated by the author’s Figure 6. 

The author’s dream of large models set up for the sole purpose of re- 
earch in the field of gate discharge relationships is heartily shared by the 
yriter, whose interests in such matters, unfortunately, are not supported by 
he means for carrying out such a program. It should be of interest to learn 
yhether the Bureau contemplates such a comprehensive program in basic 
esearch. 


Pressure Characteristics 


The results of the pressure profile studies are very well explained by the 
uthor; as such, they incite very little comment on the part of the writer. He 
7ould like to be informed as to the basis for derivation of the profile used as 

datum in the series of tests. In passing, it is of interest to note that the 
pillway at American Falls, which was designed a good many years ago and 
rithout the benefits of present-day knowledge, is in no imminent danger of 
ailure due to negative pressures on the ogee. 


Gate Seals 


It is obvious that this portion of the paper deals primarily with seals that 
re currently being developed or applied by the Bureau of Reclamation. It is 
ot improbable that there are as many different preferred combinations of 
eal types and methods of application as there are gate designers. However, 
nerely describing several popular combinations would likely exceed the 
jmitations placed on the length of a paper. 

_ Although it is universally conceded that low head tainter gates should be 
quipped with side seals, it does not necessarily hold that they must also be 
rovided with bottom seals attached to the gates themselves. In a good many 
arly installations, the bottom edge of the skin plate was arranged to bear on 
¢reosoted timber, which was bolted into the supporting structure so that its 
'p face was flush with the canal floor or weir crest. At Dix Dam in Kentucky, 
srtical-lift spillway gates sealing against timber sills were installed in the 
arly 1930’s; the original timbers continue to provide practically a drop- 
ght seal after more than 25 years of service. In a number of fairly recent 
\Stallations, the skin plate rests on a soft metal, such as babbitt or lead, 
hich has been melted and poured into a structural channel, or other suitable 
ember, embedded in the floor or weir surface. Such installations have not 
ways met with the anticipated degree of success; their shortcomings gen- 
rally stem from one of two sources: 


1) In cases concerning small gates, the weight of the gate does not provide 
equate indentation of the soft metal, with the result that a uniform degree of 
ntact is not achieved. 


2) In the case of a rather deep gate, the hydrostatic pressure at the sill of 
closed gate induces a wire-drawing action; this results in pinhole or sheet 
ts under the skin plate. pan 
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Where rubber bottom seals are necessary, the writer ordinarily does not 
subscribe to the practice of placing the seal on the downstream face of the 
skin plate; he prefers to place the seal on the upstream face, so that hydro- 
static pressure can assist in effecting final closure. In his opinion, the only 
practical difficulty that might arise from the latter practice is that an up- 
stream seal presents a local obstacle to streamline flow and to direct 
mathematical treatment. 

In the application of side seals, the writer is strongly in agreement with 
the author in that a nominal amount of seal preset is necessary to assure 
positive contact between the seal and the embedded seal plate. If no preset 
were provided, a number of factors could reduce (and have been known to 
reduce) the effectiveness of the seal: 


1) Inaccuracies in field alignment of seal plates and/or seals can result in 
excessive out-of-true humps and sags. 


2) A certain amount of lateral play accompanies the vertical movement of 
tainter gates; it may be limited to a fixed amount by rollers or skids, or toa 
variable amount by limiting forces developed within the seals themselves. 


3) Field observations indicate that the contact face of a seal acquires som 
permanent set, resulting from age and usage; the softer the contact face or 
the more flexible the seal bulb (as for “J-seals”), the greater the amount of 
permanent set. 


Seal preset does, of course, introduce frictional forces in addition to the 
amount attributable to hydrostatic pressure alone. This has led many to con- 
sider the possible application of retractible side seals, wherein all frictional 
forces would be absent while the gate is being raised or lowered; full hydro- 
static pressure would be utilized to force the seals into contact with the seal 
plates only after the gate has reached its closed position. Prominent among 
those individuals was the late L. F. Harza; several years ago, he conducted j 
a series of model studies devoted to exploring the possible application of re- 
tractible seals to tainter crest gates. A significant observation derived from 
his studies is the fact that the retractible body of the seal must have con- 
siderable inherent stiffness, so that the activating water pressure at lower 
levels can effect complete closure at headwater level, where hydrostatic 
pressure equals zero. Moreover, the required stiffness would tend to reduce 
the flexibility along the contact face to an extent that could require extensive 
modification in order to regain a portion of the contact flexibility originally — 
lost. , 

The Bureau’s angle seal (or “L-seal”) appears to represent a compromisé 
between the two most popular types of “J-seal” mountings. The L-seal incor 
porates the simple “compression” mounting of the J -seal, in which the seal — 
leaf is directly attached to the skin plate; it also takes advantage of the 
greater leaf flexibility of the J-seal’s “deflection” mounting, while circum- 
venting the “outboard” mounting requirements of the latter. Perhaps the — 
most outstanding advantages of the J-seal are its versatility and low cost; 4 
several manufacturers of industrial rubber products offer a wide standard — 
range in J-seal sizes and molded corners. The writer would be very much — 


x 
. 
a 


interested in learning how the L-seal negotiates the bottom corners of a : 
radial gate. § 


His small contribution to an engineering handbook(9) has made the writer 
aware that a comprehensive coverage of the problems besetting the designer 
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f tainter gates would be a voluminous undertaking. He also wonders whether 
1¢ limited audience would even pay the printer’s costs. 

The writer gratefully acknowledges the cooperation of Mr. Richard D. 
larza, A.M. ASCE, of the Harza Engineering Company, Chicago, IIl., in his 
onsent to make public certain observations derived from the unpublished 
rorks of the late Leroy F. Harza. 
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CALIBRATION OF A SUBMERGED BROAD-CRESTED WEIR 2 


Discussions by P. K. Kandaswamy, N. Rajaratnam, 
and T. R. Anand 


P. K. KANDASWAMY, | M. ASCE and N. RAJARATNAM2.—The author has 
presented a very interesting paper on broad-crested weirs with upstream and 
downstream slopes. The new method of approach presented therein to find 
out the modularity limit based on the Froude Number at the center of the 
‘weir crest is commendable. 

The flow over a weir of any shape can be classified as modular or non- 
‘modular as explained below. When the discharge over the weir depends only 
‘on the upstream head and is independent of the tail-water head, its perfor- 
‘mance is said to be modular. When with an increase in tail water depth, this 
independence ceases and the discharge depends upon both the upstream and 
‘downstream head, its performance is said to be non-modular. The limit be- 
‘tween these two ranges is called the “modular limit” or “modularity limit.” 

In the modular range, the discharge over a broad-crested weir with verti- 
cal upstream and downstream faces is given by the general equation. 


q@- cB H?/2 (1) 


Q is the discharge in cusecs, 

C is the coefficient of discharge, 

' B is the length of the weir in feet normal to the direction of flow and 
is the head in feet over the crest of the weir. 


Investigations(1,2,3,4) by various research workers have shown that the coef- 
ficient C varies mainly with H/L where L is the width of the crest in the 
‘direction of flow and is independent of the velocity of approach represented 


‘by the term a where P is the height of the weir crest above the approach 


channel floor. The curve given in Figure 1 shows the variation of C with 

‘#H/ L for a broad crested weir with vertical upstream and downstream faces. 

: When the broad crested weir is provided with upstream and downstream 
slopes, the flow pattern and hence the discharge characteristics are likely to 
be different. Considering the downstream slope first, it is known that when 
water flows over a broad crested weir, the critical depth occurs over the 
erest. The exact location of the critical depth section is difficult to define as 
can be seen from Figure 2. To be sure, the stage of the flow changes from 


a. Proc. Paper 1960, March, 1959, by Richard A. Smith. 


EE. Executive Engr., (Research), Irrig. Research Station, Poondi, India. 
: Junior Engr., (Research), Irrig. Research Station, Poondi, India. 
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Legend 
Pfra, L = 0-50 Foot 
Pfra, L = 1-00 Feet 
WoodburnL = 2-00 Feet 
Basin L = 6-56 Feet 


WoodburnL = 10-00 Feet 
Poondi L = 1-00 Feet 
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Values of H/L 


Fig-| STANDARD COEFFICIENT CURVE FOR BROAD CRESTED WEIRS. 
(From Ref- 4) 


the sub-critical to the supercritical condition over the crest. Hence, so long 
as the downstream slope provided is steeper than the critical Slope for the 
range of discharges involved, as it happens to be in the case considered by 
the author, any change in the super-critical downstream slope is not likely 
to affect the flow pattern over the crest. Consequently, the discharge coef- 
‘icient C in the equation Q = CBH3/2 will be independent of the down- 
stream slope. If in any case, the downstream slope is milder than critical 
slope for the discharges involved, the tail-water may back up and the 
modularity limit may occur at some particular slope which may be termed 
is the “modularity slope.” 

It has been observed in a general way, that the discharge coefficient for a 
sroad-crested weir increases as the steepness of the upstream slope de- 


sreases. To incorporate the effect of upstream slope in the general formula, 
squation (1) is modified as 


Q=Cy BH3/2 (2) 


vhere y is known as the “slope factor.” The variation of y with H/L is 
shown plotted in Figure 3 for values of upstream slope varying from 1/100 
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Fig:2- WATER-SURFACE PROFILES AND POSITIONS OF 
CRITICAL DEPTH ON WEIR WITH HORIZONTAL CREST 


From Ref: 3) 


to 1/1000. The results are taken from a series of experiments conducted in 
the Irrigation Research Station, Poondi on bye-washes which are flush escapes 
for floods provided in small tanks in the State of Madras. Further experi- 
ments have to be conducted to collect more information for steeper slopes up 

to 1/100. 
| The author has measured the depth of flow at a point located slightly 
downstream from the center of the weir crest and has calculated the cor- 
responding Froude number. This has been used for defining the modularity 
limit, that is the stage at which flow changes from modular to non-modular 
range. When conducting similar experiments at this Research Station, the 
writers have observed that the water surface over the crest is undulating and 
that sometimes the same depth occurs at more than two sections. Hence, 
the writers feel that it may not be very reliable to base any solid conclusion 
on the depth measured at the section mentioned previously. 

Considering the functioning of a broad crested weir in the modular range, 
any increase above the crest up to the modularity limit does not affect the 
upstream level. This apparent anamoly has been explained by the writers(5) 
with the conception of the critical depth section occurring as a control section 
over the crest and of the celerity of the breaking surge. In the paper referred 
above, it has been concluded that until the modularity limit is attained, the 
| ‘critical depth occurs somewhere on the crest. In fact, it can be seen that 
with a steady increase in the tailwater level, the critical depth keeps on 


“moving downstream over the crest up to the modularity limit when it falls 
into the relatively deeper downstream portion. This conclusion agrees in- 
‘directly in a qualitative manner with that of the author that the change in con- 
_trol occurs when the Froude number is approximately 0.78 at the section 


chosen. 


<a nd 


Heer 


EA 
ne 


_ 


a 


is 
Sie eel aazeesnesten Se 
se * 3 eer ttt 
+ PH 

Fo 

oon 

aan 

a 


Poo 
Beer 


rrr ni 


+ aso ag 
t ee are 


22S Se ee 


sa 
A) 

1 we 7 
ai a fe 


a 
G 


2 BR z Ha) 
- |. 


ap. 
nm 


Further, the second for a broad crested weir in the non-modular range 
is given by the relation(2) 


where Q is the actual non-modular discharge, 
Qm is the modular discharge for the given upstream head and f isa 


reduction factor which is a function of the submergence ratio a and Z, Z 
being the height of the tail water above the crest. 

The variation of the reduction factor f£ with the submergence ratio 2 fora 
weir called Kodiveri Anicut(5) is given in Figure 4. Studies on Dowlaishwaram 
Anicut(6) indicate the variation of f with H/L. For a broad crested weir with 
upstream and downstream slopes, it is likely that : 


t-9{2,3,s,, sgh (4) 
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where Su is the upstream slope and Sq is the downstream slope of the weir. 
The nature of the function can be studied by detailed experimentation. 
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T. R. ANAND, ! A.M. ASCE.—The writer heartily agrees that previous 
literature on the broad-crested weir coefficients was incomplete. The 
analytical relationships worked out by the author between Froude number, 
degree of submergence and weir coefficients (illustrated graphically in 
Fig. 3 and 4) give a sound basis for selection of proper coefficients in parti- 
cular cases. 

The writer has been concerned with selection of weir coefficients for es- 
timating discharge obtainable if reservoir dykes (rock fill or earth fill with 
adequate top and slope protection) are designed as broad-crested weirs to 
pass flood flows. From Fig. 4 it appears that for weirs with side slopes 1 
on 1 or flatter, a coefficient "C’ = 3.0 may be used for submergence up to 
70% and where there is no likelihood of tail water rising above weir crest, 
this value may be raised to 3.2. Perhaps it would be necessary to reduce 
these values somewhat for design work to allow for variation of roughness 
coefficients of the weir materials, velocity of approach etc. etc. Laboratory 
work done by L. W. Long and N. W. Conner(1) indicates that at high flow 
rates the discharge coefficient of broad-crested weirs is lower for paraffin 
surface than for glass. The writer would be interested to know to what extent 
Prof. Smith’s values were modified for application in the design office to 
determine discharge over road-beds. 

Figures 2 and 4 indicate that the effect of increasing side slopes of weirs 
from 1 on 1 to 1 on 2 on the discharge capacity is very small (e.g. if weir 
crest length 1000 feet and total head 10 feet, discharge capacity =134,500c.f.s 
for 1 on 2 slopes and 127,000 c.f.s. for 1 on 1; % difference of discharge = 5% 
only). This leads to the useful derivation that flatter side slopes do not in- 
crease discharge rating significantly and considerable savings may be made 
by building weirs to say 1-1/2 on 1 slope instead of 1 on 2, or flatter, pro- 
vided other conditions such as structural stability, hydraulic gradient, etc. 
etc. are satisfied. 
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TABULAR SOLUTION OF OPEN CHANNEL FLOW EQUATIONS 2 


Discussions by R. Silvester, Ivan M. Nelidov, 
Walter J. Tudor, and Warren W. DeLapp 


R. SILVESTER!.—The author has provided a very useful table giving 
values of specific energy readily and accurately. The author’s equation 3 has 
appeared previously.(1) The analysis of the dimensionless form has been 
restricted to the rectangular channel. It can be modified to include other 
shapes as follows: 


Spy Sant Rea oat _ de(bule 
Es Sa Ae pee healed c w/c 
d. d. 2ga*d. ay Z asin: devs 


where by, is the water surface width and b‘, is the width of a rectangular 
channel equivalent in area to the original channel; the suffix c denotes values 
at the critical depth. 
d 


thus aS =— ee de (bile 
do do 2 a’ d. (by). 


E, d (se) (bi). 


or — son +— 
de de 2\d/ Lbs (bi). 


As seen in the accompanying sketch values of by and bly can easily be 
determined for trapezoidal channels for any stage operating. For other 
shapes a scale drawing could provide values close enough for most practical 
| ‘problems. Where depths are being determined an initial estimate based on 
the original equation could be the basis of a correction applied for the shape 
factor. 

This procedure may appear a little arduous in order to obtain one or two 
| Specific values as illustrated in the problems cited. But the table could be 

| used in the step method solution of backwater profiles. It can be shown that 


E-F dx _ E2/de —E./de 


a z Se-S¢ de Se — S¢ 


a. Proc. Paper 1961, March, 1959, by H. A. Babaoek. 
1. Senior Lecturer in Civ. Eng. (Hydraulics), University of Western 


Australia, Nedlands, Australia. 
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where dx is the length of channel over which Ss is the friction gradient and 
Ej and Eg are the specific energies at the upstream and downstream ends 
respectively; S, is the bed slope. 

By knowing the stage at each end of a section of channel, length dx, the 
specific energy can be determined from the table supplied by the author. For 
a rectangular channel this is done directly. For other shapes the correction 
outlined above can be used. 

The correction in fact is only applied to the velocity head term and, as can 
be seen in the author’s Figure 1, this is small for mild slopes once d/d. > 1s 
Below this value down through the critical depth and for steep slopes the cor- 
rection assumes importance. The velocity head is found by subtracting d/d, 
from the value of E,/de as read in the table. An example should assist to 
illustrate the points raised. 


Example 


A channel of bottom width 25 ft, side slopes 1 on 1-1/2 and bed slope 
0.00079 delivers 900 cusecs into a reservoir whose depth at the junction is 
10 ft. It is required to find the distance upstream at which the depth is within 
2 percent of normal; Manning’s n = 0.025. 


Note: Graphs and a nomograph for ease in determining normal and critical 
depths and other terms in the calculation of S¢ are contained in articles by 
the writer(2,3) to which the reader could refer for checking the computations. 


K Q 900 
o= ~ /S, Me" = = = 

ry vz fa2.2 
900 = <o35 [0.00079 = 159 


gal apie de = 3.28 ft. 
-thus the normal depth - since do>de and d>do 
do = 6.0 ft. curve is mt 63)(4). 


Thus the total length of the backwater profile is 3445 x 3.28 = 11,330 ft. 
It may be remarked here that the same problem was solved by the conven- 
tional step method in reference 3 where the answer given was 14,360 ft; but 
if it were reworked using the average values of friction and kinetic energy 
terms for each channel section, as above, the answer would become 11,300 ft: 
Actually using the table without the shape factor correction results in a lengtl 


of 11,100 ft. 
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(a) value as taken from table. 
(b) value corrected for shape factor. 


(c) equal to 0.040/0.05 «10? 


IVAN M. NELIDOV, ! F. ASCE.—The writer was attracted by the method 

|presented by the author for facilitating solution of open channel problems. 
In accordance with the preference of the writer for graphical solution of 

‘hydraulic problems he wishes to present a graphical interpretation of this 


problem for trapezoidal channels. ‘ A 
Since for critical flow the velocity head is hyo = a = se and 
c 2g 

V 3/2 


ae Ag? /2 ¢ 
a= rug £5 then Q, = AcVe = 475 OY and Qe/g1/2 = 7173 = Fy- 
c 


: For trapezoidal channel: 
(B+ mdc) 3/2 3/2 (B/m + de)3/2_ 43/2 
ves ec. a Sy Ene ae 
1° (B+ 2md,) (B/m + 2d.) 
The curves representing F, in terms of dc and B/m are plotted on Fig. 2. | 
For the specific energy related to critical depth the following expression 


can be obtained: 


ove sme od ay 
d. d. 2gd, c 2gd, AZ 


" 


d+ de 9 Ac,3, d \2 
e* 4/2--)*( (4) = 
do ada! Sk ae 
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de 2 (B + mdc)3 =F 
a) (B+ amd,)(Be maj? *2» Where: 
2 (B/mde + 1) 


SURO’ 
Fo-—g> t V2) Bfmdg 2) (B/mdy = dae)? 


The curves representing Fy = Eg/d, are plotted on Fig. 3 in terms of d/d¢ 


and B/md. The application of these curves will be illustrated on the examples 


of the author. 
Example 1. Data as per author, except that the side slopes are 1.5 to 1, or 


m-=1.5. Then, Q,/ 1/2 = 400./5.67 = 70.7 and B/m = 9./1.5 = 6. From 
curves on Fig. 2 we obtain d, = 3.25 ft. and E,/d, = 7./3.25 = 2.15. With this 
value, we obtain d,/d, = .57 and dt/d, = 2.11 from curves on Fig. 3. This 


gives us d, = .57 X 3.25 = 1.85 ft. and d, = 2.11 X 3.25 = 6.85 ft. 
Example 2. Data as per author, except side slope factor m=2. Then 
Q/,1/2 = 300./5.67 = 53. and B/m = 10./2 = 5. From curves on Fig. 2 we 


obtain d.) = 2.50 ft. and dy/d,) = 6-/2.5 = 2.4. With B/md,) = 10./2x 2.5 =-2. 


we obtain from curves on Fig. 3 EsV/ do} = 2.43 and Eg) = 2.43x2.5 = 6.07 ft. 


For the narrow section, with the same side slopes, we obtain: with @/g1/ 2 


= 53. from Fig. 2 and B/m = 9/2 =4.5, dog =2.65 ft.; neglecting losses, 
Es2/dc9 = 6.07/2.65 = 2.29. Then, from Fig. 3 and with B/md,9 = 9./2x 2.65 

= 1.7, we obtain dp/d.o = 2-32 and dp = 2.32x2.65 = 6.15 ft. 

| Example 3. With data by the author and side slope factor m = .5, also with 

: Q,/gl/2 = 250./5.67 = 44.1 and B/m = 8./.5 = 16, we obtain d, = 2.95 ft. Then 
dy/d, = 5./2.95 = 1.7. From Fig. 3 with B/md, = 8./5X2.95 = 5.43 we find 


By/d, = 1.83: Eg9/dy = 1.83 - -5/2.95 = 1.83 -.17 = 1.66; with this we obtain 


| do/d, = 1.48 and do = 1.48x 2.95 = 4.39 ft. 
| For circular conduits if we consider the central angle a subtending the 
water width T, the same will have an expression T = DSin 1/2a, where D is 
‘the diameter of the conduit; the depth of water will be d = D/2(1-Cos 1/2 a) 
and the area of cross section A = D2/8 (a - Sin a). The U. S. Bureau of 
Reclamation tables contain coefficients necessary for plotting curves similar 
i to those for trapezoidal channels. (Table 12a, p. 48a, Ed. 1923). The rela- 
tionship between the critical depth ratio to diameter d,/D is given by coef- 
ficients Q,/D°/2, which we will call fy(d¢/D), so that Q¢ = £1 (d./D)D®/2; 
the relationship between d/D and the area is given by the coefficients A/D# 
(p. 48), which we will call f9(d/D), so that A = f9(d/D)D2. Then Fp = Es/d, 


2 

Q2c f1“(d,/D) 

= d/d,, + 1/2 ——> = + 1/2 which can be plotted for 
c / gd. AZ d,/D / fo” (d/D) 3 


various ratios of d/D and d,/D. " 
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WALTER J. TUDOR, ! A.M. ASCE.—The cubic equation which arises in 
ydraulic problems can be obtained from the author’s equation (3) 


2 
Bren sieh{tts) 
d. d 22\ ad 


Ny multiplying by d.d2 and substituting 


3 
Be ees ge 
2 2g ; 
to yield qo Bde re oe FO 


This cubic always has two positive roots or depths (that is, tranquil and 
shooting) when Es and the constant term are positive. A nomographic solu- 
ion for the extraction of these two positive roots is shown on Figure 4. 
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As an example illustration of the diagram consider the author’s first 
example. The cubic equation to be solved is 
a> - 7a* + 30.6 = 0. 


2 
Since the scale of ae extends only to 10, set d = 2d, to obtain 


eres: 28d,*+30.6 = 0 


d 


3 2 a 
Pe ids 3 Aye i Ss BS ae 0 


2 
A straight edge through = = 3.83 and the pivot point intersects the 


curves E, = 3.5 at dy = Jeoo ane dy eae Lp 
ae d. = 2.66 and d, ties 


The author’s second example is solved as follows: 
E: 30, 4 
Es = 6 + (“7) /2g = 6.388, 
so that the pertinent cubic equation is 
304 
2g 


a> = 6.388 d* + 


or a a By 39 A oF zee Ay 


Let d=1.5d,, so that 


S988 dys elags ay Te SD 


a d,> - 4.23 dj* +5,09 = 0. 


2 
Enter the chart with E, = -4.23 and a = 5.09 to obtain 


He 30 and d = Se 92. 


[en 
_ 
" 


J. dg =i 2595 * and dy = 5.88; 


This nomographic chart can be constructed to various sizes for specific 


2 
cases and desired accuracies, and the oa scale can be calibrated directly in 


terms of q. However, the main purpose of the discusser is to demonstrate 
the nomographic solution of the hydraulic cubic. 
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WARREN W. DeLAPP,! F. ASCE.—The author has demonstrated that the 
dimensionless specific-energy function, in either graphical or tabular form, 
can be a useful tool in the solution of certain open channel flow problems. 
The fact that the method presented permits the inclusion of energy losses, 
either as a numerical value or in terms of the velocity head, gives it a dis- 
tinct advantage over previously published procedures. This was pointed out 
to the writer by the author some time ago when the latter was a graduate 
student in a course in open channel flow given by the writer, and since that 
time the method has been included in the course material. 

The use of the diagram and table of the paper are restricted to the rec- 
tangular cross section, and the method cannot be used for the trapezoidal 
section because the ratios involved change with the geometry of the channel. 
There are, however, two other shapes, less commonly used, for which a simi- 
lar procedure may be used. These are the triangular and parabolic channel 
cross sections. Since the determination of the depth for a given energy for 
these sections is even more difficult than for the rectangular channel, con- 
siderable work can be saved through the use of the author’s methods if either 
of these shapes are being used. 

A general criterion for critical flow is found in most textbooks and is de- 
rived by the substitution of Q/A=V in the author’s Equation 1 and setting the 
Herivative of the expression equal to zero for minimum energy. Noting that 
dA/dd=T, the width at the free surface, results in the expression 


ah nae (9) 
Ga-5 57. 
which is satisfied in any open channel when the depth is critical. 
Consider now a symmetrical triangular channel, vertex down, with a ver- 
tex angle 8. Both A and T in Equation 9 may be expressed in terms of d and 
tan 8/2 and the critical depth is given by 


: 9492Q° (10) 
: ae V3 tan* G2 


The equation for specific-energy is 


fg zd t alae = d+ rant (11) 


This is reduced to dimensionless form by following the author’s procedure 
and dividing Equation 11 by Equation 10 which may be simplified to yield 
| 


Ee pad Cp (4)" (12) 
de Zed a 
ch can be shown by the same procedure to be valid for an unsymmetrical 


‘iangular section as well. 1/2 
_ The parabolic channel may be defined by the equation T=kd4/4 and 


~ Prof. of Civ. Eng., Univ. of Colo., Boulder, Colo. 
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proceeding as above from Equation 9 the critical depth is 


Ppa Aa (13) 


8Gk* 
and the dimensionless specific-energy equation 
cee et / eet (14) 
EP LE META Et 


It is evident that, as for the rectangular section, either Equation 12 or 
Equation 14 may be used in graphical or tabular form for the solution of 
problems involving their respective shapes and that the curves which they 
represent are independent of the values of 8 and k. 

It may be of interest to note that the pressure-momentum function, fre- 
quently used in connection with the hydraulic jump, may also be expressed in 
dimensionless form for each of the above channel shapes. For the rectangu- 
lar channel of width b the pressure-momentum is given by 


mee a4 : 7 Q* (15) 
P+M= z of gba 
in which -y is the specific weight of the fluid. The critical depth is 
{es 
d.=z 16 
c 6°9 ( ) 


and making the substitution Q2 = a3 b2 g and dividing by do gives the 
dimensionless form of the equation 


<PrM ted d 17 
yea? = zlg)t > “a 
In the same way the expression for the triangular channel is 
PEM SS ctsaee etl Toe 
Fin a@® = og) +z) se 


and for the parabolic section 
mM i: | 
74 a -B (E) + $ @ i a 


Equation 17 for the rectangular channel is shown graphically in Figure A, 
curves for the other channel being similar inform. The use of the curve is, 
of course, quite simple. Knowing either of the depths, say dy, for a hydraulic 


jump, the ratio d)/d, is calculated. The ratio dy/d, is then selected from — 
the point on the curve having the same abscissa value as do/d, and the value 
of dy then calculated. Although the conjugate depths for a rectangular | 
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Figure A. Dimensionless pressure-momentum diagram. 


channel may be found readily by other means, this is not true for the other 
shapes. The conjugate depths for the latter can, however, be easily deter- 
mined from a graph such as Figure A drawn to an appropriate scale, or from 
a table similar to the author’s for a higher degree of accuracy. 
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RELATIONSHIPS BETWEEN PIPE RESISTANCE FORMULAS2 


Discussions by K. Szesztay, Henry J. Miles, 
M. B. McPherson, Nicholas Bilonok, and J. V. Dawsey 


K. SZESZTAY!.—Recent investigations into the movement of fluids provide 
. common theoretical basis for the analysis of flow in closed conduits and in 
pen beds. The older exponential formulae used for the design of closed con- 
uits are compared in a very clear manner to the inclusive formula and dia- 
ram of Colebrook-White, respectively Rouse-Moody in the paper of 
fr. W. Moore. 

These formulae rest on a sounder theoretical basis and are of a more 
eneral validity so that their use in computations relating to open beds would 
desirable. The main obstacle encountered here is the difference in units 
f roughness. 

The roughness of open beds is usually expressed by the n value of Kutter- 
ianguillet, while the equivalent particle diameters k of Nikuradse are used 
n theoretical investigations and in design computations of closed conduits. 

Combining the Karman-Prandtl formula for rough pipes and the Manning 
oOrmula H. Rouse developed for the relationship between the two roughness 
arameters the formula:(1) 


pi/é 


0,093 = 2 log — + 2,35 [ft] (1) 


n 
In order to study the effects of factors which cannot be taken into account 
nalytically (kinematic viscosity y and the gradient S) and errors in observa- 
on, a statistical investigation has been carried out recently. The evaluation 
f the 116 data obtained experimentally, respectively by measurement, and 
xtending over the entire turbulent range in both open beds and closed con- 
wits, yielded the following direct relationship of two variables (Fig. 1):(3) 


| (2-a) 
: — = 19,8 log 2 [cm] , 

espectively 

|i + = 19,8 log —22 = [rt], (2-b) 


Owing to the scatter due to the simplifying assumptions and uncertainties 
entioned above, the statistical investigation failed to reveal the effective of 
lues of R. 

In his paper Mr. W. Moore points to a very important prerequisite for the 
re widespread application of the lately developed formulae, which although 


Proc. Paper 1962, March, 1959, by W. L. Moore. 
Cand. of Eng. Sc. Researche Institute for Water Resources, Budapest, 
Hungary. 
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Fig. 1 


of a more general validity, are of a more complicated form. This is the 
problem of readily usable graphs, whereby the need of successive approxi- 
mations can be eliminated. An easy solution can generally be found to this 
problem in case of conduits having a circular cross section, all characteris- 
tics thereof being the function of a single parameter, the diameter. In the 
case of the most common type of open canals, the trapezoidal cross section, 
the situation is considerably more involved. 

Graphs were recently issued by the Research Institute for the hydraulic 
design of canals and pipelines. (4) Attempts have been made therein to give 
a direct solution for the dimensioning of the optimum trapezoidal canal cross 
section for any given Q discharge, § terrain slope and n, respectively k, 
roughness conditions. 

The problem will not be unambigous, unless certain shape characteristics 
of the trapezoidal cross section are either assumed as constant, or con- 
sidered as a direct function of the water depth (or any other dimension). 
Similar basic assumptions for the construction are shown in Fig. 2. 

Variations in the shape characteristic B = a/h with the waterdepth h, hav 
been adopted after studies on stable beds (Lacey, Lindley, Girskan and others 
and on basis of empirical data. The relationship between the bank slope 9 _ 
and the waterdepth h represents a statistical summary of solutions common. 
ly used for unlined beds. It should be noted, that the method to be presented 
below can be extended by the construction of an adequate number of graphs to 
constant, respectively to any arbitrary, 8 and Q values. digs 
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Fig. 2 


Having selected the values of 8 and 9, respectively the corresponding 
ondition equations, the construction of the graphs constitutes a unambigous 
ydraulic problem. The first step is to establish the relationships 


a= Bh (3-@) 
b-hh+2gh (3-b) 
D-Bhe2h Weg? Ge) 
A-n? (8 +9) (3-0) 


ne(B--9) F 
A Bh Oh 1+9 (3-€) 


is the wetted perimeter, A the cross sectional area), expressing the 
iriation of the hydraulic parameters of the cross section with waterdepth. 
In the relationship expressing the conveying capacity of the beds 


| Qs AveAc JES (4) 
e velocity factor c is computed from the expression 

: k 0,222) 

se — emw- 132 lo 0,0676 —— sO ) (5) 
Boas gS: (050078 TAA ee 


= 981 cm/sec2, the constant of gravity) developed as a transformation of 
Colebrook-White formula by E. S. Crump.(2) 
The inclusion of the second term in Eq. (5), expressing the effect of vis- 
sity and slope, is not necessary in case of conduits and canals carrying 
r unless k < 1 cm, respectively n <o,017.(3) 
Consequently, with the substitutions according to Eqs. (3-d), (3-e) and (5) 
‘discharge expressed by Eq. (4) is - for any given roughness coefficient 
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k - actually a unambiguous function of the surface slope S and the waterdept 
h. 

~ A direct graphical solution for the roughness coefficient n = 0, 025 (respec 
tively k = 7 cm) is shown in Fig. 3. When using the graph, the mean flow 
velocity v and the waterdepth h is determined for the given discharge Q 
and surface slope Ss, whereafter the B and 9 values corresponding to the 
waterdepth h can be obtained from Fig. 2. With these values available, all 
further parameters of the cross section can be computed by Eqs. (3), or can 
be obtained from the corresponding graphs. 
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The values obtained from the graphs may be used as a first approximation 
even if local conditions as regards roughness, bank slope, or the ratio a/h 
differ from those adopted originally. 

By compiling a series consisting of an adequate number of sheets, the 
direct solution according to Figs. 2 and 3 can be extended for any combination 
of n, 9 and a/h occurring in practice. 
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HENRY J. MILES, 1 M. ASCE.—Professor Moore’s paper is a valuable 
contribution to a better understanding of the subject of pipe friction. The 
writer is particularly interested in this paper as he had independently de- 
rived a similar relationship between exponential formulas and the general 

‘resistance diagram, and for many years in his classes in hydraulics has been 

‘using this idea for introducing the overall concept of fluid friction. As the 
author points out, the large number of exponential formulas is not only con- 

fusing to the beginner but also gives the wrong impression of the inter- 

‘relationship of hs, D, Q, and V unless it is realized that each formula applies 

only to a limited range of Reynolds number. 

The writer has used a slight modification of the author’s formula (6), 

hamely 

-h nm 

is tops AGB nav’ () 

: = 2g D oe. 

In the range of laminar flow the plot of f vs. Np is a straight line on the 

: general resistance diagram with a negative slope of 1. Therefore, m = 1 and 


the formula becomes 
| m2 
| 2g D 


220 7LV (2) 


‘For high Reynolds numbers the slopes of the lines approach O, m approaches 
2, and the formula becomes 
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2 
OL ev 
he = 8 tees, (3) 
2g D 

the viscosity term having disappeared completely. For intermediate Reynolds 
numbers, as the author points out, the value of m varies from about 1.7 to 2. 

Thus it is possible to visualize the general picture of friction loss as 
Reynolds number increases: hg varying initially as V and inversely as D2, 
and gradually changing to the relationship of varying directly as V2 and in- 
versely as D; viscosity initially playing an important role but gradually de- 
creasing in importance until it finally disappears from the equation. It is also 
apparent from observing the slopes of the lines on the general resistance 
diagram that the rougher pipes approach complete turbulence at much lower 
Reynolds numbers than the smoother ones. 

It will be noted in the writer’s formula (1) that the sum of the exponents of 
D and V always totals 3. This tendency can be observed in the various ex- 
ponential formulas proposed in the past by Hazen-Williams, Scobey, Barnes, 
Lea, and others. Since the sum of the exponents of D and V in these ex- 
ponential formulas does not always total exactly 3, the question arises as to 
whether the discrepancy is caused by slight errors in the empirical formulas, 
or whether the idea of representing corresponding values of f and Np on the 


general resistance diagram by smooth curves is only approximate and intro- 
duces small errors. 

In conclusion, it is interesting to use the methods proposed by the author tc 
compute the approximate numerical relationship between the e/D ratios of 
Nikuradse and Hazen-Williams C. The latter formula can be written in the 


form 


8 
308 ti Se 


nian erent (4) 
£ oe pi-17 


If the exponents of V and D are adjusted slightly so that they total 3, the 
formula can be equated to the writer’s formula (1) 
3.02 L yi ~B4 oe yern rr ya 
he © = (5) 
ogee” pi-16 2g p72 


from which it is obvious that m must equal 1.84. If this value of m and the 
kinetic viscosity of water at 60° F. are inserted in this equation it reduces to 


©; ies (6) 


which, if substituted in the author’s formula (5) 


teal 


2-m 
N, 


reduces to 


go 85. 28 1e67, : (7) 


ry 9 DISCUSSION 143 


If corresponding values of f and Np are read from the general resistance 


iiagram for each e/D line at the point where its negative slope is 2-m or 
).16, and substituted in this last equation, Table IV is obtained. If these 
alues are plotted on semi-logarithmic coordinates, the resulting curve ap- 
roaches asymptotically Cy = 160. 


Corresponding Values of Nikuradse’s e/D and 
Hazen-Williams’ C, for Water at 60°F. 


1 
e/D ro 

.05 93 

-04 97 

03 102 

02 108 

.015 113 

01 118 

.008 121 

.006 124 

.004 128 

.002 135 

001 140 

.0008 142 

.0006 144 

.0004 146 

.0002 150 
: .0001 153 
: .00005 155 
00001 158 
: .000001 159 


Mz B. McPHERSON, ! M. ASCE.—Since the emphasis in the author’s pre- 
ntation is directed towards a precise interpolation of the ‘general resistance 
iagram,’ it is necessary to review limitations and assumptions which must 


considered in using the diagram. 
General Resistance Diagram 


Colebrook, (5) noting that the general trend for new commercial pipe data 
as a continuous transition between the characteristics of hydraulically 
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“smooth” and artificially “rough” pipes, arrived at an asymptotic “transitior 
law” of good approximation based upon integration limits for these extremes 
given by Nikuradse. The Colebrook “transition law” is graphed in Fig. D-1, 
a photo-copy of Fig. 1, p. 136 of reference (5). 
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Figure D-1 q 

Reproduction of Figure 1, p. 136, Reference (5). | 

(Note thata =f, k=e, d=D, Vy=Vvi/8 and R=Np). — 

(The expression for the “transition law” is equivalent to Equation 167 by — 
Rouse(3)), -F 

7 
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“Roughness V” represents a non-uniform artificial sand grain roughness 
in which 5% (much larger size than the remainder) was cemented in a diamond 
pattern in a test pipe by Colebrook. The original raw data for the three 
classes of new commercial pipe were selected from a variety of published 
American and private British references. Colebrook determined the equiva- 
lent sand grain roughness, e, for each individual commercial test pipe direct- 
ly and indirectly by means of the Karman-Prandtl “smooth” and “rough” pipe 
relations. From these calculations for each pipe, he then determined the 
“mean” e for each class. The extremes and mean for each class are sum- 
marized in Table D-1. Using the calculated mean e for each class together 
with the raw test flow data, Colebrook then determined empirically the “mean 
transition curves” of Fig. D-1 for each of the three classes of commercial 
roughness. Rouse(3) has compared a variety of test data (mostly small 
diameters) with the Colebrook “transition law” in his Fig. 110. It is to be 
noted that Colebrook interpreted his results in terms of a system of transi- 
tion curves, whereas only his semi-empirical “transition law” has been in- 
corporated in the general resistance diagrams presented by Rouse(6) and 
Moody(7) as being “a close approximation to the actual resistance law.” 6) 
Inasmuch as the majority of the tests supporting the validity of the transition 
region curves of the general resistance diagram were for new pipes of small 
diameter, there remains some question as to the probable accuracy in the 
use of this diagram for large pipes. 


(Dise. 1962) 
TABLE D-1 


Review of Colebrook Summaries -- Ref. (5), pp- 152-153 
New Pipe--Flow of Water 


"A sphalted (tar coated) Cast Iron Pipes" 
Diameter, inches Approximate e, inches 
4 -0090 
6 0042 
8 (3 tests) -0024, 20060, .0067 
| 12 (3 tests) .0027, .0030, .010 
| 16 (2 tests) -0025, 0084 
: 24 0055 
| 40 0074 
| Wes 0062 
| 48 20047 
él 20016 
For the Three Types Considered by Colebrook: 
Type "Mean" e Approx. Max. 6 Approx. Min, ¢ 
Galv. Iron* 0.006-in. (.0005-ft.) 0.0ll=in. 0.001-in. 
| Asph. C. I. 0.005-in. (.0004-ft.) 0.010-in. 0.0016-in. 
i Wrought Iron 0.0017=in. (.00014-ft.) 0.0038-in. 0.0004-in. 


(* 0.35 to 4-in. pipes; * asphalted) 


Note: "Approximate" means that the figures cited were read by the writer from 
' a graph of plotted points. 


Ce 
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The mean G. I. transition curve of Fig. D-1 would lead to values of f up 
to 5% higher than the “transition law.” If the range of calculated e values 
for G. I. pipe given by Colebrook (Table D-1) are applied to the general re- 
sistance diagram, Figure 1, the variation in f about that for the mean e is 
of the order of plus and minus 15% for a 4-inch pipe at a Reynolds number of 
105. This is a reminder that the first four equivalent uniform sand grain 
roughnesses for new commercial pipe tabulated in Fig. 1 are mean values. 
(Ranges of values, such as those given by Streeter, (18) would seem to be less 
misleading). Most references suggest that the accuracy of estimating f in 
the “rough” region (for which there is a paucity of supporting data) is 
probably within about plus or minus 10%: for “smooth” pipes the probable 
variation is about plus or minus 5%. Recall also that the supporting data for 
the transition region analysis was selected from laboratory and field data 
considered representative of new, clean straight pipe without fittings 
(Colebrook made arbitrary head loss deductions for bends and fittings from 
some of the field measurements). 

The universality of the “smooth” pipe characteristic for steady near- 
incompressible conditions has been demonstrated using oil, air and water 
(reference (3), p. 203). Nikuradse’s experiments as well as all the com- 
mercial pipe tests used by Colebrook and subsequent adaptations were con- 
fined to water flow. The trends of the functions for new commercial pipe with 
a random or ‘statistical’ type of roughness are indeed remarkably similar. 
However, the preciseness of the general resistance diagram implied by the 
“exponential equations” of the author does not appear to be justified. The 


author leaves the impression that the characteristics of new commercial pipe 


are fully established and completely understood. The extensive studies of 
artificial random sand grain roughness by Colebrook and White(28) and 
Square-groove machined-roughness patterns by Streeter(19) indicate that 
much more remains to be evaluated. Rouse(20) warned some time ago 


against leaning too heavily upon a single linear dimension e, in attempting to — 


describe any roughness other than uniform. In reference (6) Rouse recom- 
mended the Colebrook relation “as a close approximation to the actual re- 
sistance law.” From the standpoint of velocity distribution, Robertson(21) 
states that the effect of commercial pipe roughness “has not yet been ex- 
plained in terms of our knowledge of the action of artificial roughness.” The 
line of demarcation for m = 2 is arbitrary and by no means exact. These 
remarks are not intended to disparage the vast and invaluable knowledge ac- 
quired to date, but to warn against the false sense of exactness implied by 

a too-literal interpretation of the trend-curves as in the author’s « exponen- 
tial equations.” 


Empiricism in Design Practice 


The writer assumes that the instruction in fluid mechanics tendered 
American undergraduate civil engineers for more than a decade has left at 
least a residual impression, and is confident that contemporary designers 
(mostly young men) are conversant with sound mechanics principles and the 
limitations inherent in commonly used empirical relations. He assumes also 
that a comparison between the Hazen-Williams formula and the general re- 
sistance diagram, similar to the presentation by Vennard, (22) is offered to 
the majority of civil engineering students. The timely, comprehensive paper 
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yy Hinds(12) in 1946 coincided with the advent of widespread intensive fluid 
mechanics instruction in this country. There has certainly been a visible up- 
rrading in American design practice since that time and it is hoped that the 
withor’s pessimism does not truly reflect prevailing conditions except in the 
sxxtreme. The relatively contemporaneous heated discussions on formulae to 
which the author alludes(10,14) are specific to British practice; some of the 
ormulae in these references are peculiar to British usage. 


New Express or Feeder Mains 


Responsible designers usually appraise the expected initial, or “new,” 
sapacity of this type of main by means of the general resistance diagram. 
since designs are normally based upon ‘future’ requirements, an estimate of 
oth new and aged roughnesses is necessary. Recourse is made to well- 
locumented and carefully selected data such as presented by the U.S.B.R. (23) 
ind/or to local records in estimating pertinent or critical (design-wise) 
ralues of e. 

The author has made some valuable and timely comments under the head- 
ng “Equivalent Pipe Concept.” However, the numerical values of Table I, 
ised to support some of his arguments, pass beyond the realm of design prac- 
icality. In water works design, an overall arterial, feeder, or express main 
ead loss exceeding about twenty feet per mile (equivalent to nearly one 
ound per city block) at peak flow can seldom be tolerated. This criteria is 
‘quivalent to a flow of less than 3 cfs in the 12-inch pipe and less than 0.5 cfs 
n the 6-inch pipe of Table II. Even ignoring the unrealistically excessive 
osses cited, it is implied that the field data from which the pipes would be 
calibrated’ would be sufficiently precise to assign three significant figure 
efficients and exponents to the “exponential equation.” Such precision is 
jarely attained and should not be expected for typical distribution system 
‘- conditions. 

_ Equation (20) represents nothing more than a restatement of the equiva- 
ent pipe concept and has been used for a variety of purposes in this or simi- 
ar form for several decades. The equivalent pipe procedure is occasionally 
mployed in modern network analyses as a computation aid (for instance, to 
educe two pipes in one street to a single resistance for simplification, if 
Bette). Subsequently it will be shown that the procedures proposed under 
ne heading “Exponential Equations and Complex Pipe Systems” have little 

ace in the analysis of distribution networks. Further, the writer submits 

at when a distribution problem is as simple as the “complex” cases 

aguely mentioned by the author, there is no logical excuse for using an em- 
irical fit to the general resistance diagram except possibly for a first trial; 
D be consistent, the diagram itself should be used in determining the final 
mswer. The author states that with a twenty-fold variation in Q (quite hy- 
‘othetical) the exponential equation will agree with the diagram within 5%. 
or his first example in the Appendix, this statement might appear to be cor- 
sct: the head loss for the 0.544 coefficient given would be about 0-2% below 
for the diagram. However, there is an arithmetical error in the example 
d the head loss should be 4.88-feet and the coefficient should be 0.454. A 
enty-fold increase in Q for the corrected equation will yield a head loss 
ut 17% below that for the diagram. The C for the corrected example is 
5, which is not likely to be associated with a small, 12-inch main. 


. 
7 
£ 
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Regardless, a more pertinent comparison would be for a doubled flow for 
which the corrected equation would lead to a loss about 2% below that for the 
diagram while a C of 145 in the Hazen-Williams formula indicates a loss 
about 4% below the diagram. (Since field test data is often acquired during 
periods of near-average demand and since peak hour flows in large arterial 
systems are typically two to three times the annual average, a rate approxi- 
mately double the test rate should represent a reasonable estimate of the 
probable maximum). 


Distribution Network Design 


Network analyses commonly comprise over 90% of waterworks distribution 
design problems. At the present time there are two rather widely used me- 
chanical means for the detailed calculation of head losses in complex gas and 
water distribution networks: by direct analogy using the McIlroy Network 
Analyzer, 24) and by iteration via the Hardy Cross relaxation procedure using 
a digital computer. 25 

Calculations on a digital computer are restricted to a constant m-value 
(1.85 normally used for water service), since a large number of exponential 
equations are being solved simultaneously. Introduction of a varied relaxa- 
tion factor (differing values of m) would hopelessly prolong the iteration 
process. With the McIlroy Analyzer, the m for individual pipes can be ar- 
ranged to fall anywhere between 1.85 and 2. Thus, in an analogue calculation, 
it would not be impossible to assign different values of m to individual pipes. 

However, the purpose of these calculations is normally for the design of — 
facilities to meet peak demands under future conditions (head losses for peak 
demands under existing conditions can be estimated more easily and more 
expeditiously by means of extrapolation of direct field tests). In this connec- 
tion it would perhaps prove fruitful if some of the principal assumptions nor- 
mally implicit in network analyses were reviewed: ‘ 


~ 


a. Neither all the arterial pipes in a system nor all segments of the pipes 
‘calibrated’ are accessable for field head loss measurements. Field head 
loss measurements can seldom be made over a sufficient range of flow, even 
if induced, to cover peak demands. (Design of large arterial networks is 
generally governed by peak-hour conditions). Field head loss measurements 
over a range of significantly different rates of flow are very expensive; 
measured values normally would have to be extrapolated to greater rates for 
future peak conditions. 4 


b. After ‘skeletonizing’ the system to an arterial network, loads (or ‘draw 
offs’) must be somewhat arbitrarily consolidated and arranged as to their 
location. The proportion of the total demand assigned each load is based . 
upon annual averages. (The most reliable information, meter readings and — 


leak-and-waste-survey data, is usually in terms of an average demand rate). 


c. Due to the normal absence of adequate information, it is commonly as- 
sumed that each local load fluctuation has been and will be proportional to 
total system demand. 


d. Future total system demands, future friction coefficients and the relo- 
cation, proportion and variations of loads for the future are all necessarily © 
judgment projections. “ 
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These and additional but less important approximations would thoroughly 
scOurage a designer were it not that ‘control’ analyses of existing networks 
rfiormed under these assumptions usually check out remarkably well with 
id system-head-loss measurements. One aspect of this fortunate coinci- 
nee has been at least partially explained by Mcllroy.(26) He made a study 
the effect of random errors in assumptions of values of C or f for the 
lividual pipelines in networks and found a “pronounced tendency for rather 
vere, randomly distributed, individual errors to cancel themselves in the 
termination of total head losses.” He concluded from his study that “there 
pears to be no strong reason to insist on the use of rigorous expressions in 
solution which is not based on precise data.” And further, “For network- 
alysis work, the widely used Hazen-Williams formula is entirely adequate, 
the same degree of judgment is applied to the results as is necessary in the 
e of — (other) —methods.”(27) 

For the several reasons outlined in the preceding paragraphs of this sec- 
n, the writer is forced to reject as obviously unrealistic (and clearly not 
oven) the statement by the author: “Although the methods for analysis of 
tworks will not be considered here, it is pertinent to note that exponential 
sistance relations will be most convenient for use in this type of analysis.” 


The Hazen-Williams Formula 


Although the Hazen-Williams formula does not deserve very serious con- 
leration in the design of express or feeder mains, some relation with a con- 
int m-value is not only convenient but practically necessary for network 
alysis. This equation can be compared to any point on the general resistance 
gram in the turbulent flow region. The statement “Thus the Hazen-Williams 
‘mula with its exponent of 1.85 corresponds to a combination of relative 
ighness and Reynolds number lying along the line m = 1.85 on the Moody 
gram” is completely incorrect. Fig. D-2, a photo-copy of Fig. 5, p. 258 

the author’s reference (14) clearly shows the comparative relation between 
and f, The equivalent friction coefficient for the Hazen-Williams formula 

1 be expressed as 


| 1060 
| f .oooCo 
H.W. wget p?+015 


: Om. 


water at 60° F. Comparing Fig. 1 and Fig. D-2, it can be seen that the 
= 1.85” curve in Fig. 1 is for a variable C, ranging between about C = 140 
1 C = 150, since the points of intersection of the “m = 1.85” line in Fig. 1 
inot collectively describe a single-valued exponent. Obviously, neither can 
Hazen-Williams formula be described by a single line on the diagram. 


Lamont’s and Other Data 


The author’s Fig. 2 is a redraft of Lamont’s Fig. 3. The table which should 
e accompanied the author’s Fig. 2 is given here as Table D-2. Lamont did 
give the original data. All but twenty of the experiments cited have been 

d previously in the works of Colebrook or others. The five highest curves 

Fig. 2 happen to be for new G. I. pipes with a diameter ranging between 
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Reproduction of Figure 5, p. 258, Reference (14). 
(Note that 2 
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TABLE D-2 


Reproduction of "Table 1l- Summary of Experimental Records for Ne ; Zz 

nt w Pipes", p. 2! 
eference (14) by Lamont (Describes pipes of Lamont's Fig. 3, redrawn by the eke 
is "Fig. 2- Resistance Records on New Pipes"), 


| | Range of experiments Equivalent surface 
| roughness k: i 
Type No. | sere gs ki: inch _ 
of of | Diam: city: Reynolds ee 
pipe expts.| inches ft/sec. Number | Renee Mean rare: 
design 
value 
Uncoated peed 1.175 2.43x104 0.0060 
cast iron 3 to to to to 0.00E9 0,010 
19.68 16.168 7.66x10? (0.0130 
tines et ho | DO. SGI)G) OsbexlOe 0.0007 5. 4 = = ua 
cast iron eee’ to to to to 0.0040 0.005 
61.0 250 6.48x10° 0.008 
Galvanized 0.35C 0.556 3.40xl0" 0.0016 
iron 9 to to to to 0.0040 0.005 
4.110 20.50 515x109 __ 0.0110 
Wrought iron 0.3004 O.1844 &.02x10 0.0003 
18 to to to to 0.0020 0.002 
8.0481 22.50 9.11x10°__ 0.0040 
Coated steel 0.47 0.249 7.57x10 0.0006 
7 to to to to 0.0018 0.002 
IG.12-. 22,05 4.09x107 _ 0.0050 
0.62 0.75 3.75xL02 0.0003 
Uncoated steel 10 to to to to 0.0011 0.0015 
12.0 6.00 1.28x10© 0.0032 
Asbestos £97 0.72 1.90xL04 Smooth Snooth 
cenent 3 to to to pire pipe Smooth 
: 5.77 5.68 1.80x10? 
Spun 3.41 0.53 267X104 Snooth Smooth 
concrete 6 to to to to pipe Smooth 
Spun 10.000 0.048 1.41x10 Smooth Smooth 
bitumen 7 to to to pipe pipe Snooth 
lined 6.25 ah 1.24x10 
' 0.24204 0.145 6.30x1 Smooth 
Smooth pipe” 16 to to to to in- Smooth Smooth 
: 9995 __32.30__8.74x105 __definite pipe ____ 
: 13.10 2.08 2.38x10> Smoot OE, 
Concrete “Ff to to to to Table 
Tate noel 2.825 2.38x102 0.0069 
relined 3 to to to to 0.0130 0.015 


i 
| 
| 
: 
*gnooth pipes include smooth drawn non-ferrous pipes of aluminium, brass, copper 


pipes aha 5.160 8.30x102 0.0195 
Cr 


lead, etc. and non-metallic pipes of Alkathene, glass, Saran, etc. 
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0.35 and 2-inch. These curves, as well as several other irrelevant curves 
for diameters up to about six inches, should have been omitted in the redraft 
by the author; effects arising from differences in manufacture on roughness 
are not considered to be separable for larger sizes, but small sizes, parti- 
cularly under 3-inches, are known to exhibit dissimilar characteristics. In 
terms of distribution design, any size smaller than 6-inches would be of little 
interest. (The curve well below the “smooth” curve is for an 11-inch “spun 
bitumen-lined” pipe). With the curves for the small sizes omitted, it is pos- 
sibly reasonably to argue that “for (selected) new pipes as exponent m from 
1.80 to 1.90 is common - - - - .” In Colebrook’s(5) Fig. 2 (reproduced by 
Hinds(12) as his Fig. 8) which summarizes selected tests on sixteen new 
“centrifugally spun concrete” lined C. I. pipes 4-inch to 60-inch in diameter 
and on six new “centrifugally spun bitumen” lined pipes 11-inch to 46-inch in 
diameter, the divergence from a “smooth” pipe characteristic is slight. The 
same result is indicated in Table D-2 for spun linings (but it appears that 
much of the same data was probably reused by Lamont). 

Although the selected data discussed above is reassuring, it is not neces- 
sarily representative of typical field conditions. In Table D-3 is presented a 
summary of results from field measurements on cleaned existing pipe with 
1ew cement linings obtained under the direction of Mr. Herbert Thamer, 
resident engineer in Philadelphia for the Pitometer Associates. Not many of 


(Dise. of 1962) 


TABLE D-3 
Recent Resistance Tests on Cleaned and Cement Lined Fhilsdelphis wiater Department Mains 
NEW LININGS 
Internal Average Length Head Reynolds Friction Equiv. Hazen- 
Diameter, Velocity, Used in. Loss Number Factor Sand Williems 
Inches, fps, for Test, per 1000!, ¢ 60°F f Rough. , S 
Di Di ok. Ft. x07) Ft. for Di 
47.5 3.47 11,040 0.663 11.0 -0140 -0006 138 
35-3 4.02 2,780 1,251 9.9 0147 -0006 138 
29.8 2.31 1,270 0.475 4.8 0142 -0002 147 
29.4 3.51 10,190 1,047 ise 0134 0002 148 
29.8 2.87 3,690 0.885 5.9 0172 0011 131 
23.6 3437 5,950 1.500 505 -0168 0007 134 
23.6 3.91 3,500 2.01 6.4 -0167 0007 132 
23.6 4.60 6,240 2.43 7.6 0146 0006 140 
23.6 e515 1,770 3.13 8.5 «0150 0004 137 
19.8 2.68 8,565 1.195 3.7 0176 0007 134 
19.5 3.45 7,500 Dein Gh 4.7 0156 0003 141 
19.3 2.19 9,500 0.793 2.9 0171 0005 140 
19.5 2-71 6,470 1.4 a7, 20163 0004 
19.5 3.45 7,460 1.945 4.7 0171 0007 135 
19.3, 2069 ithe} 1.191 3.6 0170 20006 143 
19-4) 4.19 4,365 2.92 Sat 20173 20007 132 
19.4 241 2,710 1.18 3.3 0211 20015 124 
n.4 345 2,400 4.60 2.7 20236 0017 118 
8.0 1.91 940 1.96 1.1 20230 «0006 130 


*Including 6 bends and 4 reducers. 
"Including 2 bends and 6 reducers. 
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these tests represent a near-“smooth” characteristic. There is considerable 
scatter in the values of e despite the possibly consistent ‘statistical’ rough- 
ness of a machine-applied lining. The scatter can be attributed to differences 
in local losses due to alignment of pipes and bends and fittings, assuming that 
the tests were of consistent accuracy. 

The author’s Fig. 3 is a redraft of Lamont’s Fig. 4,(14) Colebrook and 
White(4) in 1937 had observed that “old pipes - - - - obey the square law,” 
and that their graphs included rather flat curves in the transition region 
(same data as used by Lamont?). In preceding comments it was argued that 
the anticipated field roughness characteristics for “new” commercial pipe 
are necessarily inexact, (viz: results to be expected from cleaning and lining). 
To plot “old” pipe friction factors on a general resistance diagram is to sug- 
gest that the increase in roughness with age indigenous to a given water sys- 
tem (i.e. incrustation due to a variety of local causes) is as generally applica- 
ble as the somewhat inexact case of new commercial pipe. It is difficult to 
believe that incrustations are necessarily of the same pattern and form, or 
yield equivalent friction results from one system to another for chemically 
different waters. Details on the data plotted in Fig. 3 are lacking in 
reference (14). The evidence presented by the author does not provide the 
reader with an adequate basis for concluding with any assurance that “for old 
pipes an exponent of nearly 2.0 is common.” A partial basis for this argu- 
ment is given by the author by means of the m-values of his Table II, in 
which the Hazen-Williams C turns out to be only 89 for a 3-year old, and 100 
for a 5-year old cement-lined 24-inch pipe. In both instances it would appear 
that the lining work was substandard, or a number of bends and fittings were 
included, or that tuberculation in Austin is extremely rapid, or that the meas- 
urements were not very accurate, or some combination of these factors. For 
comparison, data from the Pitometeg Company for some “old” cement-lined 
pipes in Philadelphia are presented in Table D-4. It is hoped that the author 


will supply similar details for the pipes listed in Table III. 


CONCLUSIONS 


Continued use of the Hazen-Williams formula in the analysis of water dis- 
tribution networks will persist until such time as the accuracy in appraising 
the several major assumptions affecting design is substantially improved. 
Consideration of detailed refinements for only one factor, pipe resistance, 
may create the illusion that by default the other variables are associated with 
a small margin of error. For estimating head losses for future new pipes or 
linings in a network, the Hazen-Williams formula is of the same order of re- 
liability as the general resistance diagram when all other factors affecting 
the design are considered. The convenience in network analysis afforded by 
a relation with a single exponent (relaxation procedures are made feasible; 
tables and other computational aids are available) outweighs the probable 
greater accuracy of the diagram for existing pipes in good condition which 
have been field calibrated at one or more rates of flow. When “old” pipes 
‘comprise part of an existing network, the expediency of a constant m-value 
for flows in all pipes appears to be justified since arterial mains in poor 
condition must often be cleaned and lined to meet expanding demands for the 
future; if rehabilitation is not indicated it follows that the line will not carry 
‘an important amount of flow and hence it’s contribution to the system head 


' 
B'.] 
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TABLE D-4 
Recent Resistance Tests on Cleaned and Cement Lined Philadelphia Water Department 
Mains. 
"OLD' LININGS 
Internal Age, Average Length Head Reynolds Friction Equiv. Hazen- 
Diameter, Years Velocity, Used in Loss Number Factor, Sand Williams 
Inches, fps, for Test, per 1000!, @ 60°F f Rough., Cc 
Di Di Ft. Ft. x1079 e; for Di 
Ft. 

5he 12 Eyes 3,900 0.813 14.0 0171 0023 122 
5he 10 3.78 6,075 0.755 14.0 0153 0014 129 
46.8 12 4.69 95450 1.447 15.0 0165 0016 124 
47.5 41. 3.75 5265 1.010 12.0 0183 20027 ng 
47.8 nh 3.62 4,750 0.952 12.0 .0186 0029 ns 
4765 10 Dati 6,045 0.895 13.0 -0160 0014 128 
47.8 10 5.62 6,085 1.740 19.0 0142 0008 132 
23.6 12 2-49 10,700 1.275 4el -0261 20055 110 
2326 5 2.58 3,400 0.796 hed 0151 0003 145 
19.4 6 2.70 440 1.530 3.6 .0218 0023 124 
19.3 p) 2.57 4,160 1.280 3.4 «0200 ‘0014 126 
19.3 5 3.41 6,030 2.310 4.6 0205 -0016 122 
15.6 4 2.04 12,060 0.976 2.2 20195 0009 133 
15.3 nA 1.81 6,400 0.980 1.9 20245 20026 ng 
15.4 7 2.20 10,625 1.200 2.4 20205 0012 132 
15.4 7 Auld 5,700 1.192 2.5 20185 0006 136 
15.4 5 2.75 43550 1.570 2.9 -O171 20004 140 
11.4 ales) 2.03 10,350 1.265 156 0188 20004 140 
11.3 13 248 5370 2.250 1.9 0221 20012 126 
11.6 10 3.95 5,415 5.010 3.2 20200 0009 128 
11.6 8 2.07 25225 Ty if 40: 1.7 20257 0021 neg 
11.6 ?) 4035 3,725 5.065 3.5 ~0166 0003 140 


loss will be relatively minor. At the present time there is insufficient evi- 
dence available to warrant adoption of further design complications through 
use of a mixed set of m-values for analysis when a network includes “old” 
pipes. 

In the design of non-network (‘express,’ ‘pumping,’ ‘feeder,’ etc.) the 
empirical curve-fitting technique offered by the author is a poor and unneces- — 
Sary substitute for the general resistance diagram itself. The author dis- | 
dains the use of empirical formulae, yet proposes a system of infinite ap- 
proximation formulae which he purports to be “correct” and “proper” because 
one value for a given pipe is made to coincide with the diagram. The author 
treats only the case where a precise value of e is known in advance, a rather ) 
atypical design situation. : 

It is quite evident that there is much remaining to be investigated from the — 
standpoint of theory and basic experiment if a higher degree of precisionin — 
design calculations is to be achieved. Periodic shuffling of warmed-over data 
is not a very satisfactory substitute for fundamental investigation. Itis sin- _ 
cerely hoped that the ASCE Hydraulics Division “Task Force on Flow in Large 
Conduits” will stimulate interest in further basic research on pipe friction 


characteristics. 
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| 
: NICHOLAS BILONOK. 1_the improved knowledge incorporated in the 
uthors general pipe resistance diagram shows the effect of the various fac- 
rs entering into the problem, which is of great importance to the civil engi- 
.: Valuable experimental data regarding resistance in existing water 
rns are recorded in this paper. 
The general pipe resistance diagram has not received wide usage among 
| hehe engineers due to its inconvenience as compared with similar 
nential formulas for water flow where viscosity has little effect and 
estricted range of Reynolds number, and the velocity of flow, also the dia- 
am in its common form, trial-and-error solution is required for finding 
e prue of flow or pipe diameter; and for given boundary material the pipe 
ameter must be known before the relative roughness may be evaluated. 
The diagram is based on the Prandtl-Karman formula for pipe resistance 
slations which are improved by the experiments of Nikuradse, Colebrook, 


Civ. Engr., Eng. Civ., Corps of Engrs., U.S. Army, Philadelphia, Pa. 
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and White in form (authors Table I) 


At eel BG 
= Sel Cee Net 7 5) (1) 


2: 
Vf 


The friction factor f expressed in Equation (1) is a combination of variables, 
which are evidently applicable to the problem, in very complicated form 
where f is dependable from itself, or friction factor must be known before 
the friction factor may be evaluated. 

It is evident that aging of the pipe under varying conditions of use will re- 
sult a varying of absolute roughness which at present are not easily pre- 
dicted. The experiments show that the surface of the pipe is changing with 
the age increasing the value of f. The factor e expressed in inches is a 
complete measure of the size, shape, and distribution of individual roughness 
element. 

The pipes that are rough in practice generally have no uniformity of the 
roughness element as assumed these by the authors of the equation. As there 
is at present no method for measuring or defining distribution of roughness 
element, the value of e is hypothetical. (1) In author’s examples illustrate 
the steps in obtaining an exponential equation with the aid of the modified di- 
diagram 


Lit 0.155 (1000) (4.5)~_ 


instead hp= 5.84 ft 


h 4.88 
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hp= 0.456 Qt+88 


Recommended use of Fig. 1 to find exponent in example 2 m = 1.95 is not in 
agreement with f and e/D. It would be helpful to show symbols how to find 
numerical variables entering problem in examples. 

The pipe resistance relations listed in the Table I as well as many other 
equations known by the writer, were based on experimental data with the floy 
of water, therefore, in common language, in development of pipe resistance - 
relations, must be used term water or liquid versus fluid. ( (@y 

“A liquid is a body in that state in which the particles move freely among 
themselves, but remain in one mass, keeping the same volume but taking _ 
always the form of the containing vessel; a liquid is an inelastic fluid; a gas 
is an elastic fluid that tends to expend to the utmost limits of the containing 
space. All liquids are fluids, but not all fluids are liquids; air and all the — 
gases are fluids, but they are not liquids under ordinary circumstances, 4 
tho capable of being reduced to a liquid form by special means, as by cold ‘ 
cat Some Water at the ordinary temperature is once a fluid and a : 

q 

The values of e, f, and exponents from the modified diagram in the Darel 

Weisbach equation are important for designers to know these values proper] 


" 
; 
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and apply it intelligently. It is probable that, if accuracy alone were the de- 
ciding factor, any of the values given from the modified diagram, as well as 
C, f, n, and exponents in other suggested formulas which are in general use, 
if properly applied would be equally satisfactory. 
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J. V. DAWSEY, ! M. ASCE.—Professor Moore has collected information 
which was heretofore scattered in many parts of the literature. Research to 
determine methods for predicting or computing resistance encountered by 
fluids flowing through a pipe is of interest in various fields of engineering 
and industry. 

Numerous experiments to determine the relative roughness of a pipe, & D? 


have been completed throughout the world in recent years. Some of these 
tests were made for the purpose of correlating the general resistance equa- 
tions with friction and roughness data collected in model and prototype tests. 
The late Andre Jorissen(1) of Belgium experimented with relative roughness 
measurements on large conduits by forcing air at high velocities through 
concrete-lined and unlined rock tunnels. Viparelli(2 of Italy recently has 
experimented with a photographic determination of roughness. 
The Waterways Experiment Station is using a method devised by Mr. Frank 
B. Campbell,(3) F. ASCE, to determine the equivalent sand roughness of con- 
duit linings. Two-in.-square plaster cast impressions are made at selected 
locations in the conduit that represent the roughness of some percentage of 
the total conduit lining. Vertical-displacement measurements of the cast sur- 
face profile yield an average surface roughness of the conduit lining repre- 
sented by the cast. 

An Ames dial is presently used for surface-profile measurements on the 
plaster casts after aninitial attempt to measure profiles electrically. The 
point of contact on the shaft of an electrical contour analyzer was modified 
by using the point from a ball-point pen. This pen point indented the surface 
of the soft plaster casts under force of the analyzer probe. The Ames dial 
shaft was also modified by the addition of a similar pen point, but did not in- 
dent the cast surface during measurements. For profile measurements, the 
casts are firmly anchored in a mechanical vise and traversed under the shaft 
of the Ames dial. If the root-mean-square deviation of the surface -profile 
measurements is interpreted as the equivalent sand-grain roughness € (one- 
half the diameter of the sand grain k in inches), then the friction factor f, at 
the Reynolds number of the test, can be computed by using the general pipe 
resistance equations. 

Roughness casts of representative areas inside the 20-ft-diameter 


it. Hydr. Engr. (Research), Prototype Section, Hydr. Analysis Branch, 
Hydraulics Div., U. S. Army Engr., Waterways Experiment Station, 
Vicksburg, Miss. 
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flood-control conduit No. 6, Denison Dam,(4) were obtained by making two- 
inch-square impressions of the conduit lining with elastic dental-impression 
material. Plaster casts were then reproduced from the impression material 
negative. Equivalent roughness € for the conduit lining was determined from 
root-mean-square deviations of the surface profile measurements on the 
casts. This value of € was then used in the Colebrook-White equation, which 
includes the effect of Reynolds number, to compute the Darcy-Weisbach fric- 
tion factor, f. The value of f computed in this manner was about 35 per cent 
greater than the value of f computed from pressure measurements with piezo- 
meters installed in the conduit during construction. Similar prototype pres- 
sure measurements and roughness casts were recently completed at Enid 
Dam.(5) The value of f determined from the casts was only 15 per cent less 
than the value computed from the slope of the friction gradient. 

Although there was considerable difference between the measured and 
computed values of roughness, €, the corresponding friction factors differ 
much less because of the exponential-logarithmic relationship between f and 
€. With some additional refinment in sampling the conduit surface and 
analyzing the plaster casts, it should be possible to determine conduit rough- 
ness accurately from plaster casts alone, including conduits where different 
types of forms and construction techniques are used. 

Additional plaster-cast roughness data are needed to aid in the analysis 
and interpretation of the equivalent sand-grain roughness relationship in 
turbulent flow. Information from anyone making casts in conduits where the 
hydraulic friction gradient is known or can be measured would be appreciated. 

It appears that parentheses were inadvertently omitted in the Colebrook- 
White equation in Table I of Professor Moore’s paper. The equation should 
read: 


ein aeng’s Hogs fen + 
Zz N, Vi 3.2 
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HYDRAULIC ANALYSIS OF SURGE TANKS BY DIGITAL COMPUTER 2 


Discussion by Charles Jaeger 


CHARLES JAEGER!.—The paper of Mr. Barbarossa deals with two most 
mportant problems. 


[he use of digital computers 


The use of digital computers for advanced research on surges in complex 
systems of surge tanks is clearly demonstrated by the author. Digital com- 
uuiters have been used in hydraulic research by several authors especially 
or solving problems of water hammer. 

The paper by McCaig and Jonker, is a typical example of what can be 
ichieved with this new tool. The writer too has used a digital computer for 
solving a simple but tedious problem. It was suggested that a small hydro 
ower station could be provided with a surge tank too narrow to cope with 
arge surges, provided the power sets were hand operated and put slowly on 
he line, the operator observing very definite rules concerning the turbine 
‘ates opening times. There was no turbine governor and the load on the sta- 
ion was supposed to remain constant for long periods. 

When a large number of possible load combinations and possible errors by 
he operator were considered and their effect on surges analysed, the pro- 
osal was abandoned. 

The approach by digital computer can be compared with more conventional 
‘raphic methods: For example the paper by McCaig and Jonker, can be com- 
ared with graphical computations carried out by Bouvard and Molbert of 
“lectricité de France. The precision given by a direct calculation on the 
ligital computer goes far beyond of what can be achieved with a graphic 
iagram. In the case described by McCaig and Yonker such high precision 
yas required for comparison of cost estimates for several design alternatives. 
Paynter has been using the electric analog method for research on surges 
don water hammer. His methods are highly suitable for very general re- 
earch. 
| When problems concerning concentrated pressure and head losses (losses 
iccurring at this junction of a surge tank with the pressure tunnel) have to be 
malysed or when turbulence is involved, experiments on model tests is often 
ne best approach. 


1e Thoma condition for surge stability 


“The author correctly remarks that the whole surge stability theory estab- 
ished during the past 40 years by Thoma, Calame and Gaden, Escande and 


. Proc. Paper 1996, April, 1959, by Mr. Nicholas Barbarossa. 
. English Electric Co., Ltd., Rugby, England. | 
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others is based on the assumption that the demand for power varies abruptly, 
then stays absolutely constant for the period under observation and that the 
governor is able to correctly adjust the power generation to the demand. 

The paper shows how trial calculations on the digital computer prove this 
basic assumption to be over-simplifying. According to Fig. 10, 11, and 12 of 
the paper, the gates continue to open to “full gate opening” beyond the point 
where better stability could be expected. * 

It should be remarked here that the points on Fig. 10 and 11 on which the 
discussion is based correspond to the very beginning of the surge: time abou 
1 second, whereas the surge period is about 200 seconds (Fig. 7.). There is 
no precise indications given in the paper why a variable surge with constant 
power generation could not be established during the rest of the period. De- 
pressed water levels will raise again. In addition, the governing equation re- 
lating speed variations to the moment of inertia of the rotating masses could 
possibly be introduced as an additional condition to the computation: it may 
show how the governor tends to stabilise the power output and additional re- 
search is suggested. The main point to be stressed here is the authors be- 
lief shared by the writer, that further research on surge stability is definitel 
wanted. A first phase of research has been concluded: simplified mathemati 
cal theories expanding Thoma’s stability creterion to all types of surge tanks 
and systems of surge tanks have been worked out. The Routh-Hurwitz mathe 
matical stability criterion and the Liapunoff theorem have been used success 
fully by skilled mathematicians. All this approach work is based on the sams 
basic assumption of abrupt change in the load demand, which then remains 
constant. : 

How does this basic assumption fit in with the behaviour of real electric 
systems? This is the question facing us for the second phase of the re- 
search. A certain amount of preparatory work has already been done: 


Chevalier and Hug (Electricité de France) have tested the Cordeac surge 
tank; They found that the governor was reacting correctly to the power de- 
mand but with a time lag of about 10 seconds. This time lag was found to 
dampen the surges quite efficiently. 


Cuenod and Gardel (Lausanne) have tested the differential surge tank of 
Oelberg (Switzerland). By modifying the actual surge tank they caused it to 
become unstable. They tested different types of feed back linking surge leve 
and line voltage to the frequency in order to re-establish the stability of tne 
surges. 

Finally Aillert (Electricité de France) has designed and built several 
surge tanks, which are definitely unstable but where the surges are con- 
trolled otherwise. ; 

These are definite efforts to approach the stability problem from a com- 
pletely new angle. More light could be thrown on the way governors react 
and surges can be controlled. Possible savings on the tank volume could be 
achieved thanks to a closer analysis of the physical conditions in large elec: 
tric systems. 

It is suggested that a comprehensive paper summarising all these new ef 


forts should analyse the up to date position of the problem and show the way 
to further research. 


: 


*Similar effects have been observed on large tanks. Too little is known for 
any explanation to be suggested. d 
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The author has to be warmly congratulated for his successful use of a new 
tool in surge computation and for his criticism of the assumptions on which 
the present surge stability theory is based. 
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TWO METHODS TO COMPUTE WATER SURFACE PROFILES# 
Discussion by R. Silvester 


R. SILVESTER!.—In their use of Manning’s formula the authors have 
julked area, hydraulic radius and the roughness coefficient “n” into the con- 
reyance factor Kg. There seems little to gain from including Manning’s 
efficient in order to determine some nebulous quantity of conveyance which 
s of little use in itself. The writer discarded this concept some time ago(1,2) 
vith the result that K(=1.486 A R 2/3) is determined solely by the geometry 
f the channel cross-section and the stage. Such values when graphed on a 
ogarithmic scale resulted in straight lines and slight curves depending upon 
he shape of the channel. Graphs were supplied for trapezoidal, parabolic, 
‘ectangular, circular and U shaped channels. Once K has been so determined 


the discharge is given by Q = Ls VSs 3; for the solution of this a nomograph 
Was included. = 
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